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I. Teil: Bestimmung von strengen Lösungen 
der Einsteinschen Feldgleichungen 


Einleitung 


Durch die Entdeckung der allgemeinen Relativitätstheorie 
war es möglich, die Frage nach der Beschaffenheit des Raumes 
im Großen von der physikalischen Seite anzufassen, und sie 
mit den Massen im Weltall in Verbindung zu bringen. An- 
fänglich schrieb man infolge nur kleiner beobachteter Stern- 
geschwindigkeiten dem Raum statische Struktur zu. Dies 
führte zu den statischen Modellen?) des Weltalls. Die Mes- 
sungen Hubbles*) an den außergalaktischen Nebeln, die eine 
systematische, zum Abstand vom Sonnensystem proportionale 
Rotverschiebung, und Radialgeschwindigkeiten bis zu 20000 
km/sec ergaben, wiesen der Theorie neue Wege. Die Annahme 
einer nichtstatischen Metrik im Weltall wird diesem Tatbestand 
gerecht. Schon 1922 hat A. Friedmann‘), ohne Kenntnis 
von Beobachtungsmaterial eine nichtstatische Lösung der Feld- 


1) Auszug aus einer Berliner Dissertation. 

2) A. Einstein, Berl. Ber. 1917, S. 142; W. De Sitter, Monthly 
Not. R. A. §. 1916, 1917; M. v. Laue, Berl. Ber. 1923, S, 27. 

3) E. Hubble,. Proc. Nat. Ac. 15. S. 168. 1929; Astrophys. Journ. 
4. S. 43. 1931. 

4) A. Friedmann, Ztschr. f. Phys. 10. S. 377. 1922. 
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gleichungen angegeben. Die den Friedmannschen Unter- 
suchungen zugrunde liegenden physikalischen und metrischen 
Voraussetzungen wurden später von verschiedenen Forschern 
erweitert, und diese Theorie zur Deutung des Hubbleschen 
Phänomens herangezogen.') Die Friedmann-Lemaitresche 
Lösung weist als Streckenraum einen Kugelraum mit zeitlich 
veränderlichem Radius und homogener Massenverteilung auf. 
Es gilt die Maßbestimmung’) 


(1) ds*= P?.S\y,,da* -dz4—d = 


P bedeutet den Radius, Sy,,-da*-da* das Linienelement 
eines dreidimensionalen Kugelraumes vom Radius 1. Zwei im 
räumlichen Koordinatensystem von (1) ruhende Punkte besitzen 
Relativgeschwindigkeiten 
1 

gegeneinander, wenn r der Abstand der Punkte zur gleichen 
Zeit ist. Wir wissen, daß man in erster Näherung das 
Dopplersche Prinzip zur Berechnung der Rotverschiebung 
anwenden?) darf. Um in vollem Einklange mit den Hubble- 
schen Beobachtungen zu sein, hat man nur anzunehmen, daß 
zur Zeit D>0. 

Die mit dem kosmologischen Gliede versehenen Feld- 
gleichungen ergeben nur eine statische Lösung zum kosmologi- 
schen Problem, die einer mit Materie erfüllten Welt entspricht, 
die Einsteinsche Zylinderlésung. Es hatte den Anschein, 
als ob die allgemeine Relativitätstheorie die Frage nach der 
Beschaffenheit des Raumes im Großen eindeutig beantworten 
könne. Das änderte sich, als man der Welt nichtstatischen 
Charakter zuschrieb. Es zeigte sich nämlich, daß man, bei 
Voraussetzung homogener Massenverteilung und sphärischer 
Welt, eindeutige Aussagen über die Art der zeitlichen Ver- 
änderung der Raummetrik nur dann machen kann, wenn man 
das kosmologische Glied unterdrückt. Die folgenden Unter- 
suchungen und inzwischen von O. Heckmann’) veröffentlichte 
Arbeiten weisen nach, daß die Feldgleichungen in genau der- 


1) G. Lemaitre, Ann. Soc. Sci. Bruxelles 47, Serie A, S. 49. 
1927; Monthly Not. R. A. S. 91. S. 483, 490. 1931; W. De Sitter, Proc. 
Nat. "Ae. 16. S. 474. 1930; Bull. Astr. Instit. of the Netherlands 6. S. 41. 
1930; R. C. Tolman, Proe. Nat. Ac. 16. S. 320, 409, 511, 582. 1930; 
Phys. Rev. 38, S. 1758. 1931: Sir A. Eddington, Monthly Not. R. A. S. 
90. S. 668. 1930. 

2) Griechische Indizes laufen immer von 1—3. 

3) M. v. Laue, Berl. Ber. 1931, S. 213. 

4) O. Heckmann, Göttinger Nachrichten 1931, S. 126; 1932, S. 97. 
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selben Weise wie mit sphärischen, auch mit hyperbolischen 
und euklidischen Räumen verträglich sind. Weiter werden 
Lösungen der Feldgleichungen angegeben, die wohl homogene 


Massenverteilungen, aber keine Geometrie konstanter (räumlich) . 


Krümmung im Streckenraum, und solche, die weder räumlich 
konstante Massendichte, noch eine Geometrie konstanter Krüm- 
mung im dreidimensionalen Raume aufweisen. 


§ 1. Vorläufige Problemstellung 

Die Vorgänger haben P als Funktion von x* allein auf- 
gefaßt. Diese Annahme soll hier abgestreift werden. Dem- 
zufolge soll also in (1) P Funktion aller vier Koordinaten sein 
dürfen. Weiter verlangen wir, daß S\y,,-da*-da* den definiten 
Maßtensor irgendeines statischen dreidimensionalen Raumes 
darstellt. 

Es werde angenommen 
3) Ti=-—d-u, 
wo u die Massendichte bedeutet. Damit wird auch 

T,,= 0, T= T = p. 

Physikalisch setzt das voraus, daß 1. man von allen im Welt- 
raum vorhandenen elektromagnetischen Feldern, und dem durch 
sie bedingten Druck absieht, 2. die Materie ruht, d. h. daß 
jeder Massenpunkt seine Werte von 2’, x”, x? dauernd beibehält. 
Dies steht nicht im Widerspruch mit dem Linienelement (1), 
denn die zweiten Differentialquotienten von z', x?, x? nach der 
Zeit verschwinden nach den Gleichungen für die geodätischen 
Linien, wenn dies für die ersten gilt, weil in ihnen kein Glied 
mit (#%)” auftritt. Geometrisch besagt das Linienelement (1), 
daß ein unendlich kleiner Koordinatenbezirk des Strecken- 
raumes nur konformen zeitlichen Änderungen unterworfen ist. 
Der Streckenraum bleibt sich im Kleinen ähnlich. Ist P nur 
von x* abhängig, besteht Ähnlichkeit auch im Großen. 

Gleichzeitig ist es durch den allgemeinen Ansatz (1) 
möglich, zu einer allgemeineren, als der zu speziell scheinenden 
gleichférmigen Massenverteilung, zu gelangen. 

Die Einführung des kosmologischen A-Termes war not- 
wendig für Einsteins statisches Modell des Weltalls. Ein- 
stein hat neuerdings darauf hingewiesen, daß er nicht nötig 
ist für nichtstatische Modelle, und setzte A gleich Null.) Im 
folgenden wird der Vollständigkeit halber das kosmologische 
Glied beibehalten. 


1) A. Einstein, Berl. Ber. 1931, 8. 236; A. Einstein u. W. De 
Sitter, Proe. Nat. Ac. 18. S. 213. 1932, 
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§ 2. Aufstellung der Feldgleichungen 
P und die y,, in (1) sind so zu bestimmen, daß das Feld- 
gleichungssystem von Einstein +4-9,,= —x-(T,,—1-T-g, ) 
erfüllt ist (« Gravitationskonstante, T = 3 T;'). Di Christoffel- 
symbole von (1) sind (die gesternten Größen beziehen sich auf 
daR): 


v y * 1 oP A oP 1 oP 


Die Berechnung der i4-Komponente des Krümmungs- 
tensors R,, liefert 


4 
(4a) 
und 

2 
(4b) 
Ox’ 

Der Materietensor ist divergenzfrei. Dies liefert = 0, 
daraus 
(5) 


Die Feldgleichungen mit den Indizes k = 1, 2,3; 1=4 ergeben 
log P 


= 0. Demnach ist allgemein 


dr 

(6) P= pie, 2’, 2°)» (zr). 
Die zu (4a) gehörige Feldgleichung wird 


2Ps 
Nach (6) wird u von den Ortskoordinaten unabhängig. Man 
wird also notwendig zu einer räumlich konstanten Massendichte 
geführt. In (1) kann nun der Faktor p von P zu den y,, 
hineingezogen, und P als nur von der Zeit abhängig angesehen 
werden. In (5) ist & eine Konstante. Einmalige Integration 
von (7) ergibt: 
dP\ i-P§—3-K-P+xca 

(8) (ga) = 
(K ist Integrationskonstante). 

Zur Bestimmung der y,, dienen die restlichen sechs Feld- 
gleichungen. Es ist nun 


dix* ( dar" 
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alle übrigen sind Null. Sind Ri, die Komponenten des 
Krümmungstensors im „y-Raume“, so ergibt die Rechnung für 
den Krümmungstensor der Raum-Zeit Mannigfaltigkeit 

» dP \? aP 


Die entsprechende Feldgleichung lautet: 


Die Werte für R,, aus obiger Gleichung, und für d?P/dx*” 
und (dP/dx* aus (7) baw. (8) eingesetzt gibt R2,=—2K-y,,). 
Bekanntlich ist im dreidimensionalen Raum durch den Rie- 
mannschen Krümmungstensor R,, auch der Riemann- 
Christoffelsche Tensor R,,,, festgelegt. Um ihn wirklich zu 
berechnen, gehen wir aus von einem Satz von Schouten}), 
der besagt, daß fiir jeden dreidimensionalen Riemannschen 
Raum die sogenannte Konformkrümmungsgröße C,,,  identisch 
verschwindet. Ist 
so wird 
Dieser Satz, auf unseren en one angewandt, gibt 
=K- (Yea Vip): 
Damit ist aber gezeigt, daß unser „y-Raum“ nur ein Raum 
konstanter Krümmung K sein kann. Die Maßbestimmung in 
einem solchen läßt sich immer auf die Riemannsche kano- 


nische Form: 
da” + + da® 


+ + a? + 


do?’= 


bringen. Als Lösung der Feldgleichungen ist also gefunden: 
ds? P2. + + dr 
” | Ki 3? F 
Ist K > 0, aber +1, so ergibt die Transformation 
das Linienelement: 


ue + dat da? 
— (a * + af? + 


1) J. A. Schouten, „Ricei-Kalkül“ Aufl. 1924, S. 170. 
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Ist K <0, aber £— 1, so folgt durch die Transformation 
die Form: 
dat * + da®’+ da” 


_ 
- + a+ af 


d= 


2 

Die Faktoren 
mit P? bezeichnet, ebenso lassen wir an den Koordinaten die 
Striche weg. Außer K=0, kann K noch Null sein. Un- 
beschadet der Allgemeinheit können wir damit K auf die 
Werte + 1,0, — 1 beschränken. 

1.K=1. Es folgt die Friedmann-Lemaitresche 
Lösung. Durch die Transformation: 


die hier auftreten, werden in Zukunft 


2.tg -cos?, 2.189 -sind-cosg, 
= 2.tg-2 ‘sind. sing, 
erhält man: 
(9a) ds?= P?. (dO? + sin?@ -(d9? + sin? — da” 
und P hat der Differentialgleichung: 


dP \?2 A-P§—3-P+2Cu 

(9b) en 3-P 
zu geniigen. 

2.K=0. Es folgt: 
(10a) ds? = P?. (da? + dy? + dz’) — dx* 

y 
und 
dP \2 A> Pe+xCa 

(10 b) 


3. K =—1. Die Transformation, die völlig analog zu 
der unter 1. ausgeführten lautet, wenn tg 2. durch Igo 
ersetzt wird, liefert: 

(11a) ds®= P?. (d0?+ Sin? O - (d9?+ sin? dg’) — 
Ks ist: 
(11b) (<3 


dit 


3-P 


Konformität der zeitlichen Veränderungen des Strecken- 
raumes führt zu homogener Massenverteilung und zu einer 
Geometrie konstanter Krümmung, d. h. zur sphärischen, eukli- 
dischen oder hyperbolischen Geometrie im Weltraum. Neben 
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den nichtstatischen Kugelraum als Modell des Weltalls tritt 
völlig gleichberechtigt der nichtstatische hyperbolische bzw. 
euklidische Raum. 


$3. Diskussion der Lösungen 
Die Zeitabhängigkeit von P in (9a) ist schon von Fried- 
mann diskutiert. In (10) und (11) kann P nur dann peri- 
odisch sein, wenn A<0; wenn A>0 ist P eine monoton 
wachsende bzw. fallende Funktion der Zeit. Sie wächst bzw. 
fällt rascher als eV4= bzw. e—V4=*. Für 2=0 erhält man 
durch Integration von (10): 


3 2 
P= (5 const) 
und von (11b): 
Sin|z -VP@+P)- + const| (#- 
In (9), (10), (11) ist die Massendichte durch „= = gegeben. 
Kiner leeren Welt entspricht «= 0. Dafür folgt speziell: 


1. Aus (9): 
dst = .Cof? 4 ot. (40°+ sin?0(49°+ sin*-dg?) — dat”. 
2. Aus (10): 


i 
ds? = (dx? + dy? + dz’) — dx*, 
3. Aus (11): 
ds =>. Ginty/4 2%. (d0°4 Sin? 4 


Diese drei Linienelemente lassen sich sämtlich durch Ko- 
ordinatentransformation aus dem De Sitterschen gewinnen.’) 

Im sphärischen, euklidischen und im hyperbolischen Falle 
läßt sich das Hubblesche Phänomen als Folge der Expansion 
des Raumes erklären, und es scheint, als ob nur Gefühlsgründe 
in der Bevorzugung der einen oder der anderen Lösung mab- 
gebend sein könnten. Ganz so ist es nicht, wie wir später 
sehen werden. 


1) Wie dies bei 1. ~~ vgl. C. Lanezos, Phys. Ztschr. 23. 
S. 530. 1922. Bei 2. H. P. Robertson, Phil. Mag. 5. S. 835. 1928: 


H. Weyl, Phil. Mag. 9. S. 936. 1930. Bei 3. analog zu 1., nur daß die 
trigonometrischen Funktionen von © durch Hyperbelfunktionen zu er- 
setzen sind. 
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Die euklidische und hyperbolische Metrik führen bei Aus. 
schluß von singulären Punkten, und Einbettung in einen höher- 
dimensionalen euklidischen Raum zur Unendlichkeit des Raumes, 
Besitzen zwei im räumlichen Koordinatensystem ruhende Massen- 


punkte einen Abstand voneinander, der größer als r = + ist, 


so übersteigt nach 2. die Relativgeschwindigkeit der Massen- 
punkte die Lichtgeschwindigkeit. Steht diese Tatsache im 
Widerspruch zur allgemeinen Relativitätstheorie, und wird 
dadurch die kosmologische Anwendungsfähigkeit der Lösungen 
(10) und (11) beeinträchtigt? Das Kausalitätsprinzip verlangt, 
daß das Weltlinienelement der Massenpunkte zeitartig ist. Da 
die Massenpunkte im räumlichen Koordinatensystem ruhen, ist 
diese Forderung unabhängig von ihrem gegenseitigen Abstand, 
und dem zugrunde gelegten Koordinatensystem erfüllt. Vom 
Standpunkt der allgemeinen Relativitätstheorie aus ist somit 
nichts gegen die oben erwähnten Überlichtgeschwindigkeiten 
einzuwenden. 

Zwecks Entscheidung der Frage nach dem Vorzeichen der 
Krümmung des Raumes, dividieren wir (8) durch P? und er- 
halten unter Berücksichtigung der Beziehung (5) 

xo K 
D? = + 3 — Dr . 

Da P reell ist P* positiv. Aus der Gültigkeit dieser Gleichung 
in der Jetztzeit, die durch ein angehängtes b charakterisiert 
sei, folgt dann: K=Z0, je nachdem 4+ xc?u,—3-D,?=0. 
Setzen wir 4=0, und verwenden den Wert 0,5-10-27 cm 
für D,, so kommt KZ0, je nachdem u, 0,4- 10727 g/cm’. 
Durch eine Dichtemessung im Weltall wäre es also prinzipiell 
möglich, wenn wir die Annahme 4 = 0 als bündig ansehen, 
das Vorzeichen der Krümmung des Raumes festzulegen. 

Das kosmologische Glied, das in die Differentialgleichungen 
für P eingeht, steht in keinem direkten Zusammenhang mit 
der Massendichte, und läßt sich auch weder nach Vorzeichen 
noch Größe durch die Massen im Weltall festlegen. So liegt 
die Annahme 2=0 nahe, und Einstein hat die sich er- 
gebende Lösung (9) und auch (10) zu kosmologischen Betrach- 
tungen herangezogen. Letztere verdient wegen ihrer Einfachheit, 
da sie ohne Annahme einer Raumkrümmung die Hubbleschen 
Rotverschiebungen und die Anwesenheit von Materie zu er- 
klären vermag, besonderes Interesse. 

Aus (11b) folgt für A = 0: 

(42)’- 
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=" a ist positiv. In welcher Phase der Veränderung von P 
wir uns auch befinden mögen, immer ist = >1. Damit D-P>1 


oder P> - Zahlenmäßig P, > 2 - 10?" cm. 


Nach Einstein liegt die größte Schwierigkeit für die 
Anwendung der Lösung (9) mit 4 =0 in der Kosmologie darin, 
daB der Zeitpunkt, fiir welchen P gleich Null war, nur etwa 
10° Jahre zuriickliegt. Wie man aus den Differentialgleichungen 
leicht ersieht, ist diese Schwierigkeit auch bei den Lösungen 
(10) und (11) mit 4 = 0 nicht beseitigt, und man entgeht der 
Schwierigkeit dadurch, daß man darauf hinweist, daß die In- 
homogenität der wirklichen Massenverteilung, die gerade bei 
kleinem P, also großer Dichte, besonders in Erscheinung 
treten wird, unsere Behandlung wertlos macht. 


$4. Erweiterte Problemstellung 


Nachdem nichtstatische Linienelemente zu kosmologischen 
Betrachtungen herangezogen wurden, wäre es von höchstem 
Interesse, eine gewisse Übersicht über die mit den physikali- 
schen Voraussetzungen des § 1 verträglichen Lösungen der 
Feldgleichungen zu bekommen. Allgemein wird das, der großen 
mathematischen Schwierigkeiten wegen, nicht gelingen. Physi- 
kalisches Interesse besitzen zunächst, infolge unserer noch 
mangelhaften Kenntnis der Dichteverteilung im Weltraum, 
nur solche Modelle, die gewisse Symmetrieeigenschaften der 
Raummetrik, oder homogene Materieverteilung aufweisen. 
Friedmann, Lemaitre, De Sitter verlangen in ihren Ar- 
beiten über nichtstatische Linienelemente höchste Symmetrie, 
d.i. sphärische Struktur des Streckenraumes. Tolman setzt 
zentrische Symmetrie voraus. Hier soll die Fragestellung 
dahin verallgemeinert werden, daß nicht der Streckenraum, 
sondern nur eine Schar von zweidimensionalen Koordinaten- 
hyperflächen, etwa z*= const nicht sphärische, sondern all- 
gemeiner eine Geometrie konstanter Krümmung haben. Unsere 
Linienelemente sollen von der Art sein: 


(12a) ds? = 4%. (d + sin? .d499) + BY. da?’ — 
oder 
(12b) ds? = A?- (da? + dy?) + B?. da?” — 
oder 
(120) ds? = 4?.(d0?+ Gin?@.d 9) + B®. da?” — da”. 
A sei Funktion von x* allein und B Funktion von ©, 
%, x*. Im allgemeinen wird ein dreidimensionaler Koordi- 
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natenbezirk im Gegensatz zu den in § 2 untersuchten Fällen 
im Laufe der Zeit verzerrt. Da die Rechnungen für die Linien- 
elemente (12) gleichartig verlaufen, möge sie nur für (12a) 
ausgeführt werden. Die g-Klammern von (12a) sind: 


= — sin - cos @; 
-4-44; . -sin20; rr =B. 
Die Divergenzfreiheit des Tensors T,, liefert analog zu (5) 
(13) p= 


Das Nichtvorhandensein von Materiestrémung driickt sich 
in den Gleichungen R,,=0 aus. Dies ergibt (i =1, 2, 3): 
1 dA 22 
(14a) 


1 dA OB 
Wenn Fr) und ao £9 lassen sich diese Beziehungen 


auch so schreiben: 


6 (6B 6 
az (36) aa ng 
oB A da* OB A dz 
OF 
Die Gleichungen sind erfüllt, wenn 
0°B 0° B 


~ 


II. B= w(O, 9.4. 

Im folgenden werden die drei Fälle einzeln untersucht. 

Fall I. Es ist B = Funktion (9, #) + Funktion (x*). Bei 
Zugrundelegung von (12a) existiert eine Lösung nur, wenn 
24> 0, entsprechend bei (12b) und (12c) nur, wenn A = 0 baw. 
2<0. Wir werden im Falle des Linienelementes (12a) 
A =}, im Falle (12b) 4 = 1, und im Falle (12) 4 =; 
wählen. Die Feldgleichung mit den Indizes 44 lautet: 
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Indessen liefert die direkte Berechnung R,, = 5 = . Damit 


und mit (13) folgt 2 —z =4-B- =< “. Das allgemeine Inte- . 
gral lautet: 

2 
(15) B= Vi + 4 


(3,- und £,-Konst). Die Feldgleichungen 0 und R,=0 
sind identisch erfüllt; dagegen kommt aus R,, = 0 die Be- 
ziehung: 


6B 
(16a) 86553 — = 0. 
Die Feldgleichungen 11, 22, 33 unter Verwendung von (13): 
®B __ da 
067 

=. 
(16¢) B sing. cos - sin? 
und ‘ 
x a 


Durch (16b) wird (16c) mit (16d) u Die Integration 
der G1.(16) ergibt = 7 - sin sin #(« Konst), oder@=y-cos 0. 
Somit 

oder 


B= + B,-e-Viw + . 


Durch geeignete Maßstabsänderung bei der Koordinate x? 
können wir die gefundenen Linienelemente wie folgt schreiben 


(7 $ 0): 
ds? = | dO? + sin?9.d492 +(ß,- e-Vü 
+ sin sin | — da", 
ds? = . EG + .d (8, + e-Viz' 
| + cos | —dz*. 
Die Massenverteilung ist gegeben durch: 


(17a) 


@- sind 2h 
und + e Vie, + sin sin nied 
(18b) cos 9 24 


8, Ba e Vit 


(5) 
ich 
3): | 
| | 
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Die Linienelemente (17) und (18) können durch eine räum- 
liche Koordinatentransformation ineinander übergeführt werden, 
Die Identität beider ist damit gezeigt, und wir beschränken 
uns im weiteren auf die Betrachtung des Linienelementes (17), 
Die analoge Rechnung führt, wenn 4 = 0, zu: 


(19a) d?= dr? + ay+[8- = 
(19b) u= 


xe re (8-Konst.) 
(«+ y?+a') 
und wenn 4 < 0 zu: 
ds?= . ja 62+ (Bsin(y—A 29) 
(20 a) 
+ Sin - sin #)?- | 


Gin O sin $ 2} 
Die Linienelemente (17), (19), (20) sind strenge Lösungen 
der Einsteinschen Feldgleichungen. 
Fallll. Bist nur eine Funktion der Zeit. Lösungen dieser Art 
sind mit den Feldgleichungen verträglich, wenn A und B und die 
Dichte, die notwendig räumlich konstantseinmuß, den Gleichungen 


K 1 d? A 1 dA\? 1 dA dB 
508, 
1 #B 2 x 
(#10) 


genügen, wobei K eine Konstante ist, die im Falle (12a) + 1, 
im Falle (12b) Null, und im Falle (12c) — 1 ist. Die erste 
resultiert aus den Feldgleichungen 11 und 22, die zweite aus 
der Feldgleichung 33, die dritte aus derjenigen mit den 
Indizes 44. Diese drei Differentialgleichungen sind nicht un- 
abhängig voneinander. Es besteht vielmehr die Relation (13) 
zwischen ihnen, was sich aus den Gleichungen heraus zeigen läßt. 
(21a) mit 2, (21 b) mit —1 multipliziert und zu (21) addiert, gibt 
K 2 @®A 1 /dA\? 
—A+ =0. 

Durch einmalige Integration folgt 

AA3-3K-A+ß 


dx 


(8-Konst.). 


3A 
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Für 3>0 und K=1 ist dies genau dieselbe Differential- 
gleichung wie (9b), der der Radius P der sich expandierenden 
Kugelwelt genügt. Für K = 0 und K =—1 wird sie formal 
mit (10b) und (11b) identisch. 
Allgemein ist B eine ziemlich verwickelte Funktion der 
Zeit. Wenn 2 = 0, lassen die Differentialgleichungen (21) eine 
Integration zu, die einen expliziten Ausdruck für B liefert. 
Aus (21b) folgt: 
A}. = = -a'+y (y-Konst.). 
Damit aus (21a): 


1A 


Bei geeigneter Wahl des Anfangspunktes der Zeit 
B dA Ixc’« 

und nach (13): 


@ 1 

K=1. A wächst von Null gegen den Maximalwert 8/3, 
um dann wieder umzukehren. Während dieses Vorganges 
schlägt das Vorzeichen von mw unbedingt um, denn für den 
44 =0, also u>0, dagegen ist — 
groß, wenn A in der Gegend des Nullpunktes ist, damit u< 0. 
Dies widerstreitet den Grundvoraussetzungen, die wir bezüg- 
lich der Massendichte machen müssen. 

K=-—1. A wächst monoton gegen Unendlich. Eben- 
falls B. Die Massendichte ist eine abnehmende Funktion 
der Zeit. 

K=0. Man kommt auf (10) zurück. 

Fall III. Dabei ist also B= A- w(9, 9). Das Linieu- 
element läßt sich so schreiben: ds? = A?-(d 0? + sin?O-d # 
+ — Dieses Linienelement ist ein Sonderfall 
des allgemein in § 2 untersuchten Linienelementes, und fiihrt 
zu den dort gefundenen Lösungen. 


sehr 


Maximalwert ist 


§ 5. Betrachtungen zu den erhaltenen Lösungen 


Die Maßbestimmungen (17), (19), (20) entsprechen Mo- 
dellen des Weltalls mit ungleichférmigen Massenverteilungen, 
die sich zeitlich ändern. 

Sind in (17) speziell %, und 8, Null, so erhalten wir die 
Einsteinsche statische Zylinderlösung mit der homogenen 
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Haben #, und A, verschiedenes Vor- 


zeichen, so ist in einem bestimmten Augenblick, der sich aus 
der Gleichung ß, - eViw 4 B,-e-V4* = 0 bestimmt, die stati- 
sche Lösung als Zwischenstadium enthalten. Grenzfälle hier- 
von sind diejenigen Lösungen, bei denen entweder #, oder ß 
Null sind. Im ersten Fall wäre der statische Zustand das 
erst nach unendlich langer Zeit erreichte Endstadium der 
Welt, im letzten wäre umgekehrt das Ausgangsstadium der 
Welt die statische Lösung. 

Analoges gilt für die Maßbestimmung (20) deswegen nicht, 
weil keine statische Lösung mit positiver Massendichte exi- 
stiert, und der Lösung eine negative Dichte zukommt. Diese 
Lösung ist deshalb wohl physikalisch unbrauchbar. Aus dem- 
selben Grunde ist es auch unwahrscheinlich, daß die Lösung 
(19) eine physikalische Bedeutung hat. Was wird aus den 
besprochenen Maßbestimmungen für verschwindende Massen- 
dichte? Aus (17): 


dst=—. sin?) .dı#? 


Massendichte u = 


+ (By + - det 
Aus (20): 
ds?’ = . [dor + Sin?O.d sin?) 


(19) geht in die pseudoeuklidische Maßbestimmung über. 

Die bei der Untersuchung des Falles II gewonnenen 
Lösungen sind deswegen interessant, weil sie wie die Fried- 
mannsche Lösung homogene Dichteverteilungen besitzen. 

Was die Verschiebung der Spektrallinien anbelangt, liegen 
die Verhältnisse, im Gegensatz zu den in $3 besprochenen 
Lösungen, komplizierter. Die Verschiebung wird von der 
Richtung abhängig, aus der das Licht kommt. Vielleicht wird 
sich einmal herausstellen, daß letzteres wirklich der Fall ist. 


II. Teil. 
Spezielle Probleme in nichtstatischen Räumen konstanter Krümmung‘) 


§ 6. Energie und Impulsverhältnisse 
in der räumlich geschlossenen Welt 


Die Erhaltungssätze des Impulses und der Energie sind 
in der allgemeinen Relativitätstheorie in den Gleichungen 


1) Darunter verstehen wir die in $2 gefundenen Lösungen. 
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(22) DS = 0 (n = 1, 2, 3, 4) 


6a™ 
(m) 


enthalten, darin bedeutet 7." die Komponenten des Energie- — 
impulstensors und ¢,™ die des Affintensors des Gravitations- 
falls. Um zur Integralform der Erhaltungssätze zu gelangen, 
integriert man diese Gleichungen über einen geschlossenen 
Raum und erhält, wenn die durch einmalige Integration ent- 
stehenden Doppelintegrale verschwinden: 


(23) zafffV- 


Für n = 1, 2,3 sprechen diese Gleichungen den Impuls- 

satz, für n = 4 den Energiesatz aus. Weiter läßt sich zeigen, 
daß die Energie und Impulsintegrale in weitgehendem Maße 
von der Wahl des räumlichen Koordinatensystems unabhän- 
ig sind. 
m Als abgeschlossenes System haben wir das räumlich ge- 
schlossene Weltall zu betrachten. In der Einsteinschen 
Zylinderwelt verschwinden obige Doppelintegrale. In der nicht- 
statischen sphärischen Welt hingegen wird sich herausstellen, 
daß jene Doppelintegrale nicht sämtlich Null sind. Dies führt 
dazu, daß die Gesamtenergie zeitlich veränderlich ist. Auch 
ist die oben erwähnte Invarianz der Energie gegenüber ge- 
wissen räumlichen Koordinatentransformationen in Frage ge- 
stellt. Der Impulserhaltungssatz ist noch erfüllt. 

Wir wollen dies allgemein für die Maßbestimmung 
d= P?. (dO? + sin?O-(d9° + sin?d - d g)) — wobei 
die Zeitabhängigkeit von P völlig offen gelassen sei, beweisen. 
Die Feldgleichungen ergeben, wenn p den Druck, und o die 
Energiedichte bedeutet: 


2 d’P 


1 fap\ı 1 

3 (dP\?, 3 
(24b) =—At + pr 


Es sind die ¢," zu berechnen. Dies geschieht nach der 
Formel: 


1 i i j 
+2) 
i,j,k,t 
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es folgt: 
1[3 IP\? etg?@ 
1 1/3 ote 
n 0 0 
ap 1 dP etg¢ 


n gibt die Zeilen, m die Kolonnen an. Durch Einsetzen 
der Werte von ¢,” in (22) folgt für n = 4, indem über den 
ganzen Raum integriert wird: 
in? dP 
[ffv- ses +1940 dg =0. 


x 
Bezeichnen wir das dreifache Integral, die negative Gesamt- 
energie mit — E, so ergibt sich: 


(25) E= er - P + const. 


Das Energieintegral werde noch direkt ausgewertet. Für die 
Gravitationsenergie E, folgt 


2n? P® 3 (dP\2 1 

Für die gesamte sönstige Energie 

E„= P® 2+ x +9 

Damit 
(26) E=E,+E,=*2-P. 
In Ubereinstimmung mit (25). Die Gesamtenergie in der nicht- 
statischen räumlich geschlossenen Welt ist somit direkt 


proportional zum Weltradius, und damit notwendig zeitlich 


veränderlich. 
P ist zur Zeit nach den Hubbleschen Messungen mindestens 


von der Größenordnung 102% cm, deshalb = x 105° Erg/sec. 


Betrachten wir die Gl. (22) für n = 1,2,3. Die durch Inte- 
gration über den ganzen Raum entstehenden Doppelintegrale 
verschwinden, und wir werden zu den entsprechenden Gl. (23) 
geführt, die den Impulserhaltungssatz ausdrücken. Die Aus- 
führung der Integrale (n = 1,2,3) 
ergibt, das der Gesamtimpuls des Weltalls gleich Null ist. 
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0 
0 
3 dP\2 etg?@ 
| 


§7. Übertragung der Hertzschen Lösung 
auf nichtstatische Räume konstanter Krümmung 

Wir untersuchen die Wellenbewegung, die von einem 
um den Ursprung eines rechtwinkligen Koordinatensystems in 
Richtung der Achse schwingenden elektrischen Dipol vom 
Moment f(t) ausgestrahlt wird. Die klassische Beschreibung 
liefert, wenn wir uns auf ein Polarkoordinatensystem beziehen, 
dessen Aquatorebene die 2y-Ebene ist, in dem die Maßbe- 
stimmung gilt 


(27) ds? = dr? + + sin? — da®” 

für die Feldkomponenten die Ausdrücke 
M,,= E,= sind - (t ( _ =) 

M,, = E, = M,, = H,, = My = H;, = 0 


Aufgabe dieses Paragraphen soll sein, zu untersuchen, wie 
sich die Beziehungen ändern, wenn wir anstatt (27) die Maß- 
bestimmungen (9a), (10a), (11a) zugrunde legen. Über die Zeit- 
abhängigkeit von P machen wir keine Voraussetzungen. Die 
Rechnungen werden nur für die Maßbestimmung (9a) des 
Kugelraumes durchgeführt. 

Außer dem ungestrichenen Koordinatensystem (9a) ver- 
wenden wir ein gestrichenes, das durch die Transformation 
der Zeitvariablen !) 


1) M.v. Laue, Berl. Ber. 1931, 124. 
Annalen der Physik. 5. Folge. 16. 
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dax* ‚ dax* 
(29) 
aus dem ungestrichenen hervorgeht. In ihm gilt die Maß- 
bestimmung 


(30) ds? = P®-[d@? + sin? O-(d # + sin? g?) — 
Die Maxwellschen Gleichungen für das ladungsfreie 
Vakuum in beliebigen Riemannschen Räumen sind 


ootM=0 und 4ivM=0. 


Setzen wir die Komponentengleichungen im gestrichenen Ko- 
ordinatensystem an, so stellt sich heraus, daß der zeitabhängige 
Faktor P gar nicht explizit auftritt.) 


a) Übertragung der Hertzschen Lösung auf einen Kugelraum 
von konstantem Radius P, 

In infinitesimalen Gebieten kann der Kugelraum als eben 
betrachtet werden. Analytisch drückt sich dies darin aus, 
daß die Maßbestimmung (9a) stetig in (27) übergeht, wenn wir 
uns immer mehr auf die nächste Umgebung des Poles 0 = 0 be- 
schränken. Dabei wird P,-© mit r identisch. Die Maxwell- 
schen Gleichungen im Kugelraum gehen dabei auch immer ge- 
nauer in die des gewöhnlichen Raumes über. Das gleiche gilt für 
entsprechende Lösungen der Maxwellschen Gleichungen. Um 
diese Grenzbeziehung zwischen (9a) und (27) deutlicher zum 
Ausdruck zu bringen, führen wir in (27) anstatt der Variablen 


r und x* neue r’ und ein durch r= >." 
b b 

Aus (27) wird jetzt 
(27a) dd + sin? d-dg?) —da*”). 

Das Koordinatensystem dieser Maßbestimmung entspricht voll- 
kommen dem gestrichenen Koordinatensystem (30). 

Wir übertragen die Hertzsche Lösung (28) auf das ge- 
strichene Koordinatensystem (27a) mit Hilfe der üblichen 
Transformationsbeziehungen. Es genügt die M,,-Komponente 
zu transformieren. Es folgt 

1 P, — 1 — r 
(31) Mi = sin a. 7(P, ==") + 5 f(P, 5 


In diesen Ausdruck muß der entsprechende für den Kugel- 
raum bei Beschränkung auf kleine © übergehen, dabei wird 0 


1) M. v. Laue, Berl. Ber. 1931, 124. 
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mit 7’ identisch. Versuchsweise setzen wir analog zur Lösung 
im ebenen Raum 


(32) M3 = M3; = My = 0. 
Damit liefern die Maxwellschen Gleichungen 
OM, OMy 


= = = 0. 
09 09 
Aus den übrigen Feldgleichungen folgt 
aM, aM; aM, 
(838) + = 93 
ö (sin - M, aM, 
(33b) 35 2) _ sin? @ = 0; 
aM, aM, 
ö (sin O- M, ö (sin. My. 
(33d) sind 5 wy es 2) _ 9, 


Diesen 4 partiellen Differentialgleichungen müssen die Feld- 
komponenten Mj4,M54, Mi, genügen. Durch Eiimination von 
Mi, und My aus (33a), (33b), (33c) ergibt sich für Mj. die 
partielle Differentialgleichung 


1 ö | 1 ö (sin + - Mj.) | aM, 6M, 

Es ist eine Lösung dieser Gleichung anzugeben, die für kleine 

6 mit (31) übereinstimmt. Dies leistet der Ausdruck 


34) sin |. - 9) + — 0))). 

Aus (33b) und (33c) bestimmen sich Mj, und M;, zu 


(34), (35), (36) in Verbindung mit (32) beschreiben den Aus- 
strahlungsvorgang im statischen Kugelraum. 


(36) 


b) Übertragung auf einen Kugelraum von zeitlich 
veränderlichem Radius 
Diese Verallgemeinerung läßt sich in einfacher Weise 
durchführen. Die Einführung des zeitabhängigen Faktors P 
ändert an den Maxwellschen Gleichungen für das gestrichene 
10* 
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System gar nichts. Als Bedingungen haben wir nur noch zu 
erfüllen, daß für P= P, die Lösung im nichtstatischen Raum 
in die des statischen übergeht. Dem wird durch den formal 
den Gl. (34), (35), (36) identischen Ansatz genügt, nur ist z# 
jetzt durch das Integral (29) gegeben. Beim Übergang zum 
ungestrichenen Koordinatensystem ergibt sich 


c 


H,, = sind [#7 - 9) + ))|; 
| = H,, = H,, =0. 
In der Wellenzone wird die Energiedichte der Strahlung 
P, »?-sin? F 

We 
Die Hertzsche Lösung im nichtstatischen euklidischen Raume 
(bezogen auf ds?= P?-[dr?+ r?.(d9?+ sin? dg’) )] — da") 
erhält man aus (37), indem man sin@ = O@=r und cosO=1 
setzt. 

Die Hertzsche Lösung im nichtstatischen hyperpolischen 
Raum mit der Maßbestimmung (11a) erhält man aus (37), 


wenn man die trigonometrischen Funktionen von © durch die 
entsprechenden Hyperbelfunktionen ersetzt. 


§ 8. Zur Lichtausbreitung 


Die Untersuchungen des vorigen Paragraphen gestatten 
uns die Lichtausbreitung genau zu verfolgen. (Intensitätsver- 
hältnisse vgl. § 9). Hier soll speziell die raumzeitliche Aus- 
breitung interessieren. Zunächst einige Definitionen. Unter 
dem Polabstand eines Nebels verstehen wir den Wert seiner 
-Koordinate, und unter der Entfernung in einem Zeitmoment 
das Produkt aus dem Polabstand und dem jeweiligen Momen- 
tanwert der Größe P. 

Beginnt die Ausstrahlung zur Zeit Null, so hat der 
Wellenkopf nach der Zeit t, wie sich aus dem Argument von 

t 


f ergibt, den Polabstand 0 = c f 5 erreicht. Ist P z. B. eine 
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t 
periodische Funktion der Zeit, so geht das Integral f = mit 


0 
wachsender oberer Grenze gegen Unendlich. Das heißt, daß - 
wenn man genügend lange wartet, schließlich auch jeder noch 
so weit entfernte Punkt erreicht wird. Dies ist auch für 
solche Punkte gültig, die so weit entfernt sind, daß die ihnen 
nach (2) zukommende Relativgeschwindigkeit die Lichtge- 
schwindigkeit um beliebiges übersteigt. Dieses Paradoxon er- 
klärt sich dadurch, daß der Lichtstrahl von einem solchen 
Punkte, solange P wächst eine größere Entfernung haben kann, 
als er anfangs hatte. In dem Teil der Periode von P, wo D 
klein, und sogar schließlich negativ wird, wo also auch die 
Relativgeschwindigkeiten ihre Richtung umkehren, hat der 
Lichtstrahl günstige Gelegenheit vorwärts zu kommen. Auch 
in den Fällen (10) und (11) mit 4 = 0, wo P nicht periodisch 
ist, sondern monoton gegen Unendlich geht, erreicht ein Licht- 
strahl schließlich jede Poldistanz. Denn es folgt für (10) 


Beide Ausdrücke werden für unendlich großes P unendlich, 
was zu zeigen war. 

Anders liegen die Verhältnisse in (10) und (11), wenn 
wir das A-Glied positiv voraussetzen. Es gilt die Ungleichung 


t — 
dt 3/1 ä 3 1 
<Vi 
d.h. die Lichterregung bleibt immer innerhalb eines gewissen 
Koordinatenbezirks. Nicht alle Stellen des Weltalls können 
durch elektromagnetische Strahlung miteinander in Wechsel- 
wirkung treten. 

Diese Erscheinung findet ihre Analogie in der De Sitter- 
schen Welt, wo auch das Licht über eine gewisse maximale 
Entfernung nicht hinauskommt. 

Trotzdem kann man am Rande eines solchen Erregungs- 
bezirkes, über den hinaus das Licht nicht kommt, bezüglich 
eines im Koordinatensystem ruhenden Beobachters, nicht von 
einem „Stehenbleiben“ der Lichterregung sprechen, da der 
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Beobachter die Welle immer noch mit der Lichtgeschwindig- 
keit sich fortpflanzen sieht. Daß die Welle im Koordinaten- 
system nicht vorwärtskommt, liegt daran, daß P sehr groß ist. 


§ 9. Rotverschiebung und scheinbare Helligkeit 


Wir vergleichen die scheinbaren Helligkeiten zweier Nebel 
von gleicher absoluter Helligkeit miteinander. @** sei die 


b 
Poldistanz des einen Nebels, = das Verhältnis der Frequen- 


zen des beim im Pole @ = 0 befindlichen Beobachter an- 
kommenden, und von diesem Nebel ausgehenden Lichtes, P* 
und P, seien die Werte des Weltradius bei der Emission und 
der Ankunft des Lichtes. Dem anderen Nebel entsprechen 


die Größen 0%, or P*. Zwecks Beantwortung der Frage 


nach dem scheinbaren Helligkeitsverhältnis, betrachten wir 
von den Nebeln sich ausbreitende Kugelwellen. Den einen 
Nebel machen wir zum Pole eines Koordinatensystems 0, 9,9, 
und untersuchen die Ausstrahlung eines senkrecht zur Ver- 
bindungslinie dieses Nebels mit dem Beobachter schwingenden 
elektrischen Dipols. Die Schwingungsrichtung wählen wir zur 
Richtung += 0 in unserem Koordinatensystem. Die zeitlich 
gemittelte Energiedichte, die ein Maß für die scheinbare Hellig- 


keit ist, wird am Beobachtungsort (9 3) nach (38) 
pr: 
we 


Entsprechend machen wir den anderen Nebel zum Pol eines 
Koordinatensystems ©, 9,9. Für die zeitlich gemittelte Energie- 
dichte am Beobachtungsort folgt ein dem obigen analoger Aus- 
druck, und für das Verhältnis der Helligkeiten ergibt sich 


Oder 
(2) 
H* P, „sin 
H* (P* sin &* 
(=) 
ee 
it den elationen — = 7 un we. olgt 


(39) 
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Im nichtstatischen euklidischen Raum ergibt sich ent- 

sprechend (9 = Ya? + y? + 2?) 
(40) > # 
Den gegenwärtigen Wert von P können wir hier gleich 1 
setzen, dann bedeuten ©* und O** die gegenwärtigen Ent- 
fernungen der Nebel. Man sieht dann, daß das klassische 
quadratische Abnahmegesetz der Helligkeit für kleine Ent- 
b 

fernungen, wo und — nahezu 1 sind, Gültigkeit hat. 

Im Raum ergibt sich 


a” 1 
Bei der Lichtausbreitung von einem Nebel zum Beob- 
achter ändert sich P (wir beschränken uns auf den Fall 4 = 0) 
Der Ausbreitungsvorgang schreitet 


d 
gemäß ( 
nach der Beziehung dO = — = vorwärts (das negative Zeichen 


deshalb, weil das Licht sich, wenn der Beobachter im Pole des 
Koordinatensystems ist, in Richtung abnehmenden © bewegt). 
Die Abhängigkeit des P vom jeweils erreichten © regelt nun 


die Differentialgleichung (35) = P.(P„—- P). Für das Licht 


vom Nebel mit den einmal gesternten Größen folgt durch Inte- 
gration 
b * 2 

Den Nebel mit den zweifach gesternten Größen (den „Hilfs- 
nebel“) nehmen wir so nahe am Beobachter, daß man sin 0* 
durch ©** ersetzen kann. Diese Voraussetzung trifft mit ge- 
nügender Genauigkeit bis zu Entfernungen von 10?* cm zu, 
da nach den Hubbleschen Messungen Entfernungen von der 
Größenordnung 102° cm existieren. Die Entfernung des Hilfs- 
nebels ist während der Zeit, die das Licht braucht, um zu 
uns zu gelangen, nahezu konstant, da »®/»** für Entfernungen 
von 10% cm nur wenig von 1 abweicht. Die Entfernung r** 
des Hilfsnebels ist conn identisch mit P,- ©*, und es folgt 


aus (39) P,- sin O* = 7-2, wo 


4) M. v. Laue, Berl. Ber. 1931. S. 124. 
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Nach der klassischen Theorie wäre z die tatsächliche Ent- 
fernung des gesternten Nebels. Hier bedeutet z lediglich eine 
der Beobachtung zugängliche Größe. Um Klarheit darüber zu 
bekommen, in welchem Intervall sich z bewegt, wenn der ge- 
sternte Nebel seine Lage im Raume ändert, eliminieren wir 
aus (42) und (43) »’/»*, woraus 

P, sin o* 

Der Nenner wird für ein mögliches ©* Null. Dem entspricht 
z=%. Aus seiner Definition folgt, daß z nur unendlich 
werden kann, wenn H* Null ist. Nebel mit solchem ©* sind 
unsichtbar. Ihre Lage im Raum ist derart, daß die Emission 
des Lichtes, das bei uns jetzt ankäme, mit dem Zeitpunkt 
zusammenfällt, wo P den Wert Null hatte. Im allgemeinen 
gehören zu irgendeinem Werte von z zwei Werte von ©*, ein 
größerer und ein kleinerer, als dem unsichtbaren Gebiet ent- 
spricht. Durch Helligkeits- und Rotverschiebungsmessungen 
ist es also nicht möglich, die Lage eines Nebels im Raume 
eindeutig festzulegen. 

Dasselbe läßt sich über z in den nichtstatischen hyper- 
bolischen und euklidischen Räumen aussagen. Im letzteren 
liegt das unsichtbare Gebiet in der Entfernung 


„4.10% em. 


== 


D; 
Elimination von ©* aus (42) und (43) ergibt 
Ay (1+ V1 +2D,-2z) 
v (1+ V1+2D,-z) 2 
4 
Mit K =1. 


Im hyperbolischen Raum führt die analoge Rechnung 
zu (44) mit K =— 1, und im euklidischen Raum zu (44) mit 
K=0. 

Als unbekannter Parameter geht in (44) die Größe P, ein. 
Der Variabilitätsbereich von z ist, wie nun festgestellt (0, + oo). 
Die graphische Darstellung der Beziehungen zwischen Av/r 
und z liefert die Fig. 1. Die Kurven zeigen übereinstimmend 
ein asymptotisches Konvergieren gegen 1, wenn z gegen Un- 
endlich geht. Die Kurven im sphärischen Falle verlaufen 
alle oberhalb, die im hyperbolischen alle unterhalb der durch 
Kenntnis von D, eindeutig festliegenden Kurve im euklidischen 
Falle. Es hat zunächst den Anschein, als ob im hyper- 
bolischen Raum ein gegen Unendlich wachsendes z zu nega- 
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tivem »°/v* führen könnte, was natürlich sinnlos wäre. Dies 
ist aber nicht der Fall, wie man sich durch Entwickeln der 
4. Potenz im Nenner und Beschränkung auf die Glieder 


höchster Ordnung in 2 leicht überzeugt, da nach § 3 P, > ix 


7 
10 "200, (Sohärischer kaum) 
7 
Ye 
/ 
05 Y 
4 
B 
/ = 
| / 4 
L L zZ 


| 
7 2 13 
Fig. 1. Der angeschriebene Abszissenmaßstab 


bezieht sich auf die gestrichelt gezeichneten Kurven. 
Die nicht gestrichelten haben einen 10mal so großen Abszissenmaßstab 


Die oben erwähnte systematische Abweichung der 4»/v- 
Kurven in den verschiedenartigen Räumen gibt prinzipiell 
eine Möglichkeit, die Frage nach dem Vorzeichen der Krüm- 
mung des Raumes zu entscheiden. Allerdings kann dies noch 
wesentlich tieferes Eindringen in den außergalaktischen Raum 
erforderlich machen. 


IIL Teil: Zur Thermodynamik der Strahlung 
in der allgemeinen Relativitätstheorie 
Einleitung 


R. ©. Tolman und P. Ehrenfest!) haben die Hohlraum- 
strahlung in statischen und nichteuklidischen Räumen mit der 
Maßbestimmung 


1) R.C. Tolman u. P. Ehrenfest, Phys. Rev. 36. S. 1791. 
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dst= dada — der = alt 


untersucht und gefunden, daß die Energiedichte infolge der 
Schwere der Energie innerhalb des Hohlraums variiert, und 
zwar wird die Abhängigkeit der natürlich gemessenen Energie- 
dichte u, durch die Beziehung 

(45) v* = const 

geregelt. Abgeleitet ist diese Gleichung mit Hilfe des Energie- 
impulssatzes, unter der Annahme eines isotropen Druckes, der 
ein Drittel der Energiedichte ist. Im folgenden wird eine 
geometrisch-optische Methode angewandt, die gestattet, die 
gesamte Thermodynamik der Strahlung auf beliebige statische 
und solche nichtstatische Maßbestimmungen, deren Strecken- 
räume nur konformen, zeitlichen Änderungen unterworfen sind, 
zu übertragen. 


$ 10. Statische Maßbestimmungen 
a) Ein dem Clausiusschen Theorem entsprechender Satz 


Die Theorie der Strahlenbündel wird beherrscht von dem 
Clausiusschen!) Satz, der besagt: Sind ds, d2, & die Be- 
stimmungsstücke eines unendlich schmalen Strahlenbündels, 
das sich im Vakuum fortpflanzt, so gilt bei der Ausbreitung 


(46) d Q-ds-cos + = const. 


Für den Beweis?) dieses Satzes ist es hinreichend zu ver- 
langen, daß die Lichtstrahlen dem Fermatschen Prinzip ge- 
nügen. Ein solches Prinzip existiert aber auch im statischen 
Gravitationsfeld. Nach Levi-Civita®) fällt dort das vierdimen- 
sionale geodätische Prinzip mit dem Prinzip von der Minimal- 
zeit, bezogen auf die universelle Zeit t, zusammen. Charakte- 
ristisch für Lichtstrahlen ist demnach das Variationsprinzip 


sf =o. 


Wir fassen drei*) Flächen ins Auge, die in dem kleinen 
Bezirk, in dem sie uns interessieren, eben sein sollen, und be- 
zeichnen sie der Reihe nach mit den Indizes 1, 3, 2. Die 
Flächen 1 und 3 seien einander parallel und so nahe bei- 
sammen, daß in dem Raum zwischen ihnen die euklidische 


1) R. Straubel, Phys. Ztschr. 4. S. 114. 

2) R. Clausius, Poggendorf. Ann. 121. S.1. 

3) Levi-Civita, ,,Absoluter Differentialkalkiil 1928, S. 223 ff. 

4) Das folgende schließt sich eng an die Ausführungen von 
Clausius. 
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Geometrie gilt, und die Lichtgeschwindigkeit, auf ¢ um- 
gerechnet, konstant gleich v, ist. In den ebenen Teilen 
unserer drei Flächen wählen wir je ein geradliniges, recht- 
winkliges Koordinatensystem als lokales Bezugssystem aus. - 
Der Ursprung des Systems in 3 sei die senkrechte Pro- 


jektion des Koordinatenursprungs in 1, und entsprechende 


Koordinatenachsen seien parallel. Bezüglich dieser lokalen 
Koordinatensysteme habe ein PunktQ, von 1 die Koordinaten 
2, ¥,- Analog ein solcher Punkt Q, von 3 die Koordinaten 
und ein Punkt Q, in 2 die Koordinaten %,-T,, 
sei die universelle Zeit, die das Licht braucht, um von @, 
nach Q, zu kommen, längs des die Punkte verbindenden 
Strahles. Entsprechende Bedeutung besitzen T,, und T,,. 
Hält man Q, und Q, fest und läßt Q, wandern, so ist der 
Punkt in der Fläche 1, der auf der Verlängerung des Strahles 
Q,, Q, liegt dadurch gekennzeichnet, daß für ihn T,—T,, 
ein Extrem ist. Die 7 sind Funktionen der 6 Koordinaten 
Yır Yo: Zz, Yg- Also folgt als Bedingung für ein 


rem: 
1s = 0, 0. 


(4) ay, 


Zwischen den 6 Koordinaten bestehen diese zwei Gleichungen, 
die, wenn 2 Punkte gegeben sind, den mit ihnen auf einem 
Strahle befindlichen dritten Punkt bestimmen lassen. Jetzt 
werde Q, als fest betrachtet. Die Gleichungen (47) bewerk- 
stelligen dann eine eindeutige Abbildung der Fläche 2 auf 3 
und umgekehrt. Ein Flächenelement ds, in 3 gehe in ds, 
in 2 über. Es gilt 
Ys) wx (2% , 9% _ , 

Die Gleichungen (47) nach x, und y, differenziert, wobei zu 
beachten ist, daß die z,, y,-Funktionen von z,, y, sind, und 
daß T,, von 2, 9, und T,, von 2,, y, unabhängig ist: 

ex, Ox, 02,0%, 0x, 0m, 5 
Drei weitere solche Gleichungen erhält man, wenn man z, 
durch y,, 2, durch y, schrittweise ersetzt. Aus diesen vier 
Gleichungen kann man die oben vorkommende Funktional- 
determinante berechnen. Es folgt dann 
Oy, Ox, G2, O2,0% OY, 


Oy, Or, Ox, OY, Ox, 0x, OY, Oy, 
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In unserem Falle ist T = — (r Abstand zwischen Q, und Q,, 


r= Vo? + (x, — 2” +(y — y,)®, @ der senkrechte Abstand 
von 1 und 3). Die Berechnung des Nenners von (48) ergibt 


Der Zähler sei abgekürzt B,. Weiter ist, wenn #, 


den Winkel zwischen Strahl und Lot auf 1 bedeutet, cos 3, = £, 
ferner ist, wenn dQ, der räumliche Öffnungswinkel des Q, aus- 
gehenden Strahlenbiindels ist d 2, = en. Damit folgt 


dQ, - cos #, = v,*- B,-ds,. In gleicher Weise betrachten wir 
ein von Q, ausgehendes Strahlenbiindel, das den räumlichen 
Offnungswinkel d2, habe, gegen das Lot von 2 den Winkel 9, 
besitzt und von der Fläche 1 ein Flächenelement ds, aus- 
scheidet. Die Fläche 3 befinde sich in derselben Lage gegen- 
über 2 wie oben gegen 1. Durch einfache Vertauschung der 
Indizes 1 und 2 folgt d 2,- cos, = v,?- B,+ds, (v, Licht- 
geschwindigkeit in Q,). Es ist T,, = T,,. Im Ausdruck für B 
sind auch die Differentialoperationen beziiglich der Indizes 1 
und 2 symmetrisch. Damit ist B, = B,, und man findet 
+ cos - ds, 


v;? 
2 
Dies ist die gesuchte Verallgemeinerung des Clausiusschen 
Satzes auf beliebige statische Maßbestimmungen. Sie spielt 
für die Thermodynamik der Strahlung in solchen Räumen die- 
selbe grundlegende Rolle wie jener ursprüngliche Satz unter 
der Herrschaft der pseudoeuklidischen Geometrie. Da die 
Schwingungszahl der universellen Zeit ¢ in statischen Räumen 
bei der Ausbreitung erhalten bleibt, ändert sich die Wellen- 
länge mit dem Orte wie v, und man kann in (49) v durch A 
ersetzen. Dann erkennt man unmittelbar, daß die Anzahl der 
Freiheitsgrade!) eines streng monochromatischen Strahlen- 
bündels sich bei der Ausbreitung nicht ändert. Darin liegt 
eine Bestätigung für (49). 
b) Energie und Entropie des Strahlenbündels 
bei der Ausbreitung. Strahlungsgesetze 
Ein Strahlenbündel mit den Bestimmungsstücken d2, 
ds, #, dem Frequenzintervall v,, »,-+dy,*), umschlieBen wir 


1) Uber den Begriff eines Freiheitsgrades eines Strahlenbiindels 
vgl. M. v. Laue, Ann. d. Phys. 44. S. 1197. 1914. 
2) », die auf t bezogene, v, die auf Ortszeit bezogene Frequenz. 
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durch eine Fläche F und betrachten die Energieverhältnisse 
des Bündels bei der Ausbreitung. Der Energiesatz lautet 


= 0; 


Die Doppelsumme fällt im PR Felde weg. Die Fläche F 
rücken wir so weit vom Flächenelement ds, daß bei dem Aus- 
breitungsvorgang kein Energiestrom durch die Oberfläche von F 
stattfindet. Integration über F liefert 


ga SS SW: - = = . Vg )- 
Da die Energiedichte 
W=-V-94' T*=--—:T#, 
gibt das dreifache Integral die negative REN E,, des Strahlen- 
bündels auf t bezogen. Auf die Ortszeit = 2 -¢ transfor- 


mieren wir die Komponente T,*. Es folgt TY=T,? Die 
örtlich gemessene Energiedichte ist W, =— T*’, somit 


Ein Beobachter, der die Energie des Bündels, am Orte von ds 
mißt, findet für diese den Betrag E,, der zu E, in derselben 
Beziehung steht, wie die Energiedichten W, und W,: 


(50) E, e=E,'v= const. 
Das Bündel werde in der universellen Zeit At emittiert. Der 
Schwingungsvorgang in einem Punkt des Flächenelements kann 
durch eine Fourierreihe 
u= tar, 
A, cos ( 

u=»,. At 
dargestellt werden. Die Frequenz jeder einzelnen Partial- 
schwingung bleibt bei der Ausbreitung erhalten. Das ist nur 
möglich, wenn 
(51) At= const. 


Ein bei ds befindlicher Beobachter wird eine auf seine Ortszeit r 
bezogene spezifische K, so definieren, daß 


E,=K,'d2.cos &-ds.—- At. 
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Entsprechend kann man eine auf ¢ bezogene Strahlungs- 
intensität K, definieren, daß 

E,=K,-d2-ds-cos#- At. 
Nach (49), (50), (51) wird: 


(52a) K, - v* = const 
und 
(52b) K+ v? = const 


K,, und &,, seien die spektralen Strahlungsintensitäten. Da 
„.c=9 


ist auch 
(53) er = const und dv, = const. 
0 
Aus (52a) und (52b) folgt 
(54a) Boo = const 
und 
(54b) RK. -v? = const. 


Die Gl. (49) lehrt, daß die Anzahl der Freiheitsgrade eines 
streng monochromatischen Strahlenbiindels erhalten bleibt. 
Die Freiheitsgrade eines Biindels von der spektralen Breite d», 
erhält man durch Multiplikation der Freiheitsgrade des streng 
monochromatischen mit 4t-dy,. Nach (51) und (53) bleiben 
also auch die Freiheitsgrade eines beliebigen Bündels erhalten. 


Die Gl. (53) in der Form geschrieben i; = const driickt 
aus, daB die Anzahl der Lichtquanten konstant bleibt. Die 


Entropie ist nun nach Planck ebenfalls eine Invariante bei 
der Ausbreitung, und die Lichtausbreitung im statischen 
Schwerefeld ist ein reversibler Vorgang. 

S, und S, seien die auf Orts- bzw. universelle Zeit be- 
zogenen Entropien des Strahlenbündels, Z, und L, die ent- 
sprechenden spezifischen Strahlungsintensitäten. Es ist dann 


=L,-dQ-cos-ds-. At 
c 


und 
S,=L,-d2.cos# ds. At. 


Da S, = S8,, folgt L,-v = L,-¢ und daraus L,, = L,,. Weiter 
nach dem Ulausiusschen Satz wegen S = const, und (51) 


(55a) L,,-v? = const oder = const, 
(55b) L,, = const. 
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Daher nach (54) 


u 


wie 
67a) =» se), also = const, 


(7b) 


also T,= const. 


Umkehrung 8, und & = 25 (2). 


ce) Anwendung auf Hohlraumstrahlung 


Soll die Strahlung im Gleichgewicht sein, so muß zu jedem 
Strahlenbündel im Hohlraum ein reziprokes existieren, das das 
erstere gerade kompensiert. Dies ist der Fall, wenn in jedem 
Punkt die Gl. (52) erfüllt sind. Die Strahlung ist damit isotrop. 
Da auf Ortszeit bezogen immer c, und auf t bezogen v die 
Lichtgeschwindigkeit ist, so ergibt sich 


(58a) u, + v* = const 
und 
(58b) u, = const. 


Dies ist das Ergebnis von Tolman und Ehrenfest. 

Uber die Temperaturverteilung im Hohlraum geben die 
Gl. (57) AufschluB.') Die natürlich gemessene Temperatur T, 
ändert sich umgekehrt proportional zu v, während die uni- 
verselle Temperatur 7’, im ganzen System konstant ist. Das- 
selbe Gesetz hat R. (. Tolman?) auch für andere thermo- 
dynamische Systeme abgeleitet. Die Gültigkeit läßt sich auch 
im Falle eines Gases im Schwerefeld nachweisen. Die uni- 
verselle Temperatur 7, haben wir mit der klassischen Tem- 
peratur zu identifizieren. Diese ist proportional zur mittleren 
kinetischen Energie der Moleküle, und transformiert sich beim 

bergang von universellen Einheiten zu Ortseinheiten wie die 

Energie, die sich nach (50) umrechnet. 

In üblicher Weise kann das Stefan-Boltzmannsche 
Gesetz abgeleitet werden: 


ud 


1) Vgl. A. Einstein, Ann. d. Phys. 38. S. 443, 1912. 
2) R. C. Tolman, Phys. Rev. 15 S. 904, 1930. 
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Ebenso das Wiensche Verschiebungsgesetz des Strahlungs- 
maximums: 


4,°T,=b und 7, = +b. 


§ 11. Nichtstatische Maßbestimmungen 


Es gelte die Maßbestimmung (1) (P sei jedoch nur von 
der Zeit abhängig) 

ds? = P?. Sy, ‚da* dat —da® = P2.do? — de”. 
Führen wir hier nach (29) die gestrichene Zeitvariable ein, und 
stellen wir im gestrichenen Koordinatensystem die Maxwell 
schen Gleichungen auf, so geht in diese, wie wir wissen, deı 
zeitlich veränderliche Faktor P nicht ein. Die Lichtstrahlen 
sind also dieselben, wie die der Maßbestimmung 


ds? = do? —dzx*’’. 


In dieser statischen Maßbestimmung gelten die Überlegungen 
des vorigen Paragraphen. Den Clausiusschen Satz erhält 
man aus (49), wenn man v durch die auf die gestrichene Zeit 
umgerechnete Lichtgeschwindigkeit ersetzt. Die gestrichene 
Zeiteinheit ist P-mal so groß, als die ungestrichene. Damit 
ist die Lichtgeschwindigkeit auf die gestrichene Zeit P-c. 
Damit folgt als Clausiusscher Satz 


Ein Bündel ändert seine Energie gemäß ') 
(60) P.E= const. 


Die Schwingungszahl bezüglich der gestrichenen Zeit bleibt 
erhalten. Nach ähnlichen Überlegungen wie im vorigen 
Paragraphen folgt daraus, daß die auf die gestrichene Zeit 
bezogene Emissionszeit des Bündels At’ konstant ist, also 
(61) At’ = const. 

Die jetzt folgenden Betrachtungen sind genau analog denen 
des vorigen Paragraphen. Der universellen und Ortszeit dort 
entspricht hier die gestrichene und ungestrichene Zeit. An die 
Stelle von v dort tritt hier P. So folgt analog zu (54) 


(62) Fr = const 


(8, die auf die ungestrichene Zeit t bezogene Strahlungs- 
intensitat). Für die Frequenzänderung gilt das Gesetz 


P-v = const. 


1) M. v. Laue, Berl. Ber. 1931. S. 128. 
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Die Wellenlänge ändert sich demnach proportional zu P, und 
man erkennt aus (59) unmittelbar, daß die Anzahl der Freiheits- 
grade des streng monochromatischen Strahlenbündels sich nicht 
ändert. Die Freiheitsgrade des Bündels von der spektralen. 
Breite dv’ (auf die gestrichene Zeit bezogen) erhält man durch 
Multiplikation der Anzahl des streng monochromatischen mit 
dv - At’. Dieser Faktor bleibt konstant. Demnach ist die 
Anzahl der Freiheitsgrade eines beliebigen Bündels eine In- 
variante bei der Ausbreitung. Die Beziehung (60) kann auch 


= const, und sagt aus, daß die 


so geschrieben werden: 


- Anzahl der Lichtquanten erhalten bleibt. Damit bleibt auch 
die Entropie des Bündels unverändert, und die freie Licht- 


ausbreitung ist ein reversibler Vorgang. Damit folgt + = const. 
Mit (62) schließen wir auf 
L, 1 
Das Wiensche Verschiebungsgesetz gilt in unveränderter Form 
auch in nichtstatischen Räumen. 
Die Temperatur ergibt sich wie bei Planck zu 
7-1. also T.P= const. 
Die so definierte Temperatur des Bündels ändert sich um- 
gekehrt proportional zu P. 


Es ist mir ein Bedürfnis, auch an dieser Stelle meinem 
hochverehrten Lehrer, Hrn. Prof. v. Laue für die Anregung 
und stetige Förderung der vorliegenden Arbeit ergebenst zu 
danken. 


Berlin, Institut für theoretische Physik. 


(Eingegangen 1. November 1932) 
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Magnetostriktion ferromagnetischer Kristalle 
Von Richard Gans und Jürgen von Harlem’) 
(Mit 8 Figuren) 


Die Längenänderung von Eisen- und Nickel-Einkristallen 
im Magnetfelde kann nach denselben Prinzipien behandelt 
werden, nach denen wir die Widerstandsänderungen berechnet 
haben.?) Dasselbe Problem ist für Sättigung bereits von 
Akulov?) und für ungesättigte Kristalle von Heisenberg‘) 
behandelt worden, jedoch ist ihr aus dem Symmetrieprinzip 
folgender Ansatz nur eine erste Näherung, worauf wir bereits 
in unserer vorigen Mitteilung (a. a. O.) hinwiesen. Diese Tat- 
sache und der Wunsch, das Beobachtungsmaterial möglichst voll- 
ständig zu diskutieren, veranlaßt uns, hier auf die Frage noch 
einmal zurückzukommen. 

Es handelt sich also — ebenso wie bei der Darstellung 
der Widerstandsänderungen — zunächst um die Berechnung 
gewisser linearer Kombinationen von Materialkonstanten h aus 
Messungen an magnetisch gesättigten Kristallen. Dann hat 
man mittels dieser unter Hinzuziehung der Heisenbergschen 
Hypothese die Erscheinungen für ungesättigte Ferromagnetika 
als Funktion der Magnetisierung vorauszuberechnen und mit 
der Erfahrung zu vergleichen. 

Trotzdem nun der allgemeine Ansatz für die Berechnung 
der Längenänderungen unter dem Einfluß der Sättigungs- 
magnetisierung ganz analog dem für die Widerstandsänderung 
ist, sind trotzdem die Formeln für die Transversaleffekte nicht 
einfach aus der früheren Veröffentlichung zu übernehmen, weil 
die Versuche in verschiedener Weise angeordnet waren. 

Während nämlich bei den Widerstandsmessungen der 
Transversaleffekt derart beobachtet worden ist, daß die Richtung, 
in der die Widerstandsänderung beobachtet wurde, fest blieb, 
die Magnetisierung dagegen in der dazu senkrechten Ebene 


1) Für die numerischen Rechnungen trägt v. Harlem die Ver- 
antwortung. 
2) R. Pikens u. J. v. Harlem, Ann. d. Phys. [5] 15. S. 516. 1932. 
3) N.S. Akulov, Ztschr. f. Phys. 52. S. 389. 1929 
4) W. Heisenberg, Ztschr. f. Phys. 69. S. 287. 1931. 


Gan: 

| alle 
Lan 
selb 
dreh 
| rech 
Mag 
so V 
aber 
Län; 
eine 

(1) 

die 
For: 
men 
drüc 
Mag 
| gem 
| sind 
geni 
rech 
| sind 
sich 
stril 
es ¢ 

stan 

f jede 
| gebi 
ders 
! Da 

zieh 
(N) 
[vgl. 
im | 


ctor 


- ı ~— 


Gans u. v. Harlem. Magnetostriktion ferromagnet. Kristalle 163 


alle méglichen Azimute annahm, wurde bei den transversalen 
Längenänderungen so verfahren, daß die Verlängerung in der- 
selben Ebene gemessen wurde, in der die Magnetisierung ge- 
dreht wurde, natürlich so, daß die Verlängerung immer senk- 
recht zur Magnetisierung beobachtet wurde. Nahm also die 
Magnetisierung J alle möglichen Azimute in dieser Ebene an, 


so wurde also die Verlängerung A= = stets senkrecht zu J, 


aber auch in dieser Ebene gemessen, d.h. die transversale 

Längenänderung wurde nicht, wie die Widerstandsänderung, in 

einer festen, konstant bleibenden Richtung beobachtet. 
Natürlich läßt sich Öl/l aus der Formel 


A= = hy + hy By? + +++) + hy By Bg + 
(1) + hs (@,* + + hy By? + ose) 


die aus dem Symmetrieprinzip folgt und genau analog der 
Formel fiir die Widerstandsänderung ist, für jeden experi- 
mentellen Spezialfall durch die Materialkonstanten h aus- 
drücken, da ja stets die Richtungskosinus @,, @,, «, der 
Magnetisierung sowie die Richtungskosinus ß,, ß,, 6, der 
gemessenen Längenänderung gegen die Hauptachsen angebbar 
sind. Die hierbei auftretenden linearen Kombinationen der h 
genügen, wie bei den Widerstandsänderungen, nicht zur Be- 
rechnung der einzelnen h-Werte, weil zur Messung stets be- 
sonders ausgezeichnete Richtungen und Ebenen benutzt worden 
sind. Es zeigt sich aber, daß dieselben h-Kombinationen, die 
sich aus den Sättigungsmessungen ergeben, für die Magneto- 
striktion im nicht gesättigten Zustande maßgebend sind, so daß 
es doch möglich ist, die Erscheinungen im ungesättigten Zu- 
stande vorherzusagen und mit der Erfahrung zu vergleichen. 

Im makroskopisch unmagnetischen Zustande ist bei Nickel 
jede räumliche Diagonale für die Magnetisierung der Elementar- 
gebiete gleich wahrscheinlich, d. h. die « in (1) haben alle mit 


derselben Wahrscheinlichkeit den Wert + —— oder — Mr a 


v3 
Da dann aber definitionsmäßig A= 0 ist, hd folgt die Be- 
ziehung 
(1) h, + Ath no 


[vgl. die analoge Gl. (10”) unserer früheren Mitteilung], die wir 
im folgenden häufig gebrauchen werden. 


it” 


| 


164  Annalen der Physik. 5. Folge. Band 16. 1933 


Dagegen gilt für Eisen, bei dem die Hauptachsen die 
energetisch bevorzugten Richtungen sind 


(1”) hy +h, + 
I. Nickel 


1. Longitudinaleffekte bei Sättigung 
a) Liegt die Längenänderung in der (001)-Ebene, so folgt 
aus (1) 
(2) A= a, + -sin? 29 
mit 
| a, =h, +h, +h, +h, = — 50,0-10-*, 
a, =— 2h, +h, — 2h, — 3h, = 68,7 - 10. 

Die angegebenen Zahlenwerte sind dabei so gewählt, daß 
die Ergebnisse der Beobachtungen möglichst gut dargestellt 
werden. 

b) Liegt die Längenänderung in der (110)-Ebene, so ist 
(3) A=a, +a, (sin? } — 3 sint F) + Za, (sin? $ — sin® 9) 
mit 
63) a,=h,+h +h,. 

Setzt man in dieser Formel sin? + = 2, d.h. 9 = 50°44, 
das ist die [111]-Richtung, so ergibt sich 

Aus den Messungen in der (110)-Ebene fiir diese Richtung 
folgt aber A = — 27,1-10-* Somit folgt unter Beachtung 
der Zahlenwerte für a, und a, [Formel (2’)] 

(3”) a, =h,+h, +h, =0, 

eine Beziehung, die wir im folgenden benutzen werden, und 
die auch in ganz analoger Weise bei den Widerstands- 
änderungen gilt. 

c) Liegt die Längenänderung in der (111)-Ebene, so gilt 
mit Berücksichtigung von (3”) 

(4) A= a + 4 
Fig. 1 gibt die mittels der beiden Konstanten [vgl. (2)] 


berechneten theoretischen Kurven sowie die von Mashiyama’) 
beobachteten Punkte. 


1) Y.Mashiyama, Sc. Rep. Töhoku 17. S. 955. 1928. 
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Fig. 1. Nickel, Longitudinaleffekte bei Sättigung 


2. Transversaleffekte bei Sättigung 


a) Längenänderung und Magnetisierung liegen in der 


(100)-Ebene. Dann ist 
- 2 
(5) A =a, + — sin 29. 
y6 
8 
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Fig. 2. Nickel, Transversaleffekte bei Sättigung 
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Mit 
5) a, =h, +h, = 24 10%, 
é | a, = 2h, —h, — 2h, +h, =— 39,9-10-*. 
b) Liegen sie in der (110)-Ebene, so wird 
(6) A=a,—4(3a, + a, + 6a, — 2a,)sin? — ja, sin* +. 

c) Liegen sie schlieBlich in der (111)-Ebene, so wird 

a, a a 

(7) A=—($ +444): 

In Fig. 2 sind die Beobachtungen dargestellt. Sie stimmen 
in der (100)-Ebene recht gut mit der Theorie überein, dagegen 
zeigen sich in den anderen beiden Ebenen doch merkliche 


Abweichungen, ein Beweis für die Unvollkommenheit der be- 
nutzten Kristalle. 


3. Longitudinaleffekte als Funktion der Magnetisierung 
Die Formeln für die Longitudinaleffekte in ihrer Ab- 
hängigkeit von der relativen Magnetisierung 
J 
8) n= 


lassen sich wegen der Analogie zwischen Magnetostriktion und 
Widerstandsänderung direkt aus unserer früheren Arbeit über- 
nehmen. 
Für die[001]-Richtung gilt demnach nach a.a. O. Formel (9), 
(10) und (10) 
1 
fü 0 
rV<n< V3 
(9) A=0, 
1 
fir —— 2: 
r V3 <1< 
| A =— +a, + 2a, + + ($a, +a, + 2a, + 4,) 9° 
—i@+4, + 2a, + a,)n*. 
Für die [110]-Richtung nach a. a. O. Formel (7) 
fir 0<n< v3: 


3 


für Vz 


(a, + 1”, 


(13) 


(13) 


[110 
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Schließlich ist nach a. a. O. Formel (6) und (6’) für die 
[111]-Richtung 


a, a 
(11) A= + 1-1). 
06 
Mr 
Dr = 
Hr a 
L100] £110] L111) 

Wr = = 

-I. a 
aH 


Fig. 3. Nickel, Longitudinaleffekte als Funktion der Magnetisierung 


Fig. 3 zeigt die gute Übereinstimmung zwischen Theorie 
und Erfahrung. Es ist allerdings zu bemerken, daß diese 
Übereinstimmung bei einer anderen Probe schlechter ist. 


4. Transversaleffekte als Funktion der Magnetisierung 


a) J sei parallel der [111]-Richtung, A parallel [112]. 
Dann gilt 


(12) A =— (fa, + + 1 — 1). 
b) Ist J parallel der [001]-Richtung, A parallel [010], so wird 


(13) A=0, 
f 
ür y3 <n<l 
A = 1 (4a, + a, +4a,+a,) — 3(4a,+a,+5a,+a,)n? 


as) | 
+ (4a, + a, + 8a, + a,)7*. 


Die gleiche Formel gilt, wenn J parallel [001], A parallel 
[110] ist. 
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Fig. 4 gibt den Vergleich der Theorie mit der Erfahrung. 
Man sieht, daß die Übereinstimmung für die Magnetisierung 
in Richtung der leichtesten Magnetisierbarkeit am schlech- 
testen ist. 
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Fig. 4. Nickel, Transversaleffekte als Funktion der Magnetisierung 
c) Ist J parallel [110], A parallel [001], so wird 
fir O0<n< 2: 
(14) A=0, 
fiir 
| A=a,+4(— 2a, + a, + 2a, + a,) 7? 
— 3(a, + a, + 2a, + a,) 7°. 
d) Ist J parallel [110], A parallel [110], so gilt 
fir 0O<n< 
(15) A =— 1(3a, + 
fiir <n<l: 
| A =a, — 3(3a,+a, +6a,+3a,)n? 
+3(a,+a, + 6a, + a,)n*. 


Fig. 5 zeigt, daß die Theorie den erfahrungsmäßig fest- 
gestellten Verlauf einigermaßen gut wiedergibt. 
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Fig. 5. Nickel, Transversaleffekte als Funktion der Magnetisierung 


II. Eisen 
1. Longitudinaleffekte bei Sättigung 


Hier gelten die gleichen theoretischen Formeln wie für 
Nickel, trotzdem die Richtungen leichtester Magnetisierbarkeit 
bei Fe — im Gegensatz zu Ni — die Hauptachsen sind. Die 
in Frage kommenden Zahlenwerte sind 


a,=16-10-*; a, =—87-10-*; a, =0. 


Es ist bemerkenswert, daB man sich dem gesamten Beob- 
achtungsmaterial am besten anpaßt, wenn man — ebenso wie 
bei Nickel — auch beim Eisen 


a,=h+h,+h,=0 


setzt. Diese Beziehung, welche innerhalb unserer phänomeno- 
logischen Behandlungsweise durchaus nicht begründet werden 
kann, verdient vielleicht besondere Beachtung, wenn man das 
Problem einmal als gittertheoretisches, magnetoelastisches 
angreift. 

Fig. 6 zeigt, daß man mit obigen Zahlenwerten guten 
— an die Beobachtungen von Honda und Mashiyama}?) 
erhält. 


1) K.Hondau. Y. Mashiyama, Sc. Rep. Tohoku (1) 15. S.757. 1926. 
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° 
(100) (10) (711) 


Fig. 6. Eisen, Longitudinaleffekte bei Sattigung 


2. Transversaleffekte bei Sättigung 
Liegen Längenänderung und Magnetisierung in der (100) 
Ebene, so gilt 


(16) A=— + sin? 2¢ 
mit 
a, = 88 . 10-#; 
liegen sie in der (110)-Ebene, so wird 
(17) A=- 5 + a, sin? — sin* 
mit 


a,’ = 81-10-*. 
Liegen sie schlieBlich in der (111)-Ebene, so ergibt sich 
(18) A =— ,\, (6a, — 8a,’ + 5a,). 
Hier gelten fiir a,, a,, a,, a, dieselben Ausdrücke in 
den h, wie bei Nickel [vgl. (2'), (3'), (5)). Dagegen ist 
(19) a, = $h, 
Macht man von der Beziehung (1”) Gebrauch, so wird 
a =F 
a, =— 2h, +h, — 2h, — 3h,, 
a, th, th, 2h, + 3h,+ ıh,, 
a, = 2h, —h, — 2h, +h,. 


(20) 
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Diese 4 Gleichungen genügen natürlich nicht, um die h 
zu bestimmen, es kommen aber in den Formeln für die 
Magnetostriktion ungesättigter Kristalle als Funktion der 


7” 


° 


° 
| am 


Fig. 7. Eisen, Transversaleffekte bei Sättigung 


(100) 


Magnetisierung wiederum nur die in (20) auftretenden linearen 
Kombinationen der h vor, die sich aus den Sättigungsmessungen 
ergeben. 

Fig. 7 zeigt recht große Abweichungen zwischen den Beob- 
achtungen von Honda und Mashiyama und der Theorie, 
ein Beweis für die Unvollkommenheit der benutzten Kristalle 
oder der Messungen. 


3. Longitudinaleffekte als Funktion der Magnetisierung 


Unter Heranziehung der von Heisenberg’) ermittelten 
Verteilungsfunktionen berechnet sich für die Längenänderung 
in der [100]-Richtung 
(21) A =a, (V3? + 1 —1); 
in der [110]-Richtung 

1 
fir 0 ——: 
<n< 
(22) 


1) W. Heisenberg, Ztschr. f. Phys. 69. S. 287. 1931. 


0 oo ke | 
)- 
= 


172 Annalen der Physik. 5. Folge. Band 16. 1933 


für ye 


(22) A= >. + —a,)« VI-e:, 
wo « eine Abkürzung für 
— 
v2 
bedeutet. 
In der [111]-Richtung ist 
für 0 
rV<n< y3 
(23) A=0, 


A= (6a, + 2a, + 4a,’—3a,)+ (—8a, +a, +2a,)ß? 
(23))  — jy (4a, +40, +a, (124, + 8a, — 8a, +50) 


1-6? 1, 1 — 


wo # eine Abkürzung ist für 
+ 


© Honda und Mashiyame 
x Webster 


Fig. 8. Eisen, Longitudinaleffekt als Funktion der Magnetisierung 
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In Fig. 8 bedeuten die Kreuze Beobachtungen von 
Webster’), die Kreise Messungen von Honda und Ma- 
shiyama.?) Weitere Beobachtungen der japanischen Forscher 
an anders geschnittenen Proben sind fortgelassen worden; diese - 
schließen sich den Beobachtungen weniger gut an. Es sei aber 
besonders betont, daß bei der Berechnung der in Fig. 8 dar- 
gestellten theoretischen Kurven die Konstanten a nicht etwa 
so gewählt worden sind, daß die Beobachtungen möglichst gut 
wiedergegeben werden, wie Heisenberg (a. a. O.) das bei Aus- 
wertung seiner Näherungsformel getan hat, sondern die a sind 
aus den Messungen bei Sättigung ermittelt, es stand also keine 
Konstante frei zur Verfügung. 


Zusammenfassend kann man sagen, daß die Heisenberg- 
sche Hypothese sich auch bei der Magnetostriktion so gut be- 
währt, wie man erwarten kann. 


1) W. Webster, Proc. Roy. Soc. (A) 109. S. 570. 1925. 

2) K. Honda u. Y. Mashiyama, a. a. O. 

Königsbergi. Pr., II. Physikalisches Institut, 24. No- 
vember 1932. 


(Eingegangen 27. November 1932) 
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Anderung der Gitterkonstante 
und der Leitfähigkeit des Palladiums 
durch elektrolytische Wasserstoffbeladung 


Von F. Krüger und G. Gehm 
(Mit 5 Figuren) 
Einleitung 


Der Zustand des im Palladium gelösten Wasserstoffs ist 
trotz zahlreicher Untersuchungen!) keineswegs völlig geklärt, 
Zu den älteren und neueren Bestimmungen des Absorptions- 
gleichgewichts des Wasserstoffs im Palladium als Funktion 
der Temperatur oder des Druckes?) sind in den letzten Jahren 
auch Untersuchungen der Leitfähigkeit?) des Palladiums in 
Abhängigkeit von der Wasserstoffbeladung, dann vor allem 
aber auch röntgenographische Untersuchungen‘) der Beeinflussung 
der Gitterstruktur bzw. der Gitterkonstanten des Palladiums 
durch die Wasserstoffbeladung hinzugekommen. 

Die Beladung des Palladiums mit Wasserstoff kann ent- 
weder rein thermisch in einer Wasserstoffatmosphire erfolgen 
oder elektrolytisch. Die thermische Beladung erfordert, wenn 
sie bei den betreffenden Drucken jeweils bis zur Sättigung, 
also bis zum Gleichgewicht getrieben werden soll, bei tieferen 
Temperaturen sehr lange Zeiten, und auch dann ist ein Kri- 


terium, ob auch wirklich das Gleichgewicht eingetreten ist, 


1) 1834—1900, E. Bose, Ztschr. f. phys. Chem. 34. S. 701. 1900; 
1900—1921, L. W. Me. Keehan, Phys. Rev. 21. S. 334. 1923. 

2) L. Troost u. D. Hautefeuille, Compt. rend. 78. S. 668. 1874; 
Chem. phys. 2. S. 273. 1874; C. Hoitsema, Ztschr. f. phys. Chem. 
17. S.1. 1895; A. Holt, E.C. Edgar u. I. B. Firth, Ztschr. f. phys. 
Chem. 82. 8.513. 1913; A.Sieverts, Ztschr. f. phys. Chem. 88. 
S.103 u. 451. 1914; B. Lambert u. S. F. Gates, Proc. Roy. Soe. A. 
108. 8.456. 1925; L.J. Gillespie u. F. P. Hall, Journ. Amer. Chem. 
Soc. 48. 8.1207. 1926. 

3) F. Fischer, Ann. d. Phys. 20. 8.503. 1906; G. Wolf, Zitschr. 
f. phys. Chem. 87. S.575. 1914; B. Beckmann, Ann. d. Phys. 46. 
8.481. 1915; A. Puod2iukynas, Ztschr. f. Phys. 46. S.253. 1927; 
A.Coehn u. H. Jürgens, Ztschr. f. Phys. 71. 179. 1931. 

4) Me. Keehan, Phys. Rev. 21. 8.334. 1923; M. Yamada, 
Phil. Mag. 45. S.241. 1923; F. Krüger u. A. Sacklowsky, Ann. d. 
Phys. 78. 8.72. 1925; I. D. Hanawalt, Phys. Rev. 33. S. 444. 1929; 
E. Nähring, Diss, Greifswald 1928. 
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schwer zu finden, denn die Grenze, bei der keine Absorption 
mehr stattfindet, ist wenig scharf ausgeprägt. Deshalb sind 
die meisten Absorptionsgleichgewichte bei höheren Temperaturen 
in der Nähe von 100°C und darüber angestellt; so hat 
A.Sieverts!) in seinen eingehenden Untersuchungen einer- 
seits die Absorption des Wasserstofis von Atmosphärendruck 
in Abhängigkeit von der Temperatur und zwar sowohl für 
Palladiumdraht (auch für flüssiges Palladium bei 1600°) und 
für Palladiumschwamm und -Mohr bis zu etwa 900° gemessen, 
ferner die Wasserstoffabsorption durch Palladiumdraht wie 
auch durch Palladiumschwamm und -Mohr in Abhängigkeit 
vom Gasdruck (bis etwa 1 Atmosphäre) bei verschiedenen kon- 
stanten Temperaturen zwischen 183—821° untersucht. Genaue 
Messungen der Druckkonzentrationsisothermen für 75, 103 und 
120° haben B. Lambert und p 
8. F. Gates?) bis zu Drucken von 
etwa 2 Atmosphiiren vorgenommen e 
und gezeigt, wie schwer es ist, 
richtige Gleichgewichte zu erhalten. b 
Isothermen fiir ein weiteres Tem- 
peraturgebiet von 0—200 bzw. 180° a 
hatten schon früher C. Hoitsema’) 
(und Roozeboom) und ferner Fig.1. Druck-Konzentrations- 
später L. I. Gillespie und F. P. Diagramm des Wasserstoffs im 
Hall!) aufgenommen. Auf die ver- Palladium 
schiedenen Schlüsse, die diese Autoren aus ihren Kurven ziehen, 
soll hier nicht eingegangen werden. Hier sei nur bemerkt, daß 
die Druckkonzentrationsisothermen die typische Gestalt besitzen, 
die die nebenstehende Fig. 1 zeigt. Bei geringer Konzentration 
(Kurventeil a) gilt die Formel c = k Vp (c = Wasserstoffkonzen- 
tration, k = Konstante, p = Druck des Wasserstoffgases'. Dann 
folgt in dem c-p-Diagramm meist ein fast horizontaler Teil 
(Kurventeil b), der von einer bestimmten Konzentration (Kurven- 
teil c) rasch wieder ansteigt. Aus der Verschiebung des End- 
punktes der horizontalen Kurve mit der Temperatur ergibt 
sich mit Sicherheit der Schluß, den auch Hoitsema mit aller 
Klarheit gezogen hat, daß die in dem dem horizontalen Kurven- 
teil entsprechenden Gebiet vorhandenen zwei Phasen feste 
Lösungen sein müssen. 

Für die Beladung auf elektrolytischem Wege sind vor 
allem die maximalen absorbierten Mengen bestimmt, die von 
Palladium nach Überladung bei Zimmertemperatur freiwillig 


1) a.a. O. auf $. 174. 
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Tabelle 1 

Autoren | H/Pa 
F. Krüger und A. Saklowsky*) . . . . 954 
A. Coehn und H. Jürgens) . . . . . 1030 


behalten werden. Die folgende Tabelle 1 gibt einen Überblick 
dieser von verschiedenen Autoren für Palladium in Form von 
Blechen, Drähten oder Schwamm erhaltenen Maximalbeladungen. 

Röntgenographische Bestimmungen der Gitterstruktur von 
wasserstoffbeladenem Palladium sind zuerst für die Maximal- 
beladung vorgenommen worden. Das wesentlichste Resultat 
dieser Messungen war, daß durch die Wasserstoffbeladung 
keine Änderung des Gittertypus eintrat, sondern nur eine 
Gitterdehnung; sowohl das unbeladene wie das wasserstofi- 
beladene Palladium besitzt ein flächenzentriertes, kubisches 
Gitter. Auch die Stärke der Gitterdehnung bei der Maximal- 
beladung wurde von verschiedenen Autoren fast gleich gefunden, 
wie die folgende Tabelle 2 zeigt, in der die erste Spalte die 
Gitterkonstante des unbeladenen Palladiums, die zweite die 
bei der maximalen Beladung anzeigt. 


Tabelle 2 
| Gitterkonst. Gitterkonst. 
Autoren N d. unbel. d. bel. 
|  Palladiums in AE 
F. Krüger und A. Saklowsky®) . . 3,89 4,00 


Bei Linde und Borelius’?) schwankt die Gitterkonstante 
des beladenen Gitters von 4,032—4,065 ÄE. 


1) H. Lisenko, Ber. d. dtsch. chem. Ges. 5. S. 29. 1872. 

2) M. Thoma, Ztschr. f. phys. Chem. 3. S. 69. 1889. 

3) L. Nowack, Ztschr. f. anorg. u. allg. Chem. 113. $.1. 1920. 

4) a.a. O. auf S. 174. 

5) H.O. v. Samson-Himmelstjerna, Ztschr. f. anorg. u. allg. 
Chem. 186. S. 337. 1930. 

6) a.2.0. auf S. 174. 
7) J. O. Linde u. G. Borelius, Ann. d. Phys. 84. S. 747. 1927. 
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Die Abweichungen der Werte voneinander erklären sich 
teils daraus, daß ja, wie die vorhergehende Tabelle 1 zeigt, 
die Maximalbeladungen nicht immer genau gleich ausfallen, 
teils aus der Ungenauigkeit der älteren röntgenographischen - 
Messungen. 

Bei verschiedenen Wasserstoffdrucken und damit im Gleich- 
gewicht stehenden Wasserstoffkonzentrationen im Palladium 
haben zuerst Linde und Borelius!) röntgenographische Be- 
stimmungen der Gitterkonstante vorgenommen. Die Beladung 
fand thermisch in einer Wasserstoffatmosphäre bei Drucken 
zwischen 0,17 und 4 Atmosphären und den Temperaturen 100, 
150 und 200° © statt. Ihre Messungen bei 100° gibt die 
nebenstehende Fig. 2 wieder, 
in der (umgekehrt wie in ihrer 4 
Arbeit) als Abszisse die Wurzel %%@r ° 
aus den Drucken und als Or- 
dinate die Gitterkonstanten 
aufgetragen sind. 397 ° 

Aus der Fig. 2 ersieht 
man, daß bei mittleren Drucken 
zwei Phasen existieren mit je 391 2 
einer charakteristischen Gitter- sa 
konstanten, eine bei geringen 
Drucken mit einem kleineren pn 
Gehalt an Wasserstoff undeine * V7 YF 
zweite bei höheren Drucken Fig. 2. Gitterkonstanten als Funk- 
miteinemgrößeren Wasserstoff- tion der Wurzel aus dem Druck 
gehalt; erst bei einem Gasdruck 
von etwa 2 Atm. tritt ein Anwachsen der zweiten Gitterkon- 
stante ein. Welche Wasserstoffkonzentration diesen Drucken 
entspricht, haben Linde und Borelius nicht angegeben. 

Diese Autoren haben ferner auch die Gitterkonstanten 
des beladenen Palladiums bei elektrolytischer Wasserstofi- 
beladung bestimmt und zwar bei 18 und bei 90° C. Sie haben 
jedoch nur die Gitterkonstante jeweils nach Wasserstoff beladung 
einer gemessenen Zahl von Stunden und gegebener Stromstärke 
gemessen, nicht aber die dann vorhandene Wasserstoffkonzen- 
trationen bestimmt. So finden sie zwar auch hier, daß die 
höheren Wasserstoff konzentrationen sich über einen zweiphasigen 
Zwischenzustand ausbilden, sie können aber die zu den Mittel- 
werten der Gitterkonstanten der beiden Phasen zugehörigen 
Konzentrationen nur aus einer immerhin gewagten Extrapolation 


1) a.a. 0. auf S. 176. 
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und dem bei hoher Wasserstoffkonzentration gemessenen Zu- 
sammenhang zwischen Konzentration und Gitterkonstante er- 
schließen: danach würde bei 18° der wasserstoffärmeren Phase 
die Konzentration c, = 0,03 Wasserstoffatome auf ein Atom 
Palladium mit der Gitterkonstante 3,894 AE und der zweiten 
Phase die Konzentration c, = 0,70 mit der Gitterkonstante 
4,034 ÄE entsprechen. Bei 90° war die Gitterkonstante der 
wasserstoffärmeren Phase von der des reinen Palladiums nicht 
zu unterscheiden; die Konzentration der wasserstoffreicheren 
Phase entsprach danach c, = 0,65. 

Die vorliegenden Messungen bezweckten einmal die zu 
direkt gemessenen Wasserstoffkonzentrationen gehörigen Gitter- 
konstanten zu bestimmen und ferner einen strengen Beweis 
dafür zu erbringen, daß durch die Wasserstoffbeladung nur 
die Gitterkonstante vergrößert, also das Gitter gedehnt wird, 
daß aber der Gittertypus unverändert bleibt. Davon war dann 
eine klarere Deutung der komplizierten Verhältnisse des Systems 
Palladium—Wasserstoff zu erhoffen. 


1, Bestimmung der Gitterstruktur und der Gitterkonstanten des 
Palladiums in Abhängigkeit von der elektrolytischen 
Wasserstoff beladung 


a) Versuchsanordnung und Meßverfahren 


In der vorliegenden Arbeit wurde diese röntgenographische 
Untersuchung an Palladium’) vorgenommen, das bei Zimmer- 
temperatur (20°) elektrolytisch mit Wasserstoff beladen wurde. 
Dabei wurde die Wasserstoffbeladung allmählich gesteigert und 
bei je einer definierten Beladung jedesmal eine Röntgenauf- 
nahme gemacht. 

Das verwandte Palladium hatte die Form eines Drahtes 
mit quadratischem Querschnitt bei einer Kantenlänge dieses 
Querschnittes von 1 mm. Die Drahtlänge betrug 5—7 cm. 
Die Einzelbeladungen dauerten im Anfangsstadium etwa !/, Std. 
und stiegen mit zunehmender Beladung bis zu etwa 10 bis 
12 Std. Bei der hohen Diffusionsgeschwindigkeit des Wasser- 
stoffs im Palladium ist anzunehmen, daß sich eine gleichmäßige 
Beladung des Palladiums über den ganzen Querschnitt herstellt. 

Die Stäbchen wurden vor der Beladung zwecks Homo- 
genisierung und Entfernung von Gasen etwa 6 Std. bei einer 
Temperatur von 750° im Vakuum geglüht; hierauf wurde Luft 


1) Für die freundliche Überlassung des Palladiums sind wir der 
Firma G. Siebert G. m. b. H. in Hanau zu Dank verpflichtet. 
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in das Quarzrohr eingelassen und dann dieses zur schnellen 
Abkühlung in Wasser getaucht. 

Die Bestimmung der Menge des absorbierten Wasserstofts 
fand wie bei Berry’), in folgender Weise statt: Die zu unter- 
suchenden Stäbchen wurden elektrolytisch bei Zimmertemperatur 
beladen unter Verwendung von 2 Voltametern, die hinter- 
einander geschaltet waren. Der an den Kathoden sich ent- 
wickelnde Wasserstoff wurde in 2 Büretten aufgefangen, die 
eine Ablesegenauigkeit von 0,1 cm? besaßen. Die eine Kathode 
bestand aus dem Palladiumstab, dessen Wasserstoffabsorption 
gemessen werden sollte. Die andere Kathode und die > 
bestanden aus Platin. Der Elektrolyt war 5°/,ige Schwefel- 
säure, und die Stromstärke betrug 0,006 Amp./cm?. Die Difte- 
renz der entwickelten Wasserstofimenge in den beiden Büretten 
ergab dann die absorbierte Menge Wasserstoff, die noch auf 
0° und 760 mm Druck umgerechnet wurde. 

Um den Wasserstoff wieder aus dem Metall zu entfernen, 
standen zwei Möglichkeiten zur Verfügung. Einmal konnte 
das Palladium in dem Voltameter als Anode geschaltet werden, 
so daß sich der sich entwickelnde Sauerstoff gleich mit dem 
Wasserstoff zu Wasser verband. Zweitens konnte der Wasser- 
stoff durch Glühen im Hochvakuum wieder aus der Legierung 
herausgetrieben werden. Bei der ersten Methode ist man nicht 
ganz sicher, ob nicht geringe Spuren Wasserstoff zurückbehalten 
werden, die der Messung entgehen. Darum wurde die zweite 
Methode bevorzugt. 

Als Röntgenröhre wurde ein Siemens-Reiniger-Veifa- 
Strahlenschutzrohr mit Kupferantikathode benutzt. Die Röhre 
wurde bei 21 kV eff. mit 20 mA belastet. Die Aufnahmen 
wurden in zwei verschiedenen Kameras gemacht. Mit der 
einen, die nach Angabe von E. Nähring in der Werkstatt 
des hiesigen Physikalischen Institutes angefertigt war, konnten 
sowohl sämtliche Glanzwinkel zwischen 0 und 90° auf einem 
zylindrischen Film als auch nur die Glanzwinkel nahe bei 90 
oder bei 0° auf zwei Platten aufgenommen werden. Da die 
Platten einen größeren Abstand vom Präparat als der zylin- 
drische Filmhalter hatten, konnten mit den Platten genauere 
Aufnahmen erzielt werden. Die Belichtungszeit betrug bei dem 
zylindrischen Film 3 Std. und bei den Platten 4 Std. 

Die andere Kamera war ebenfalls in der Institutswerkstatt 
nach einer Beschreibung von G. Sachs und J. Weerts?) an- 


1) A. J. Berry, Journ. Chem. Soc. London. 99. S. 463. 1911. 
2) G. Sachs u. J. Weerts, Ztschr. f. Phys. 60. S. 481. 1930. 
12* 
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gefertigt. Nur wurde sie, um die Belichtungszeit herabzusetzen, 
nach Angaben von E. Nähring insoweit abgeändert, daß statt 
der Bleiblende ein Spalt, bestehend aus zwei mit Platin be- 
‚stäubten Glaswänden, benutzt wurde. Die Breite dieses Spalts 
betrug 1 mm und die Höhe 5 mm. Die Länge der bestäubten 
Glaswände war 12cm und reichte sehr dicht an das Linde 
mannfenster heran. Durch Totalreflektion an den Innenwänden 
dieser Spiegelblende gelang es, eine größere Intensität an das 
Präparat zu bringen, und so die Belichtungsdauer, die sonst 
wegen der großen Entfernung Antikathode—Priparat 9—12 Std. 
fordert hätte, auf 4 Std. herabzudrücken. Sämtliche Präparate 
mußten bei der Aufnahme gedreht werden (eine Umdrehung 
pro Minute), da die Kriställchen sehr groß waren. 

Die Auswertung geschah derart, daß auf den Glanzwinkel 
90° extrapoliert wurde, wie es in den Arbeiten von G. Kett- 
mann!) und G. Sachs und J. Weerts?) angegeben wird. Mit 
dieser Kamera war es leicht möglich, eine Genauigkeit von 
0,001 AE und noch mehr zu erreichen; für die vorliegende Unter- 
suchung genügte aber die Genauigkeit von 0,001 ÄE. Die 
benutzten Wellenlängen für die Auswertung wurden mit 
Ka, = 1,5373 und Ke, = 1,5430 AE angesetzt. 


b) Meßresultate 


Mit dieser Anordnung wurden nun die Gitterkonstanten 
des Palladiums gemessen, nachdem es in Intervallen elektro- 
lytisch mit einer Wasserstoffmenge beladen war, die sich, wie 
oben auseinandergesetzt, aus der Differenz der beiden Volta- 
meterablesungen dabei ergab. Das Resultat gibt die unten- 
stehende Tabelle 3 wieder. In der ersten Kolumne stehen die 
Wasserstoffkonzentrationen umgerechnet in H/Pd, in der zweiten 
die aus den Röntgenogrammen berechneten Gitterkonstanten. 


Tabelle 3 
Absorb. H, Gitterk. Absorb. H, | Gitterk. 
in H/Pd in AE in H/Pd in AE 
0,0 | 3,884 0,56 3,886 
0,11 3,886 4,019 
023 3.886 0,66 | 4.020 
| 4.020 0,79 | 4.024 
0,36 | 3,887 0,85 | 4,037 
4,020 0, | 4.045 


'1) G. Kettmann, Ztschr. f. Phys. 53. S. 198. 1929. 
2) a.a.0. auf S. 179. 
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In Fig. 3 sind diese Ergebnisse graphisch dargestellt, auf 
der Abszisse sind die Wasserstoffkonzentrationen, auf der 
Ordinate die Gitterkonstanten aufgetragen; in der Mitte sind 
in der Kurve C noch 


ie weiter unten be- A 
405 


sprochenen Ergeb- 

nisse der Leitfähig- B DA w 
keitsmessungen darge- 
stellt. Die Figur zeigt, 400+ 17,80 
daß in einem großen 160 
Konzentrationsgebiet 
zwei Phasen mit je 391 
einer charakteristi- 41420 
schen Gitterkonstante 1400 
vorhanden sind, die 3o0L 

in den Kurven A und rc 

B aufgetragen sind. 335 05 7 


Sie entsprechen den 
beiden von Linde Fig.3. Gitterkonstanten und Widerstands- 
und Borelius ge- verhältnis als Funktion der Konzentration 


fundenen Phasen. 

Die Röntgenogramme des wasserstoffbeladenen Palladiums 
zeigen bei beiden Phasen dieselben Linien wie die des reinen 
Palladiums, nur sind sie mit zunehmender Wasserstoff beladung 
weiter auseinander gerückt, d.h. das Gitter wird durch die 
Beladung nur gedehnt, aber in seiner Struktur nicht geändert. 
Die der Gitterdehnung entsprechende Linienverschiebung zeigt 
deutlich Fig. 4, in der das Photogramm a dem unbeladenen 
Zustand des Palladiums, Photogramm b dem des maximal be- 
ladenen entspricht. Man ersieht die Linienverschiebung vor 
allem bei den großen Glanzwinkeln am Rande (bei Cu, 420). 
Der Linienreichtum erklärt sich daraus, daß die Röntgenröhre 
neben der Kupferstrahlung auch Eisenstrahlung lieferte. 

Um zu sehen, ob nicht etwa eine Änderung der Inten- 
sitätsverhältnisse der Linien oder gar ein Ausfall bestimmter 
Linien doch auf eine Gitteränderung hindeutet (Übergang des 
flächenzentrierten kubischen Gitters in ein raumzentriertes 
Gitter) wurden die Röntgenfilme mit dem Zeissschen re- 
gistrierenden Photometer ausgemessen. Um eine Ungleich- 
heit der Intensität infolge ungleicher Absorption wegen 
des quadratischen Querschnitts, der zur Aufnahme be- 
nutzten Stäbchen auszuschließen, wurden die Stäbe während 
= Aufnahme gedreht. Die erhaltenen Kurven zeigt Fig. 5a 
und b. 
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Da nach der Definition der Schwärzung S = lg- + ist, wo 


J, die einfallende und J, die durchgelassene Intensität be- 
deutet, wurde der Logarithmus des Verhältnisses des Elektro- 
meterausschlages des kontinuierlichen Untergrundes durch den 
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Fig. 4a. Pd unbeladen 
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Fig. 4b. Pd beladen 


Elektrometerausschlag des Maximums der Schwärzung gebildet 
und die Schwärzung der intensivsten Linie gleich 10 gesetzt 
und alle anderen darauf bezogen. Da es sich um absolute 
Schwärzungen unter 0,6 handelt, kann man die so gemessene 
Schwärzung proportional der Intensität setzen, und man er- 
hält in dieser Weise die Intensitätsverhältnisse der einzelnen 
Linien. Die so gewonnenen Zahlen gibt Tab. 4 wieder. In 
Kolumne I stehen die relativen Intensitäten des unbeladenen 
Palladiums, in Kolumne II die des beladenen. 
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Der Vergleich der beiden Spalten liefert die Verhält- 
nisse der Intensitäten der Linien, die im ganzen nicht wesent- 
lich verändert sind. Die Abweichungen dürften hauptsächlich 
auf Ungenauigkeiten der Ausmessung beruhen. Jedenfalls - 
dürften auch diese Messungen den Schluß bestätigen, daß durch 


Fig. 5a. Pd unbeladen 


| | 


Fig.5b. Pd beladen 


Tabelle 4 
Relative Intensitäten der Linien 
I | II 
Spektrall 
| (anbel. Pd) | (bel. PA) 
Cu 111 5,8 5,3 
Cu 200 6,8 | 6,0 
Cu 20 . 10,0 10,0 
Cu 311 9,3 8,5 
Cu 331 | 8,2 7, 


Cu 40 7,2 6,7 
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die Wasserstoffbeladung keine Änderung des Gittertypus, son- 
dern nur eine Dehnung des Gitters stattfindet. 

Die Phase mit der geringeren Wasserstoffbeladung 
(Kurve A) zeigt nur eine sehr geringe Gitterdehnung, nämlich 
um 0,002 AE; bei der angegebenen Genauigkeit von 0,001 AE 
ist die Dehnung zwar nur eben noch, aber doch mit Sicher- 
heit, nachweisbar. Linde und Borelius haben für diese 
Phase eine Dehnung um 0,004 ÄE gefunden, was bei der 
Grenze ihrer Meßgenauigkeit um + 0,002 ÄE als hinreichende 
Übereinstimmung anzusehen sein dürfte. Bei höheren Tempe- 
raturen scheint die dieser Phase entsprechende Dehnung zu- 
zunehmen, denn für bei 100°C thermisch beladenes Palla- 
dium geben Linde und Borelius eine Gitterdehnung um 
0,008 AE an. 

Die Gitterdehnung bei der Phase mit großem Wasser- 
stoffgehalt beträgt 0,136 AE, Linde und Borelius finden für 
diese Phase elektrolytisch beladenen Palladiums bei 18° C 
0,144 ÄE. 

Die Intensität der Linien der -wasserstoffarmen Phase 
nimmt mit zunehmender Beladung immer mehr ab, während 
gleichzeitig die der wasserstoffreichen Phase zunimmt; d.h. 
also, daß sich die zweite Phase auf Kosten der ersten bildet. 
An der in Fig.3 mit D bezeichneten Stelle tritt dann ein 
Knick in der Kurve für die wasserstoffreichere Phase auf, die 
Gitterkonstante steigt von da ab nahezu linear an, d.h. von 
diesem Punkt an erfährt das Gitter eine der Beladung pro- 
portionale Dehnung. Dieser Knick liegt bei den vorliegenden 
Messungen bei einer Beladung von c = 0,77 H/Pd. 

Der Existenzbereich der beiden Phasen mit je einer be- 
stimmten Gitterkonstante entspricht dem horizontalen Teil der 
Druckkonzentrationsisotherme in Fig. 1. Der erste Teil der 
allmählichen Gitterdehnung von dem Wert für reines Palla- 
dium bis zu dem der wasserstoffarmen Phase (Kurventeil a in 
Fig. 1) kann röntgenographisch nicht verfolgt werden. Die zu 
dem Beginn des horizontalen Teils der Isotherme gehörige 
Konzentration gibt natürlich die Zusammensetzung der wasser- 
stoffarmen Phase an. Aus den von verschiedenen Autoren 
gemessenen Isothermen liegt dieser Beginn bei Zimmertempe- 
ratur etwa zwischen 0,03—0,1 H/Pd. Von dem Punkt an, in 
dem die Isotherme horizontal wird, müssen natürlich beide 
Phasen vorhanden sein; es wurde versucht, diese Grenze, ge- 
kennzeichnet durch das Auftreten der zweiten wasserstof- 
reicheren Phase (da dann ja die erste Phase vollständig aus- 
gebildet sein muß) durch röntgenographische Messungen zu 
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bestimmen. Damit der Wasserstoff möglichst gleichmäßig im 
Palladium verteilt war, wurde dazu eine sehr dünne Palladium- 
folie von etwa 0,01 mm Dicke benutzt. Sie wurde bei einer 
sehr kleinen Stromdichte (0,0001 Amp./cm?) sehr langsam be- 
laden. Bei einer Konzentration von 0,03 H/Pd wurden noch 
keine Linien der zweiten Phase gefunden; dies war dagegen 
der Fall bei einer Konzentration von c= 0,11 H/Pd. Das 
ist natürlich nur eine obere Grenze, da unterhalb dieser Kon- 
zentration die Linien mangels ausreichender Stärke sich dem 
Nachweis entziehen konnten. Daß nicht etwa eine ungleich- 
mäßige Verteilung des Wasserstoffs im Palladium (etwa vor- 
wiegend an der Oberfläche) das Resultat beeinflußt hat, geht 
daraus hervor, daß ein Palladiumstab von 1 mm Dicke den- 
selben Grenzwert ergab. Eine Konzentration von 0,03 H/Pd 
würde also ein Verhältnis von etwa 30 Pd-Atomen zu einem 
H-Atom, eine solche von 0,11 H/Pd etwa 9 Pd-Atome zu einem 
H-Atom entsprechen. Zwischen diesen Werten, wahrscheinlich 
näher dem ersten, dürfte die Zusammensetzung der wasser- 
stoffärmeren Phase liegen. 

Die Zusammensetzung der wasserstoffreicheren Phase ist 
gegeben durch den Knickpunkt, bei dem eine Zunahme des 
Wertes der Gitterkonstanten einsetzt. Dieser Knick liegt bei 
den vorliegenden Versuchen bei der Konzentration c = 0,77 H/Pd, 
bei Linde und Borelius ebenfalls bei Zimmertemperatur bei 
c=0,70. Einem Verhältnis von c= 0,5 H/Pd würde einer 
Zusammensetzung der zweiten Phase im Verhältnis 2Pd:H 
entsprechen. Da bei einer Temperatur um 100° C herum sich 
dies Verhältnis für den genannten Knickpunkt findet, hat man 
daher öfters eine Verbindung Pd,H annehmen zu müssen ge- 
glaubt. Bei tieferen Temperaturen aber ist das Verhältnis, 
wie oben erwähnt, sicher nicht vorhanden, nach Hoitsema’) 
auch bei höheren Temperaturen nicht, so daß die Identifizie- 
rung der wasserstoffreicheren Phase mit der Verbindung Pd,H 
abzulehnen ist. (Näheres vgl. in einer späteren Arbeit von 
F. Krüger.) 

Linde und Borelius haben für den letzten ansteigenden 
Ast der Kurve die Abhängigkeit der Gitterkonstante von der 
Wasserstoffkonzentration bei elektrolytischer Beladung be- 
stimmt. Mittels der durch vier so bestimmte Punkte (der 
fünfte bei hoher Temperatur erhaltene Punkt ist nach den 
Autoren selbst nicht unbedingt einwandfrei) und den Koordi- 
natenanfangspunkt gelegten Kurve berechnen die genannten 


1) Vgl. a. a. O. auf S. 174. 
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Forscher dann die Konzentrationen für die erste und zweite 
Phase und finden bei Zimmertemperatur für jene c= 0,03 
für diese c= 0,70 H/Pd. Obwohl die Extrapolation aus der 
für den ansteigenden Ast abgeleiteten Beziehung zwischen 
Gitterkonstante und Konzentration recht hypothetisch ist, da 
keineswegs feststeht, ja vielleicht nicht einmal wahrscheinlich 
ist, daß die Art der Löslichkeit des Wasserstoffs im Palladium 
in dem ansteigenden Ast und in den beiden Phasen dieselbe 
ist, stimmen doch die für die Konzentration in den beiden 
Phasen gefundenen Werte für die Zusammensetzung der beiden 
Phasen mit den oben angegebenen Werten recht gut überein. 


2. Leitfähigkeit des wasserstoffbeladenen Palladiums 
in Abhängigkeit von der Wasserstoffbeladung 


Die Abhängigkeit der Leitfähigkeit des wasserstoffbeladenen 
Palladiums als Funktion des Wasserstofigehalts ist schon 
wiederholt exakt gemessen. Wenn diese Messung hier für die 
elektrolytische Beladung des Palladiums mit Wasserstoff noch- 
mals angestellt wurde, so geschah das nicht in der Absicht, 
die Genauigkeit zu steigern, sondern um die Leitfähigkeit an 
demselben Material und bei derselben Konzentration zu be- 
stimmen, für die auch die Gitterkonstanten bestimmt wurden. 


a) Versuchsanordung 


Die zu untersuchenden Palladiumstiibchen hatten einen 
quadratischen Querschnitt von etwa 1 mm Kantenliinge und 
eine Länge von 6 cm. Die Stäbchen wurden an ihren Enden 
in zwei Klemmschrauben eingespannt und zur Potential- 
abnahme waren Stücke von Rasierklingen so montiert, daß die 
Schneiden sich in einem Abstand von 41,5 mm fest an das 
Stäbchen anpreßten. Um eine konstante Temperatur zu er- 
halten, wurde das zu messende Stäbchen in ein Gefäß mit 
Petroleum getaucht, das von außen mit Eis gekühlt wurde; 
die Messungen sind also bei 0°C angestellt. Ein Rührer 
sorgte für gute Zirkulation des Petroleums. Auf diese Weise 
war es möglich, die Temperatur bis zu 4 Std. auf 1/,,° C kon- 
stant zu halten. 

Mit Hilfe eines thermokraftfreien Diesselhorstschen 
Kompensationsapparates!) wurden dann die Palladiumstäbe 
durchgemessen. Obgleich die Widerstandsmessungen am Kom- 


1) Für die freundliche Überlassung dieses Apparates sprechen wir 
Hrn. Dr. H. J. Seemann unsern verbindlichsten Dank aus. 
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pensationsapparat sich mit einer Genauigkeit von 0,1°/, leicht 
durchführen ließen, ist die Genauigkeit der weiter unten mit- 
geteilten spezifischen Widerstände wegen des unsicheren 
Querschnitts weit geringer. 


b) Meßresultate 
Die Meßergebnisse sind in Tab. 5 wiedergegeben: 


Tabelle 5 
Absorb. H, | Spez. Leitfähigk. Spez. Widerst. | ww 
in H/Ppd | -10~* 
0,0 8,61 11,6 1,00 
0,11 7,90 12,7 1,10 
0.23 7,00 14,3 1,23 
0,36 6,55 | 15,3 1,32 
0,56 5,72 17,5 1,51 
0,66 5,45 18,4 1,58 
0,79 5,08 | 19,6 1,69 
0,85 5,80 | 19,6 1,69 
0,91 5,09 | 19,5 1,68 


Die letzte Kolumne, die das Verhältnis des bei 0° C ge- 
messenen Widerstands des wasserstoffbeladenen Palladiums zu 
dem des wasserstofffreien wiedergibt, ist in Fig. 3 in Kurve C 
wiedergegeben. 

Der gefundene Kurvenverlauf entspricht genau dem vor 
allem von F. Fischer?) und von B. Beckmann’) gefundenen, 
d.h. es gibt zwei Unstetigkeiten, die erste bei etwa 30 Volum- 
teilen Wasserstoff, die zweite bei etwa 925 Volumteilen Wasser- 
stoff auf 1 Volumen Palladium, d. h. umgerechnet bei etwa 
0,03 H/Pd und 0,77 H/Pd. Zwischen diesen beiden Punkten 
ist der Kurvenverlauf linear und proportional der aufge- 
nommenen Menge Wasserstoff. Ein solcher Verlauf ist zu er- 
warten, wenn man, wie oben auseinandergesetzt, zwei Phasen, 
bestehend aus zwei gesättigten Lösungen von Wasserstoff in 
Palladium, annimmt; die erste Phase ist gebildet bei 0,03 H/Pd 
(vorher löst sich der Wasserstoff nach dem Quadratwurzel- 
gesetz bis zur Bildung dieser gesättigten Lösung), die zweite 
gesättiste Lösung hat sich ganz ausgebildet bei 0,77 H,/ Pd. 
Dazwischen ist ein Gemenge von Kristallen beider gesättigten 


1) Vgl. a. a. O. S. 174. 
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Lösungen vorhanden, das erfahrungsgemäß einen linearen 
Aufstieg der Leitfähigkeit als Funktion der Zusammensetzung 
aus den beiden Komponenten liefert, hier also die Proportio- 
nalität mit der Wasserstoffbeladung ergibt. 

H. Wolf!) glaubt bei einem Volumgehalt von 600 Volumen 
Wasserstoff eine Diskontinuität in der Leitfähigkeitskurve ge- 
funden zu haben, die einer Verbindung Pd,H entsprechen 
würde. Wie aber B. Beckmann!) gezeigt hat, ergibt eine 
genaue Prüfung, daß eine solche Diskontinuität nicht existiert. 
Wenn man nach Beckmann die Wolfschen Messungen der 
Leitfähigkeit als Funktion der Atomkonzentration aufträgt, so 
wird das besonders deutlich; der mittlere Kurvenverlauf ist 
dann linear, die Abweichungen von der Linearität betragen 
höchstens 2,5°/,, eine Diskontinuität bei 600 Volumteilen tritt 
nicht auf, vielmehr sagt Beckmann: „Der Verlauf der Kurve 
ist derselbe wie bei den Metallegierungen, die aus Konglo- 
meraten der beiden Komponenten bestehen (z. B. Cadmium— 
Bleilegierungen).“ Der Anstieg jenseits des zweiten Knick- 
punktes entspricht dem Anstieg der Gitterkonstante propor- 
tional der Wasserstoffbeladung. Der Anstieg des Widerstandes 
in diesem Gebiet hält jedoch nicht an, sondern geht dann, 
wie A. Coehn und H..Jürgens') zeigten, in eine Wider- 
standserniedrigung über, wie es etwa bei der Bildung einer 
festen Lösung mit zunehmender Konzentration zu erwarten ist. 

In der vorliegenden Arbeit ist der Widerstand stets an 
den gleichen wasserstoffbeladenen Stäbchen gemessen, an denen 
auch die Gitterkonstanten gemessen wurden. Die Schlüsse, die 
aus beiden Meßreihen sich ergeben, stimmen überein: es existieren 
in dem mittleren Gebiet der Beladung des Palladiums mit Wasser- 
stoff zwei Phasen, gebildet aus zwei gesättigten Lösungen. 


Zusammenfassung 


Es wurden gleichzeitig die Gitterkonstante und der 
Widerstand von auf elektrolytischem Wege wasserstoffbela- 
denem Palladium als Funktion der Konzentration des Wasser- 
stoffs im Palladium gemessen. Die röntgenographischen Auf- 
nahmen und ihre photometrische Auswertung zeigten in dem 
mittleren Konzentrationsintervall von etwa c, = 0,03 H/Pd an 
bis zu etwa c, = 0,77 H/Pd die Existenz von zwei Phasen an, 
die beide ein gedehntes Palladiumgitter, aber keinen neuen 
Gittertypus zeigen; denselben Schluß auf die Existenz zweier 


1) Vgl. a. a. O. 8. 174. 


F. | 
solc 
stal 
J.( 
beic 
suc 
stot 
bes 
Dr. 
22. 


F. Krüger u. G. Gehm. Änderung der Gitterkonstanten usw. 189 


solchen Phasen in dem mittleren Intervall ergaben die Wider- 
standsmessungen. Diese Ergebnisse stimmen mit denen von 
J.0. Linde und E. Borelius in vieler Hinsicht überein. Die 
beiden Phasen entsprechen jedoch nach unseren Unter- 
suchungen zwei verschiedenen gesättigten Lösungen des Wasser- 
stoffs im Palladium, nicht aber chemischen Verbindungen in 
bestimmten Verhältnissen. 


Für ständige Beratung und Unterstützung sind wir Hrn, 
Dr. E. Nähring sehr zu Dank verpflichtet. 


Greifswald, Physikalisches Institut der Universität, 
22. August 1932. 


(Eingegangen 26. Oktober 1932) 
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Gitterkonstanten und 
elektrische Leitfähigkeit elektrolytisch beladener 
Palladium-Silberlegierungen in Abhängigkeit 
von der Wasserstoffbeladung 


Von F, Krüger und G. Gehm 
(Mit 7 Figuren) 


Die Löslichkeit des Wasserstoffes in den Palladium—Silber- 
legierungen ist bei geringer Silberkonzentration sehr groß, z.T. 
sogar größer als im reinen Palladium, um dann bei einer Silber- 
konzentration von, etwa 70°/, auf Null herabzusinken. Die 
Untersuchung der Änderung der Gitterkonstanten der Palladium— 
Silberlegierungen durch Wasserstoffaufnahme, verbunden mit 
gleichzeitigen Widerstandsmessungen, verspricht daher weitere 
Aufklärung über den Zustand des gelösten Wasserstofies in 
diesen Legierungen. 


1. Die Gitterkonstanten der unbeladenen 

Palladium-Silberlegierungen und das Vegardsche Gesetz 

An der Mischkristallreihe Palladium-Silber hatten F.Krüger 
und A. Sacklowski!) 1925 das Vegardsche Gesetz bestätigt. 
Da damals aber die Meßgenauigkeit noch gering war, schien 
eine Nachprüfung dieses Befundes erwünscht. Die Versuchs- 
anordnung der Gitterkonstanten war dieselbe wie die in der vor- 
stehenden Arbeit für das wasserstoff beladene Palladium benutzte, 
deren Genauigkeit 0,001 ÄE betrug. Vor der Messung wurden 
die Legierungen?) zwecks Homogenisierung und Entfernung von 
Gasen im Hochvakuum mehrere Stunden zwischen 750 und 
800° getempert. Das Tempern mußte sehr vorsichtig geschehen, 
da bei einer zu hohen Temperatur wegen des relativ hohen 
Dampfdruckes des Silbers sich leicht die Zusammensetzung 
der Legierungen ändert. Nach dem Glühen wurden die 
Legierungen in Wasser von Zimmertemperatur abgeschreckt 
und gemessen. 


1) F.Krüger u. A. Sacklowski, Ann. d. Phys. 78. S. 72. 1925. 

2) Für die liebenswürdige Überlassung der Palladium-Silberlegie- 
rung sprechen wir der Firma G. Siebert G. m. b. H. in Hanau unseren 
verbindlichsten Dank aus. 
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Die gefundenen Werte gibt die folgende Tab. 1 an, die in der 
I. Kolumne die Gewichtsprozente Silber, in der II. Kolumne die 
Atomprozente Silber, in der III. Kolumne die gefundenen Gitter- 
konstanten, in der IV. Kolumne die nach Vegard berechneten 
Gitterkonstanten, in der V. Kolumne die Differenz zwischen 

er beiden enthält. 


t Tabelle 1 
I Iv | Vv 
Gewichts-°/, Atom-°/, Gemess. Berech. _ Differenz 
Silber Silber Giitterkonst. Gitterkonst. der Gitterkonst. 
er- 0 tae 00 3,884 3,884 | 
T 10 990. . 3899 3,9084 | 0,004 - 
15 14,86 3908 | 3,181 | 0,005 
oa 20 19,83 3,916 3,9228 0,007 
die 30 29,77 3,933 3,9422 0,009 
40 39,73 3,950 3,9616 0,012 
nit 50 49,73 | 3970 | 3,9810 | 0,011 
55 64,73 3,980 3,9907 0,011 
Te 60 59,74 3,90 40004 0,010 
in 70 69,78 4,011 4,0198 0,009 
80 79,84 4,030 4,0392 0,009 
90 89,92 4,052 40586 0,007 
100 100,0 4078 | 20780 | a 
Die graphische Darstellung zeigt Fig. 1; sie zeigt ebenso 
er wie die Tabelle, daß zwar kleine aber deutliche Abweichungen 
st. vom Vegardschen Gesetz vorliegen. Diese Abweichungen haben 
en bei Mischkristallen mit _ 
- 50 zu 50°/, und in deren 4 
- Nihe die gréBten Werte 405+ 
8, mit 0,01 AE. 
Während die Gitter- 
konstanten aller Legie- yool 
id rungen eine Abweichung 
a, vom Vegardschen Ge- 
setz zeigten, hatten die 
g Gitterkonstanten der Le- "| 
1e gierungen mit 20 und 
kt 60 °/, Silber anfangs 
ziemlich die Größe, wie 33 
sie das Vegardsche Ge- 338 l > 
5. setz verlangt, ergeben, 
é- 
2 obw ohl die Kurve der Fig. 1. Gitterkonstanten der 
Leitfähigkeit (vgl. unten, Pd-Ag-Legierungen 
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S. 199) an diesen Stellen keine Unstetigkeit zeig. Diese Un-! 
regelmäßigkeiten hätten vielleicht von einer fälschen Angabe 
der Zusammensetzung der Legierungen herrühren können. Eine 
von der Firma G. Siebert G. m. b. H. in Hanau dürchgeführte 
Analyse ergab jedoch eine auf wenige Hundertstel Prozent 
genaue Übereinstimmung der Zusammensetzung mit dem an- 
gegebenen Werte. Die Unregelmäßigkeit mußte also eine andere 
Ursache haben. Nun zeigte sich, daß nur eine dünne Schicht 
an der Oberfläche eine größere Gitterkonstante hatte; entfernte 
man die Oberflächenschicht durch Schmirgeln, so ergaben sich 
Werte der Gitterkonstante, die sich der gezeichneten Kurve 
völlig einreihten. So ergab die Untersuchung der Legierung 
mit 60°/, Silber z. B. folgende Werte: 


An der Oberfliche . 2... 3,999 AE. 
Nach Abschmirgeln von etwa 0,03 mm 3,996 AE. 
Nach Abschmirgeln von O,lmm. . . 3,988 AE. 


Also schon in einer Tiefe von 0,1 mm stimmte die Abweichung 
mit den übrigen Werten überein. Diese Abweichung für 
20 und 60°/, Silber rührt wahrscheinlich von Verunreinigungen 
an der Oberfläche her, eine Beobachtung, die beim Palladium 
auch schon von I. D. Hanawalt!) gemacht wurde. Die K,- 
Linien der Kupferstrahlung dringen nur sehr wenig in die 
Palladium—Silberlegierungen ein. Eine einfache Rechnung 
über die Abhängigkeit der Intensität von der durchsetzten 
Dicke, ausgeführt für eine Legierung von 20°/, Silber ergibt, 
daß praktisch nur eine Schichtdicke von 0,01 mm für die 
Reflexion in Frage kommt. Da die Dichten der Metalle Silber 
und Palladium beinahe gleich sind, gilt dies, abgesehen von 
minimalen Abweichungen, auch für alle übrigen Mischkristalle 
dieser Reihe. Die verunreinigte Oberflächenschicht, die die 
abweichende Gitterkonstante hatte, braucht also nur etwa 
3—4-10—* cm zu betragen, um ein völlig falsches Bild der 
wahren Gitterkonstante zu ergeben. 

Als diese Messungen schon fertig waren, erschien in der 
Festschrift von G. Siebert G. m. b. H. eine Abhandlung von 
W.Stenzel und J. Weerts?), in der sie bei einigen Palladium— 
Silberlegierungen für die Gitterkonstanten ebenfalls eine Ab- 
weichung im Sinne niedrigerer Werte vom Vegardschen Gesetz 
feststellten. Die größte Abweichung fand sich bei der Legierung 


1) I. D. Hanawalt, Phys. Rev. 33. S. 444. 1929. 
2) W. Stenzel u. J. Weerts, Festschr. z. 50j. Bestehen d. Platin- 
schmelze G. Siebert G. m. b. H., Hanau-Main, 1931, S. 288. 
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nit 52,9°/, 1 adium und betrug 0,007 AE. Die hier gefun- 
denen Abweic., ngen vom Ve gardschen Gesetz stimmen also in 
Richtung und Größe mit den von W. Stenzel und J. Weerts 
gefundenen 'iberein. 


2. Die Gitterkonstanten der wasserstoffbeladenen Palladium- 
Silberlegierungen als Funktion der Wasserstoffkonzentration 
Die Löslichkeit des Wasserstoffes in den Palladium-Silber- 
legierungen nimmt mit zunehmendem Silbergehalt mehr und 
mehr ab, um bei einer Konzentration der Legierungen zwischen 
60— 70°), Silber auf Null zu sinken. Es wurden hier die 
Gitterkonstanten bei verschiedener Wasserstoffbeladung von 
drei Legierungen ge- 
messen, nämlich sol- A 
chen mit 90, 70 und ~ 


mit 60° Palladium. 
Die Legierungen 409\- B 


wurden bei Zimmer- w 
temperatur bei der- Wo 
selben Stromdichte 395 £ 4120 
10 
00 


wie beim reinen Pal- 4 
ladium verschieden 4 
lange mit Wasser- A ‘ 


weilige Konzentration 34 l 
wurde ebenfalls mit- 
2 erhältnis der 90°/,igen Pd- erun 

Die Messungen an der 90°/, Pd-Legierung bei verschie- 
dener Wasserstoffkonzentration uy die Tab. 2 wieder; sie sind 
in Fig. 2 graphisch dargestellt. 


Tabelle 2 
Absorb. H Gitterkonst. Absorb. H Gitterkonst. 
in H/Pd in ÄE in H/Pd in ÄE 
0,00 3,899 0,54 3,998 
018 3,906 0,62 3,999 
, 3,997 0,69 4,002 
0,20 3,998 0,76 4,017 
0,34 3,997 0,80 4,024 
0,45 3,998 


1) Vgl. F.Krüger u. G.Gehm, Ann. d. Phys. [5] = 8. 179. 1933. 
Annalen der Physik. 5. Folge. 16. 
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Auch hier wird, wie beim reinen Palladium, das Gitter 
nur gedehnt, der Gittertypus bleibt unverändert. Man sieht, 
daß auch hier in einem mittleren Intervall zwei Gitter- 
konstanten vorhanden sind und daß die Linien der ersten mit 
zunehmender Wasserstoffbeladung an Intensität schnell ab- 
nehmen, während die der zweiten schnell wachsen. Es existieren 
hier also auch zwei Phasen, denen zwei verschiedene gesättigte 
Lösungen entsprechen. Ein Unterschied gegenüber dem reinen 
Palladium ist insofern vorhanden, als die Gitterdehnung, die 
der wasserstoffarmen Phase entspricht, relativ größer ist als 
beim reinen Palladium, nämlich 0,007, und umgekehrt die 
Dehnung des Gitters der wasserstoffreichen Phase kleiner ist 
als die der entsprechenden Phase beim reinen Palladium, 
nämlich 0,099 ÄE beträgt. Die Gitterkonstanten der beiden 


Ä Phasen sind also hier 
weniger voneinander 
verschieden als beim 
reinen Palladium. Die 


erste Phase ist voll 


“T ausgebildet bei einer 
7.4 Konzentration von 
0 05 10 A etwa c=0,1, die 


Fig. der zweite bei einer Kon- 
Pd-Ag-Legierung als Funktion der H- 0,70 H/Pd; danach 


Konzentration 
tritt dann bei einem 
ziemlich schroffen Knick ein der Wasserstoffbeladung propor- 
tionales Wachsen der Gitterkonstante ein. 
Das Resultat der entsprechenden Messungen an der 
Legierung mit 70°/, Pd und 30°/, Ag gibt die Tab. 3, ihre 
graphische Darstellung Fig. 3. 


Tabelle 3 
Absorb. H Gitterkonst. Absorb. H Gitterkonst. 
in H/Pd in AE in H/Pd in AE 
0,00 3,924 0,57 3,975 
0,13 3,943 0,90 3,977 
0,30 3,943 0,95 3,994 
0,45 3,974 1,08 4,011 


Die auch hier im mittleren Konzentrationsgebiet vor- 
handenen, zwei verschiedenen Phasen entsprechenden Gitter- 
konstanten erscheinen hier noch mehr einander genähert. Die 
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itter | Dehnung der wasserstoffarmen Phase ist gestiegen auf 0,019 AK, 
ieht, | die der wasserstoffreichen ist gesunken auf 0,050 AE. 
ter- Tab.4 enthält das Resultat der Messungen an der wasser- 
mit | stoffbeladenen Legierung mit 60°/, Pd und 40°/, Ag. 
ab- 
ren Tabelle 4 
igte Absorb. H Gitterkonst. Absorb. H Gitterkonst. 
‘die in H/Pd in AE in H Pd in ÄE 
als | 0,00 3,946 0,52 3,976 
die 0,18 3,961 0,64 3,988 
ist 0,36 3,972 0,82 4,002 
Pan. Diese Ergebnisse sind in Fig. 4 graphisch dargestellt. 
ie N Man sieht, daß die Gitterkonstanten der beiden Phasen 
rw sich wieder einander genähert haben; die Dehnung der wasser- 
= stoffarmen Phase betriigt jetzt 0,015, die der wasserstoffreichen 
Die | Phase 0,026 AE. Gleichzeitig ist 
voll das Existenzgebiet beider Phasen a 
= kleiner geworden und auf ein Kon- 
4 zentrationsintervall von 0,2 bis etwa 
0,5H/Pd gesunken. 
on- I 
005 > 
ach 400 
em I 
or- ia 000 Z 
25 Pd 25 30 %Ag 
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Fig. 5. Gitterkonstanten wasser- 
stoffbeladener und unbeladener 
Pd-Ag-Legierungen als Funk- 


.4. Gitterkonstanten der 60°/,igen 
d-Ag-Legierung als Funktion der 
H-Konzentration 


tion der Zusammensetzung 


Um einen Überblick zu gewinnen, wie sich die verschie- 
denen Phasen bei wachsendem Silbergehalt ändern, ist in Fig. 5 
als Abszisse der Silbergehalt in Gewichtsprozenten und als 
Ordinate die Gitterdehnung in Ängströmeinheiten gewählt worden. 
Die Kurven stellen dann die Gitterdehnungen in Ängström- 
einheiten der beiden Phasen dar, und zwar II die der wasser- 
stoffarmen Phase, III die der wasserstoffreichen Phase. Die 


or Gitterkonstanten der Legierungen sind in der Bezugskurve I 
er gleich Null gesetzt. Mit wachsendem Silbergehalt wird die 
ie Differenz zwischen den Gitterkonstanten der wasserstoffarmen 


13* 
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Phase und dem unbeladenen Zustand größer (Kurve II in Fig. 5 
steigt), während die Differenz zwischen dem Werte der wasser- 
stoffreichen Phase und dem unbeladenen Zustand kleiner wird 
(Kurve III fällt. Die Differenz der Gitterdehnungen der beiden 
Phasen sinkt also mit zunehmendem Silbergehalt der Legie- 
rungen. Zwischen 40 und 50°/, Silbergehalt treffen sich die 
beiden Kurven, so daß von da ab nur noch ein einphasiges 
Gebiet existieren kann. 

Bei Beginn der Wasserstoffbeladung bildet sich zunächst 
unter Dehnung des Gitters die erste wasserstoffarme Phase aus, 
bis bei einer bestimmten Konzentration die zweite Phase auf- 
tritt; ist dann die ganze Legierung bis zur Konzentration der 
zweiten Phase gesättigt, so tritt wieder eine Dehnung des 
Gitters auf, und zwar der Wasserstoffbeladung proportional, 
während die erste Dehnung von der unbeladenen Legierung 
bis zur ersten Phase der Wurzel aus dem Druck proportional 
ist. Indem mit steigendem Silbergehalt der Legierung sich 
die Gitterkonstanten der beiden Phasen einander nähern, wird 
gleichzeitig das Konzentrationsintervall, in dem beide Phasen 
koexistent sind, wie besonders ein Vergleich der Fig.2 und 4 
zeigt, immer enger, bis es gleichzeitig mit dem Zusammenfallen 
der beiden Gitterkonstanten gleich Null wird. 


3. Die maximale Wasserstoffabsorption der 
Palladium-Silberlegierungen 

Beim reinen Palladium erfordert die Herstellung des Druck- 
konzentrationsgleichgewichtes bei thermischer Beladung, selbst 
bei Temperaturen von etwa 100°, mehrere Stunden, bei tiefen 
Temperaturen mehrere Tage bei nicht zu dicken Blechen oder 
Drähten; bei elektrolytischer Beladung, die schon bei geringen 
Stromdichten sehr hohen Drucken entspricht (nach Linde und 
Borelius einer Stromdichte von 1 Amp./cm? ein Druck von 
etwa 12000 Atm.), erfolgt die Einstellung des Gleichgewichtes 
schneller. Jedenfalls kann man bei reinem Palladium bei 
hinreichend langsamer Beladung sicher sein, das Gleichgewicht 
erreicht zu haben, weil eben die Diffusionsgeschwindigkeit') 
des Wasserstoffes im Palladium hinreichend groß ist. 

In den Palladium-Silberlegierungen nimmt die Diffusions- 
geschwindigkeit mit zunehmendem Silbergehalt der Legierungen, 
wie A. Coehn und H. Baumgarten?) für die elektrolytische 


1) W. Nernst u. A. Lessing, Nachr. d. Ges. d. Wiss. Göttingen 
S. 146. 1902. 

2) A. Coehn u. H. Baumgarten, Ztschr. f. phys. Chem. 130. 
S. 545. 1927. 
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Beladung bei Zimmertemperatur gezeigt haben, ab und wird 
bei Legierungen mit mehr als 50°/, Silber praktisch gleich Null. 

Um hier eine gleichmäßige Wasserstoff’beladung im ganzen 
beladenen Volumen zu haben, muß man daher in höherem Maße 
die Oberfläche groß im Vergleich zum Volumen halten und die 
Beladungszeiten lang genug wählen. Daß besonders die erstere 
dieser Voraussetzungen nicht immer hinreichend beachtet ist, 
hat die widersprechenden Resultate für die Maximallöslichkeit 
des Wasserstoffes in den Palladium—Silberlegierungen als 
Funktion des Silbergehaltes ergeben, die sich in der Literatur 
finden. 

Bei höheren Temperaturen hat A. Sieverts!) die Aufnahme 
des Wasserstoffs aus der Gasphase in den Palladium—Silber- 
legierungen gemessen. Seine Untersuchungen erstrecken sich 
auf Temperaturen zwischen 138° und 820°; er findet ein aus- 
gesprochenes Maximum der Wasserstoflabsorption bei 40°/, 
Silber bei 270°. Dann beträgt das absorbierte Volumen über 
das Vierfache gegenüber dem beim reinen Palladium. Nach 
Zimmertemperatur extrapoliert wird wahrscheinlich relativ die 
Löslichkeit beim Maximum für Wasserstoff wieder geringer. 

Für die elektrolytische Wasserstoffbeladung haben dagegen 
verschiedene Autoren ein sehr verschiedenes Verhalten der 
maximalen Wasserstofflöslichkeit in den Palladium—Silber- 
legierungen als Funktion des Silbergehalts gefunden: L. Nowack?) 
fand in einer auf Anregung von G. Tammann angestellten 
Untersuchung ein gut ausgeprägtes Maximum der Löslichkeit 
zwischen 20°/, und 30°/, Silbergehalt der Legierung, und 
ähnlich erhielten A. Coehn und H. Jürgens?) ein solches 
Maximum bei etwa 20°/, bis 30°/, Silbergehalt; im Gegensatz 
dazu konstatierten F. Krüger und A.Sacklowski*) und neuer- 
dings wieder J. Schniedermann?) eine nahezu lineare Ab- 
nahme der Löslichkeit des Wasserstoffs in den genannten Legie- 
rungen mit zunehmendem Silbergehalt, jedenfalls keine An- 
deutung eines Maximums. 

Die Wiederholung dieser Absorptionsversuche ergab zu- 
nächst auch eine lineare Abnahme der Löslichkeit mit zu- 
nehmendem Silbergehalt und kein Maximum derselben. Eine 
einfache Berechnung aber zeigt, daß dieser Unterschied in der 
verschiedenen Gestalt der Löslichkeitskurve in Meßfehlern 


1) A. Sieverts, Ztschr. f. anorg. Chem. 92. S. 329. 1915. 
2) L. Nowack, Zitschr. f. anorg. u. allg. Chem. 113. S. 1. 1920. 
3) A. Coehn u. H. Jürgens, Ztschr. f. Phys. 71. S. 179. 1931. 
4) F.Krüger u. A. Sacklowski, a. a. O. auf S. 190 

5) J. Schniedermann, Ann. d. Phys. [5] 13. 8. 761. 1932. 
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nicht begründet sein kann. Dagegen legte die, wie oben er- 
wähnt, mit steigendem Silbergehalt der Legierungen immer 
stärkere Abnahme der Diffusionsgeschwindigkeit des Wasser- 
stoffs in den Legierungen den Verdacht nahe, daß die Ver- 
schiedenheit der Resultate darin begründet sein könnte, daß 
die verschiedenen Autoren Legierungen mit verschiedenem Ver- 
hältnis von Oberfläche zu Volumen benutzten und bei einem 
ungünstigen Wert dieses Verhältnisses vielleicht keine voll- 
ständige Beladung, also kein Gleichgewicht erreicht hatten, 
Auf diese Möglichkeit weist die Tatsache hin, daß auch hier 
zunächst bei den stabförmigen Legierungen mit quadratischem 
Querschnitt von 1 mm Kantenlänge der lineare Abfall der 
Wasserstofflöslichkeit gefunden wurde, daß F. Krüger und 
A.Sacklowski Bleche von mindestens 0,5 mm Dicke und 
J. Schniedermann Drähte von 0,5 mm Querschnitt benutzten, 
während A. Coehn und H. Jürgens Drähte von nur 0,1 mm 
für ihre Versuche verwandt hatten. Es war daher wahr- 
scheinlich, daß sich bei den dickeren Blechen und Drähten 
nur die äußeren Schichten, aber nicht das ganze Volumen mit 
Wasserstoff gesättigt hatte; dann war in diesem Fall natürlich 
nur ein linearer Abfall der Löslichkeit 
mit zunehmendem Silbergehalt zu 
erwarten. 

Um diese Frage zu klären, wur- 
den zunächst die Stäbe der Legie- 
rungen mit 20 und 30°/, Silber zu 
Drähten von 0,1 mm Durchmesser 
ausgezogen und mit diesen die Lös- 
0 =»  lichkeit bei elektrolytischer Beladung 
si er %%# in der in der vorstehenden Arbeit 

Fig. 6. Löslichkeit des geschilderten Weise bestimmt. Die 

a 3 in Pd-Ag- erhaltenen Resultate sind in Fig. 6 

oe eingetragen. Man sieht, obwohl hier 

nur fiir zwei Lösungen in dem in Betracht kommenden Kon- 

zentrationsgebiet die Löslichkeit an den dünnen Drähten be- 

stimmt wurde, daß sich jetzt wie bei L.Nowack und A. Coehn 
und H. Jürgens ein ausgesprochenes Maximum ergibt. 

Hier dürfte also das Lösungsgleichgewicht erreicht sein, 
während das bei den dicken Blechen und Drähten nicht der 
Fall war, weshalb irrtümlich ein linearer Abfall der Löslichkeit 
gefunden war. Eine nähere Untersuchung dieser Verhältnisse 
nicht nur bei den Palladium—Silberlegierungen, sondern auch 
bei den Palladium--Goldlegierungen ist im Gange. 

Da die Bestimmungen der Gitterkonstanten, wie oben er- 
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wähnt, an den dicken Stäben mit quadratischem Querschnitt 
von 1 mm Kantenlänge angestellt waren, mußten die ihnen bei 
der jeweiligen Wasserstoffbehandlung entsprechenden Wasser- 
stoffkonzentrationen korrigiert werden, indem die Konzentration 
jeweils mit dem Verhältnis der Ordinate der beiden Kurven 
in Fig. 6 multipliziert wurde; die so korrigierten Werte sind 
in den Tabellen 2—4 und in den Fig. 2—4 eingetragen. 


4. Leitfähigkeitsmessungen der wasserstoffbeladenen 
Palladium-Silberlegierungen 

Wie in der vorstehenden Arbeit wurden auch hier die 
Leitfähigkeitsmessungen an denselben Stäben wasserstoff- 
beladener Palladium—Silberlegierungen vorgenommen, an dem 
bei der betreffenden Wasserstoffbeladung die Gitterkonstante 
bestimmt, war. So konnten diese Messungen zur Stiitze der 
aus den Anderungen der Gitterkonstanten fiir die Konstitution 
des Wasserstoffs in den 
Legierungen gezogenen 
Schliisse dienen. % 

Die Messungen 
wurden genau so an- 
gestellt, wie die in der 
vorhergehenden Arbeit 
an wasserstoff belade- 
nem Palladium vorge- 
nommenen. Der Quer- 
schnitt der Stäbe der 
Legierungen wurde an 
zehn verschiedenen 100% 50% 100% 
Stellen gemessen und 
daraus das Mittel ge- Fig. 7. Leitfähigkeit der Pd-Ag-Legierungen 
nommen. 

Zunächst wurden die bereits von einer größeren Zahl von 
Forschern gemessenen Leitfähigkeiten der unbeladenen Pal- 
ladium—Silberlegierungen wiederholt. Die Resultate sind in 
Tab. 5 aufgeführt und in Fig. 7 graphisch dargestellt. 

Es ergibt sich für die Leitfähigkeitskurve als Funktion 
des Prozentgehalts der Legierungen die bekannte Kettenlinien- 
form, wie sie für Mischkristalle charakteristisch ist, und wie sie 
für diese Legierungen besonders in den neueren Arbeiten von 
W.Geibel!) und von A. Coehn und H. Jürgens?) gefunden 


en. 


1) W. Geibel, Ztschr. f. anorg. Chem. 70. 8. 240. 1911. 
2) A. Coehn u. H. Jürgens, a. a. 0. 
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Tabelle 5 
Gewichts-°/, | Leitfähigk. Gewichts-°/, | Leitfähigk. 

Silber ‘| x-10-* | nd | Silber "| | Stand 
0 | 8,61 1,16 60 4,05 2,47 
10 3,66 2,73 70 5,98 1,67 
20 | 2,70 3,70 80 9,10 1,10 
30 2,30 4,35 90 15,2 0,66 
40 | 2,20 4,55 100 59,5 0,17 

3,03 3,31 


Das Verhältnis des Widerstandes der maximal mit Wasser- 
stoff beladenen Legierungen zu dem der unbeladenen ist von 
A.Coehn und H. Jürgens gemessen; die von ihnen erhaltenen 
Kurven zeigen nach einem anfänglich steilen Abfall eine fast 
lineare Abnahme mit zunehmendem Silbergehalt der Legie- 
rungen; bei einer Legierung mit 60°/, Silber ist keine Wider- 
standserhöhung mehr nachzuweisen, da hier auch die Wasser- 
stoffaufnahme schon sehr klein geworden ist. 

Hier kam es darauf an, die Änderung des Widerstandes 
bei verschiedenen Wasserstoffkonzentrationen bei derselben 
Legierung zu bestimmen. Das geschah bei einer Legierung 
von 90°/, Pd und 10°/, Ag; die MeBresultate sind in Tab. 6 
wiedergegeben als W/W,-Werte, d.h. die Verhältnisse der 
Widerstände der wasserstoffbeladenen zu den Widerstiinden 
der unbeladenen Legierung; in Kurve C in Fig. 2 sind diese 
Messungen graphisch dargestellt. 


Tabelle 6 
Absorb. H | Leitfähigkeit Widerstand Ww 
in H/Pd “10-4 W, 
0,00 3,66 2,73 1,00 
0,13 3,50 2,86 1,05 
0,20 3,41 2.93 1.08 
0,34 3,35 2.99 1.10 
0,45 3.27 3,06 1.12 
0,54 3,10 3.13 1.15 
0,62 3,13 3,19 117 
0,69 308 | 3% 1.19 
0,76 302 | 330 1.20 
0,80 | 121 


Man sieht, daB diese Kurven der relativen Widerstands- 
änderung in dem Gebiet, in dem nach der Bestimmung der 
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Gitterkonstante zwei Phasen existieren, einen linearen Verlauf 
hat, analog wie die entsprechende Kurve bei dem wasserstoff- 
beladenen Palladium. Die Linearität der Widerstandskurve 
bestätigt also das Vorhandensein zweier Phasen in diesem 
Gebiet. 

Bei den übrigen Legierungen, bei denen die Gitterkonstanten 
gemessen wurden, läßt sich wohl, wie Coehn und Jürgens 
gezeigt haben, die Widerstandsänderung bei maximaler Be- 
ladung bis zu der Legierung zu 60°/, Ag noch nachweisen, 
aber die Abhängigkeit von der Wasserstoffkonzentration wäre 
doch nur schwer genau zu verfolgen gewesen, da der Maximal- 
effekt bei der Legierung mit 70°/, Palladium nach A. Coehn 
und H. Jürgens hier nur noch etwa ein Drittel so groß ist, 
wie bei der Legierung mit 90°, Pailadium. Diese Kurven 
wurden daher nicht aufgenommen. 


Zusammenfassung 


Aus den Messungen der Gitterkonstanten und des elek- 
trischen Widerstandes in Abhängigkeit von der Stärke der 
Wasserstoffbeladung bei den Palladium—Silberlegierungen ver- 
schiedener Zusammensetzung ergab sich, daß sich diese Legie- 
rungen, soweit sie Wasserstoff aufnehmen, d.h. bis zu einem 
Silbergehalt von 60—70°/,, ganz so verhalten wie das Palladium 
selbst bei verschiedenen Wasserstoffbeladungen, d. h., zunächst 
tritt eine Dehnung des Kristallgitters auf, die aber schon bei 
einem geringen Wasserstofigehalt von etwa c = 0,1 H/Pd ab- 
geschlossen ist; bei stärkerer Wasserstoffbeladung beobachtet 
man das gleichzeitige Vorhandensein von zwei Gittern, ein 
weniger und ein stärker mit Wasserstoff beladenes und ge- 
dehntes Gitter. In diesem Gebiet existieren also zwei Phasen 
nebeneinander, die nach Analogie der Verhältnisse beim Pal- 
ladium als zwei verschiedene gesättigte Lösungen von Wasser- 
stoff in den Legierungen aufzufassen sind. Die Kurve der mit 
zunehmender Wasserstoffbeladung auftretenden elektrischen 
Widerstandszunahme, die in diesem Bereich, wie eine Meßreihe 
bei der Legierung bei 90°/, Palladium zeigt, linear verläuft, 
bestätigt diese Auffassung. Mit zunehmendem Silbergehalt der 
Legierungen nähern sich die beiden Gitterkonstanten der beiden 
Phasen mehr und mehr und entsprechend verkürzt sich der 
Konzentrationsbereich, in dem beide koexistent sind. Durch 
Extrapolation ergibt sich, daß von einem Silbergehalt der Legie- 
rungen von mehr als etwa 45°/, Silber an nur noch eine Phase 
vorhanden sein kann. 
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Die Diskrepanz in den Messungen der Löslichkeitskurve 
des Wasserstoffs bei den Palladium—Silberlegierungen in Ab- 
hängigkeit vom Silbergehalt, für die einige Forscher bei elektro- 
lytischer Beladung ein Maximum bei einem Silbergehalt von 
etwa 30°/, Silber, andere dagegen ein lineares Abnehmen mit 
zunehmendem Silbergehalt gefunden hatten, wurde aufgeklärt. 
Es zeigte sich nämlich, daß wegen der mit zunehmendem Silber- 
gehalt der Legierungen stark abnehmenden Diffusionsgeschwindig- 
keit des Wasserstoffs bei dicken Blechen und Drähten wesent- 
lich nur eine Beladung der Oberfliichenschichten eintrat, daß 
sich dagegen bei dünnen Blechen und Drähten eine gleich- 
mäßige Beladung über den ganzen Querschnitt, also ein wirk- 
liches Gleichgewicht ausbildet und dann auch das, von der 
thermischen Beladung bei hohen Temperaturen her bekannte 
Maximum der Löslichkeit des Wasserstofis bei einem Silber- 
gehalt der Legierungen von etwa 40°/, auch bei elektrolytischer 
Beladung bei Zimmertemperatur auftritt. 


Hrn. Dr. E. Nähring sind wir für ständige Unterstützung 
bei den röntgenographischen Aufnahmen zu Dank verpflichtet. 


Greifswald, Physikalisches Institut der Universität, 
22. August 1932. 


(Eingegangen 26. Oktober 1932) 
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Die Beeinflussung einer Telephonmembran durch 
ein Luftpolster mit absorbierender Abschlußwand 


Von K, Schuster und A. Hohberg 
(Mit 6 Figuren) 


1. Fragestellung 


Bildet eine Telephonmembran die eine Grundfläche eines 
Luftzylinders, der im übrigen von vollkommen starren, glatten 
Wänden begrenzt ist und dessen Abmessungen klein gegen- 
über den in Frage kommenden Wellenlängen sind, so liegt 
der Eigenton des Schwingungsgebildes höher als der Eigenton 
der „bloßen“ Telephonmembran. Denn die im abgeschlossenen 
Hohlraum befindliche Luft, das „Luftpolster“, erhöht die 
Elastizität der Telephonmembran. Die „bloße“ Telephon- 
membran soll dabei dadurch charakterisiert sein, daß auf 
ihre strahlende Fläche von der Luft her gar keine Kräfte 
wirken. Dieser Zustand ist beispielsweise nahezu erreicht, 
wenn sich die Telephonmembran in einer starren Wand be- 
findet und in den freien unbegrenzten Halbraum strahlt. Die 
Kräfte, die dann von der Luft her noch wirksam sind, äußern 
sich bekanntlich in einer Vergrößerung der Masse (mit- 
schwingende Mediummasse) und in einer Zunahme der Dämpfung 
(Strahlungsdämpfung). Die gegenüber der idealen bloßen Tele- 
phonmembran hervorgerufene Eigentonerniedrigung und Dämp- 
fungsvergrößerung ist aber sehr gering und kann, falls sie 
doch berücksichtigt werden muß, leicht rechnerisch bestimmt 
werden. 

Wird die der Telephonmembran gegenüberliegende Grund- 
fläche des Luftzylinders nicht mehr von einer starren, glatten, 
sondern von einer teilweise absorbierenden (porösen) Wand 
eingenommen, so ist die gegenüber der bloßen Telephon- 
membran hervorgerufene Eigentonerhöhung geringer als vorher; 
denn die an der Abschlußwand befindlichen Luftteilchen 
werden jetzt nicht mehr zwangsläufig in Ruhe gehalten, die 
Rlastizität des Luftpolsters hat infolgedessen abgenommen. 
Außerdem tritt jetzt gegenüber der bloßen Telephomembran 
eine Dämpfungszunahme auf; denn die durch die Abschluß- 
wand hindurchtretende Schallenergie wird dem Schwingungs- 
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system entzogen. Im folgenden soll nun der Zusammenhang 
zwischen den Absorptionseigenschaften der Abschlußwand und 
der auftretenden Änderung von Eigenton und Dämpfung unter- 
sucht werden, und es soll die Frage geprüft werden, ob sich 
auf diesem Zusammenhang eine Methode zur Bestimmung der 
akustischen Absorptionskoeffizienten aufbauen läßt. 


2. Schallabsorption und akustische Impedanz 


Für den Einfluß der Abschlußwand auf das Schwingungs- 
system ist nicht unmittelbar der Absorptionskoeffizient a maß- 
gebend, sondern die „akustische Impedanz pro Flächeneinheit“ 
3=w-+igqg. Durch a allein ist die Abschlußwand nicht voll- 
kommen charakterisiert; vielmehr muß außer a noch der bei 

einer Schallreflexion auftre- 


tende Phasensprung # ge- 
“88% 0% geben sein. Es ist zweckmäßig, 
sich den Schwingungsvorgang 


vor der Wand in eine rein 
fortschreitende und eine rein 
stehende Welle zerlegt zu 
denken, wokei die Druck- 
knoten der stehenden Welle 
mit den Druckminima des Ge- 
samtvorganges zusammenfallen 
sollen. Durch den Phasen- 
; ; sprung # ist dann die gegen- 
> seitige Lage von stehender 
flektierenden Wand und der Lage Welle und Wand bestimmt. In 
der stehenden Welle Fig. 1 ist schematisch der 
Druckverlauf der stehenden 

Welle angegeben. Ist #= 0, so befindet sich an der Wand ein 
Druckbauch; ist «+ =-+ 2, so befindet sich an der Wand ein 
Druckknoten. Das Wertepaar a, steht mit den beiden Kom- 
ponenten w und q der akustischen Impedanz in einem ein- 
deutigen Zusammenhang. Bezeichnet man die auf Schall- 


widerstand der Luft bezogene Impedanz = +2 


mit } = +17, so ist bekanntlich: 
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Man gelangt zu diesen Formeln, wenn man untersucht, wie 
eine ebene Schallwelle, die senkrecht auf eine Wand von der 
Impedanz z=w-+ig auffällt, reflektiert wird (akustische 
Leitungstheorie). Ist 5* =» — ig der konjugiert komplexe 
Wert von 3, so läßt sich (la) auch in der Form schreiben: 


(2) 1 (1 + 3)-(1 + 3%) 

Die Formeln (1) gestatten eine sehr einfache graphische Dar- 
stellung: Trägt man # und 7 als Koordinaten in einem recht- 
winkligen Koordi- 
natensystem auf, so 
stellen die Kurven 
+ = const diejenigen 
Kreise dar, die durch 
die Punkte A (+ 1,0) 
und B (— 1,0) gehen; 
die Kurven a = const 
bilden ebenfalls eine 
Schar von Kreisen, 
nämlich die ortho- 
gonalen Trajektorien 
zu den Kreisen %# 
= const (Fig. 2). Um 
aus der Figur die 
einem bestimmten 


F 


Fig. 2. Zusammenhan 
Wertepaar #,@ zu- zwischen dem Absorptionskoeffizient a 
geordneten Werte und der akustischen Impedanz 

von a und # zu er- 

halten, muß man den Punkt C (,7) mit den beiden Punkten A 


und B verbinden. Dann ist 


CA\? 
und 1-a= (75) 
Der Punkt (0,0) stellt den vollkommen schallweichen Abschluß 
dar, wie er nahezu bei der Reflexion am offenen Ende einer 
engen Rohrleitung oder bei der Reflexion in Wasser verwirk- 
licht ist; unmittelbar an der abschließenden „Wand“ liegt 
dann ein Druckknoten. Der unendlich ferne Halbkreis ent- 
spricht dem vollkommen schallharten Abschluß; an der Wand 
liegt ein Geschwindigkeitsknoten. Der Punkt (1,0) stellt die 
vollkommen absorbierende Wand dar; vor ihr bildet sich eine 
rein fortschreitende Welle aus. 
Die durch das Wertepaar a, bzw. durch das Werte- 
paar :, 7 charakterisierten Absorptionseigenschaften einer 
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Wand hängen nicht nur von dem Material ab, aus dem die 
absorbierende Oberfläche besteht, sondern von der gesamten 
Anordnung hinter der Oberfläche, soweit sie noch von Schall 
erfüllt wird. Der Absorptionskoeffizient a stellt daher nur 
bei unendlicher Schichtdicke des absorbierenden Mediums 
eine Materialkonstante dar (a). Benutzt man Absorptions- 
platten von endlicher Dicke, so fallen die Verhältnisse an 
der Plattenrückseite um so mehr ins Gewicht, je dünner die 
Platte ist. Es läßt sich leicht zeigen, daß a. durch zwei 
Messungen an einer Platte von endlicher Dicke bestimmt 
werden kann. Es muß dann einmal die Impedanz an der 
Plattenrückseite unendlich groß gemacht werden, d.h. un- 
mittelbar hinter der Platte muß eine starre, glatte Wand an- 
geordnet werden; das andere Mal muß die Impedanz an der 
Plattenrückseite gleich Null sein, d.h. im Abstand 2/4 von 
der Plattenrückseite muß eine starre, glatte Wand an- 
geordnet werden. Aus den in beiden Fällen bestimmten 
Impedanzen 3, und 3, der Plattenvorderseite ergibt sich 3% 
nach der Formel: 30 = Y3,3,. In der komplexen Ebene er- 
gibt sich daher 3. durch eine einfache geometrische Kon- 
struktion: 3. liegt auf der Winkelhalbierenden des Winkels 
zwischen den Vektoren 3, und 3, in der Entfernung J |3,|-|3, | 
vom Nullpunkt. 

Befindet sich eine verhältnismäßig dünne absorbierende 
Schicht unmittelbar vor einer schallharten Wand, so ist an 
der Vorderseite des Absorptionsmaterials die reflektierte Welle 
gegenüber der einfallenden um einen kleinen Betrag verzögert. 


&# liegt daher zwischen den Werten 0 und — T- In Fig. 2 


ist das zwischen # = 0 und & =— + liegende Gebiet, sowei 


es einem unter 70°/, liegenden Absorptionskoeffizienten ent- 
spricht, schraffiert gezeichnet. In diesem Gebiet sind die 
Impedanzen der meisten technischen Absorptionsmaterialien 
enthalten, wenn sie unmittelbar auf eine harte Wand auf- 
gebracht sind. 


3. Die Beeinflussung des Schwingungssystems 


Wie eine absorbierende Wand, die parallel zu einer Membran 
in geringem Abstand von ihr angeordnet ist, den Eigenton und 
die Dämpfung der Membran beeinflußt, läßt sich theoretisch 
am einfachsten im Fall einer Kolbenmembran zeigen. Sind 
die Durchmesser der Membran und der Absorptionsplatte 
groß gegenüber der Wellenlänge, so besteht der Schwingungs- 
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vorgang sowohl im Luftpolster wie im Absorptionsmaterial aus 
ebenen Wellen ganz unabhängig davon, wie Luftpolster und 
Absorptionsmaterial seitlich begrenzt sind. Handelt es sich 
dagegen um eine gegenüber der Wellenlänge kleine Membran, 
so muß nicht nur das Luftpolster, sondern auch das Absorp- 
tionsmaterial der Membrangröße entsprechend seitlich von 
harten Wänden begrenzt werden, wenn die Absorptionsverhält- 
nisse die gleichen wie bei großen Flächen sein sollen. Denn 
ohne die seitliche Begrenzung würde die Schallausbreitung 
im Absorptionsmaterial nicht mehr in der Form von ebenen 
Wellen vor sich gehen. 

Im Fall rein ebener Wellen besteht bekanntlich zwischen 
der Impedanz pro Flächeneinheit 3 an der Abschlußwand und 
der Impedanz pro Flächeneinheit 3, an der Kolbenmembran 


die Beziehung: 
—j,.coskh+isinkh 


(3) en ij-sinkh—coskh 

Dabei ist h die Höhe des Luftzylinders und k= * = °* 
die Wellenzahl. Da die Höhe h klein gegenüber der Wellen- 
länge sein soll (kh < 1), wird 


(4) tix 


wobei kh = x gesetzt ist. 
Durch Einsetzen von (4) in Gl. (2) erhält man 


6) 
oder 


4%, 
a hängt also in der gleichen Weise von 2,4, ab wie von 
w, 5. Es empfiehlt sich daher eine zur Fig. 2 analoge 
graphische Darstellung, bei der 2, und 7, als rechtwinklige 
Koordinaten aufgetragen sind (vgl. Fig. 6). 

Die Komponenten w, und q, der Impedanz an der Kolben- 
membran stehen in unmittelbarer Beziehung zu der experi- 
mentell zu beobachtenden Änderung von Eigenton und Dämpfung 
des Schwingungssystems. Ist & die Elongation der Kolben- 
membran, so stellt 3,-£ den Schalldruck an ihrer Oberfläche 
dar. Infolgedessen ergibt sich als Schwingungsgleichung für die 
unter der Einwirkung einer äußeren Kraft A-e'®' stehende 
Kolbenmembran 


(7) (B+ E=A- 
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wobei W=iam+r+ — s die Impedanz der bloBen Kolben- 
membran, F ihre Fläche ist. Für die natürliche Eigenfrequenz » 
des aus Kolbenmembran und Luftpolster bestehenden Schwin- 


gungssystems gilt daher, wenn die Eigenfrequenz Ys/m der 
bloßen Kolbenmembran mit w, bezeichnet wird: 


== 


@ m m (1 — x + (x 


(8) 


Wird die Dämpfung der bloßen Kolbenmembran „—— = 0, ge 


setzt, so ergibt sich fiir die Dämpfung 0 des gesamten Schwin- 
gungssystems: 
F-co 
(9) 
Die Gleichungen (8) und (9) liefern den gesuchten Zusammen- 
hang zwischen der Beeinflussung des Schwingungssystems und 
den akustischen Eigenschaften der Abschlußwand. Umgekehrt 
kann aus (6) in Kombination mit (8) und (9) der Absorptions- 
koeffizient der Abschlußwand aus Eigentonänderung und 
Dämpfungsänderung sofort ermittelt werden. 

Handelt es sich beispielsweise um eine vollkommen starre, 
glatte Abschlußwand (© = 0, 7=00), so wird nach (8) und (9) 


2 


Für absorbierende Wände ist die Frequenzerhöhung in jedem 
Fall kleiner. Im Fall der völlig.absorbierenden Abschlußwand 
@=1, 7=0) wird 


F.ce, 


w* — w,? = 0; 


Da in diesem Fall von der Kolbenmembran eine rein fort- 
schreitende Welle ausgeht, deren Strahlungswiderstand pro 
Flächeneinheit gleich co ist, ist das angegebene Resultat 
auch direkt vorauszusehen. Die Figg. 3 und 4 zeigen den 
Verlauf der Kurven gleicher Eigentonänderung und gleicher 
Dämpfungsänderung in der w-g-Ebene für den speziellen Fall 
x = 0,05, wie er beispielsweise ungefähr durch h = 0,3 cm und 
v= = = 900 realisiert ist. Die kleinen Pfeile geben die 
Richtung größerer Dämpfungszunahme an. Die Kurven gleicher 
Eigentonänderung (Fig. 3) sind Kreise, deren Mittelpunkte 


auf der g-Achse liegen. Der stark ausgezogene Kreis ent- 
hält diejenigen Punkte, für die die Eigentonänderung ver- 
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schwindet. Innerhalb des Kreises tritt eine Eigenton- 
erniedrigung auf, der Luftzylinder wirkt als mitschwingende 


Mediummasse; außerhalb tritt eine Eigentonerhöhung auf, der 
Luftzylinder wirkt als elastisches Luftpolster. Die Kurven 


Fig. 3. Kurven Fig. 4. Kurven 
gleicher Eigentonänderung gleicher Dämpfungszunahme 


x 

Ss 


gleicher Dämpfungszunahme (Fig. 4) sind diejenigen Kreise, 
die die 7-Achse im Punkte 7 = ı berühren. Für eine dünne 


Absorptionsplatte, die unmittelbar auf einer harten Wand auf- 
gebracht ist, ist im allgemeinen # ziemlich groß. Aus Fig. 4 
geht hervor, daß dann eine starke Dämpfungszunahme der 
Kolbenmembran resultiert. Beispielsweise wird für @ = 15, 
g=— 10, x = 0,05 die Dämpfungszunahme 


F+-co 16 

5’ 

also etwa fünfmal so groß als im Fall einer völlig ab- 
sorbierenden Wand. Man könnte daher bei Telephonen eine 
sehr starke Dämpfung dadurch erzielen, daß man im Gehäuse 
eine dünne Absorptionsplatte anbringt, durch die die Magnet- 
pole hindurchtreten. 

Da sich die Verhältnisse im Fall einer Kolbenmembran 
theoretisch restlos übersehen lassen, was für die eingespannte 
Membran nicht gilt (vgl. unten), würde man am besten auch 
die experimentelle Untersuchung an einem Schallgeber mit 
Kolbenmembran durchführen. Aus rein technischen Gründen 
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wurde aber von uns bei den Messungen nur ein gewöhnliches 
elektromagnetisches Telephon benutzt. Es war also notwendig, 
die Theorie auf den Fall der eingespannten Membran auszu- 
dehnen. Das ist aber in einfacher Weise nur unter gewissen 
einschränkenden Bedingungen möglich. 

Bekanntlich verhält sich eine eingespannte Membran von 
der Fläche F in vielen Fällen wie eine Kolbenmembran von 


K 
der kleineren Fläche a,?- F, wobei a,? = > . (2 n f f (r)rar), 
0 


R der Membranradius, r der Abstand eines beliebigen Membran- 
punktes vom Membranmittelpunkt und f(r) die Durchbiegungs- 
funktion ist. Beispielsweise ist fiir eine eingespannte Membran 
die Eigentonerhöhung durch ein hart abgeschlossenes Luft- 
polster vom Volumen V durch die Gleichung gegeben: 


2 
2 2_ “Ff 2 
— = -a,?F, 


mh 
während sich fiir eine Kolbenmembran die Gleichung ergibt: 


Cor 


Geht von der Membran eine rein fortschreitende Welle aus, 
ist also die Luftkammer durch eine völlig absorbierende Wand 
abgeschlossen, so beträgt die Dämpfungsänderung für eine ein- 


gespannte Membran: 6 — 0, = a ; fir eine Kolben- 


membran beträgt die Dämpfungsänderung dagegen: 
= 


2m 
In beiden Fällen geht die Formel für die eingespannte Membran 
aus der Formel für die Kolbenmembran durch Multiplikation 
mit a,” hervor. Es liegt daher von vornherein nahe, ganz 
allgemein für die Eigentonänderung und Dämpfungsänderung 
einer eingespannten Membran die Gleichungen zu erwarten: 


o m m 

Fce _ 1 (4a,?Fee 
(11) d-%, 2 ( m ) (l—x@? + «a 


Die Konstante mete = b läßt sich aus der Eigentonerhöhung 


bei hartem Abschluß bestimmen. Denn es wird dann: 
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Um die Gl. (10) und (11) abzuleiten, nehmen wir an, daß 
auch jetzt die akustische Impedanz 3 längs der Abschlußwand 
konstant ist und daß die Durchbiegungsform f(r) in dem be- 
trachteten Frequenzbereich frequenzunabhängig ist. Der Schwin- - 
gungszustand in der Luftkammer setzt sich wieder aus einer 
hingehenden und einer zurückkommenden Welle zusammen; 
in jeder Welle treten jetzt auch radiale Schwingungen auf. 
Wir führen Zylinderkoordinaten r,z ein, so daß r den Ab- 
stand eines Luftteilchens von der Symmetrieachse des Luft- 
zylinders, z den Abstand eines Luftteilchens von der Membran 
bedeutet. Die Durchbiegungsform f(r) der Membran soll in eine 
Fourierreihe nach Besselschen Funktionen entwickelt sein: 


0, oo 


f= Da, J,(#,r), wobei die 8, als Wurzeln der Gleichungen 


J, (8, R)= 0 (speziell 8, = 0) bestimmt sind. Das Geschwindig- 
keitspotential in der Luftkammer läßt sich dann in folgender 
Form schreiben: 

0, oo 


(12) a,J,(ß,r)-[A,- + B,.etio:]. 


£ ist die Mittelpunktselongation der Membran; die Konstanten A, 
und B, sind durch die Gleichungen bestimmt: 


ct tayh a+ a, 
= 
v ta, et (3 - — (3+ 
v v 
@@ 
B = . ° 
v ia, en -(s- etiayh , (s+ 
v v 


«, ist gleich Yk?— ß,?, so daB a, =k ist. Durch den Ansatz 

ist auch die Randbedingung erfüllt, daß am Mantel des Luft- 

zylinders keine Radialschwingungen stattfinden können. Die 

m Luftpolster auf die Membran wirkende Kraft ergibt sich 
on zu 


R 
X = 2no- | 
J 0 


(13) 


0, co 
aR? 9g (4, + B,)a,? J,? (8, R). 
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Entspricht die Funktion f(r) einer einfachen Membrandurch- 
biegung ohne Bildung von Knotenlinien, so sind die höheren 
Glieder der Reihe vom dritten Glied ab zu vernachlässigen. 
Unter bestimmten einschränkenden Bedingungen ist auch das 
zweite Glied der Reihe für X klein gegenüber dem ersten: 
1. muß 45 < Tq 3,83 + sein, d. h. das Verhältnis von Höhe 
zu Radius des Luftzylinders darf nicht zu klein werden und 
der Radius muß klein gegenüber der Wellenlänge sein; 2. muß 
7=1 sein. Unter diesen Bedingungen wird die Rückwirkung: 

3+i 

X = Ro - (4, + By) = cg - 

Aus Gl. (14) ergeben sich sofort für die Eigentonverschiebung 
und die Dämpfungsänderung die Gl. (10) und (11). 


4. Meßergebnisse und Diskussion 


Es wurden Versuche mit einer Telephonmembran vom 
Durchmesser 2R = 5,25 cm und von der Dicke D = 0,0195 cm 


ausgeführt. Der Eigenton der bloßen Telephonmembran », = —1 


2a 
lag in der Gegend von 700 Hertz; infolgedessen hatte iz etwa 
den Wert 0,09. Die Dämpfung 6, der Telephonmembran 
betrug etwa 45 [sec=!] entsprechend einer Halbwertsbreite 
von 15 Hertz. Die abnorm geringe Dämpfung war anscheinend 
nur eine Folge davon, daß die Membran sehr sorgfältig ein- 
gespannt war und sehr oft in ihrer Eigenfrequenz angeregt 
worden war. Eigenton und Dämpfung des Schwingungssystems 
wurden aus der Resonanzkurve des Telephons bestimmt, die 
nach der Hahnemann- und Hechtschen Methode in einer 
Wheatstoneschen Brücke aufgenommen wurde. Die genaue 
Größe von », und J, wurde vor jeder Messung aus der Resonanz- 
kurve des in den Halbraum strahlenden Telephons ermittelt. 
Fig. 5 zeigt ein Beispiel eines Kurvenpaares. Die aus der 
Resonanzkurve entnommenen Werte wurden mit Hilfe der 
Rayleighschen Formeln für mitschwingende Mediummasse 
und Strahlungswiderstand einer Kolbenmembran korrigiert. 
Die Konstante b wurde aus einer Messung mit einer dicken 
Messingplatte als Abschlußwand bestimmt. Sie ergab sich zu 


b = 151. Würde man b nach der Gleichung b = a os Soe 
m u-eD 


unter Zugrundelegung der Schwingungsform 


Schu: 


berec 
Luftz 
dann 
grun 
wurd 
kR 
3,83 
her 
Wie 
ursa 
zwisc 
perim 
nicht 
sunge 
wolle 
den ] 
mäßi 
Reso: 
der | 
scher 
men 
Kiger 
und 
eben! 
genai 
Größ 
würd 
sten 
8a 5s 
4a 
= 
Mem 
resul 
ordin 
Kury 
Achs 
Fehl. 
ecke, 
Fehl. 


Schuster u. Hohberg. Beeinfluss. einer Telephonmembran usw. 213 


2 4 

fy=1-2(5) 

berechnen, so würde sich b= 163 ergeben. Die Höhe des . 
Luftzylinders h wurde ungefähr gleich 0,25 cm gewählt. Da 
dann kh etwa gleich 0,03 wurde, war die der Theorie zu- 


grundeliegende Forderung kh<1 gut erfüllt. Tg 3,83 h/R 
wurde dann rund 0,4; die in der Theorie gemachte Annahme 
< 3g 3,88 war da- 
her nur schlecht erfüllt. 
Wie groß die dadurch ver- 
ursachten Diskrepanzen 
zwischen Theorie und Ex- 
periment sind, wurde noch 
nicht näher untersucht. 


Ehe wir auf die Mes- 
sungen selbst eingehen, 
wollen wir die zu erwarten- 
den Meßfehler überschlags- 
mäßig bestimmen. Aus der 
Resonanzkurve, wie sie nach 
der Hahnemann-Hecht- 
schen Methode aufgenom- 
men wird, könnte man den 
Eigenton auf etwa +2Hertz 
und die Halbwertsbreite 
ebenfalls auf etwa + 2 Hertz 
genau entnehmen. Die 


Summe aerin Lisen- 
verbrauchlen und 
der mechanisch- 
akustischen Leistung. 


o* — @? Fig. 5. Resonanzkurven des Telephons 
GréBe = 2(@—,) a) freistrahlend, b) unter dem Einfluß 


@ 
würde dann im ungiinstig- mit 


sten Falle einen Fehler von 
8asec—!, die Dämpfungsänderung ö — 6, einen Fehler von 
4a sec"! aufweisen. Nun ist nach Gl. (10) und (11) 
= 2 Für die Komponenten der 
Membranimpedanz #, und 7, würden dann Fehler von 0,15 
resultieren. Trägt man #, und 7, als rechtwinklige Ko- 
ordinaten auf, so bilden, wie schon gesagt wurde, die 
Kurven a= const Kreise, deren Mittelpunkte auf der =,- 
Achse liegen (Fig. 6). Nimmt man für #, und 7, bestimmte 
Fehlergrenzen an, so sind in der Fig.6 dadurch Recht- 
ecke, deren Seiten parallel zu den Achsen liegen, als 
Fehlerbereiche festgelegt. Aus der Figur geht hervor, daß die 


und q, = — 
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für a resultierenden Fehlergrenzen entsprechend der Dichte 
der Kurven a = const sehr wesentlich von der Lage in der 
%,-Q,-Ebene, d. h. von dem speziellen Absorptionsmaterial, 
abhängen. In der Umgebung des Nullpunktes werden daher 
die Meßfehler ungeheuer groß. Sind dagegen @, und 7, groß 
gegenüber 1, so werden die Fehler minimal. Auch durch die 
Voraussetzungen der Theo- 
rie war die Umgebung des 
Nullpunktes von der Be- 
trachtung ausgeschlossen 
worden. 

Die Eigentoninde- 
rung und die Dämpfungs- 
änderung enthalten außer 
den Ablesefehlern noch 
viel größere Fehler da- 
durch, daß sich während 
einer Meßreihe Eigen- 
ton und Dämpfung des 
Schwingungssystems in un- 
systematischer Weise än- 
dern. Trotz aller Vorsichts- 
maßregeln, wie außer- 
ordentlich kleiner Be- 
triebsstrom, sehr stabile 
Einspannung usw. ließen 
sich derartige Änderungen 
nicht vermeiden. Sie be- 
trugen oft 10 Hertz und 
übertrafen daher die Ab- 


Fig. 6. Lage der Meßpunkte lesefehler beträchtlich. 
in der #% — ü,- Ebene. Die Untersuchungen 
a) Celotex Normal auf starrer Wand, erstreckten sich vor allem 
b) Celotex B auf starrer Wand, auf die drei technischen 
ce) Celotex BB auf starrer Wand Absorptionsmaterialien 


Celotex Normal, CelotexB 
und Celotex BB. Die akustischen Bedingungen an der Rück- 
seite des Absorptionsmaterials wurden variiert: hinter der ab- 
sorbierenden Platte wurde entweder eine Messingplatte, ein 
hartabgeschlossener Luftzylinder von der Höhe 4/4 oder ein 
Schalltrichter angebracht. In diesen 3 Fällen beträgt die 
Impedanz an der Plattenrückseite 3=00,53=0 bzw. j=1. 

Die Tabellen 1—3 geben eine Anzahl Meßresultate für 
die am genauesten untersuchten Fälle Celotex Normal, Celotex B 
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Tabelle 1 
Celotex Normal auf starrer Wand; « = 0,04 


— a 

4,23 8,25 17,7%, 

4,98 8,65 18,1), 

4,70 9,50 15,4°/ 

4,75 8,80 17,1%, 

Tabelle 2 Tabelle 3 
Celotex B auf starrer Wand; Celotex BB auf starrer Wand; 
x = 0,03 x = 0,03 

wy | — % a a 
4,20 4,75 34,6%), 2,70 | 3,55 41°, 
4.20 4,80 34,3%), 2,30 3,25 52°), 
3,70 4,80 32,7%, 255 | 23,50 54°, 
3,60 5,00 31,49), 2,55 | 2,75 50°) 
3,46 5,53 28,3%), 2,60 | 2,85 49°), 
3,34 5,08 30,0°/, 2,75 3,30 44%), 
3,92 5,22 30,4%, 225 | 275 50°), 
3,34 4,95 30,8%), 270 | 3,20 45°), 
4.05 4,55 343°) 250 | 2,65 43%), 
4,18 4,50 35,8% 275 | 3,70 40°), 


und Celotex BB auf harter Wand an. Die Werte für a 
schwanken bei Celotex Normal zwischen 15,4°/, und N ä 
bei Celotex B zwischen 28,3°/, und 35,8°/,, bei Celotex B 
awischen 39,5°/, und 52°/,. Die gemessenen Werte sind in 
die Fig. 6 eingetragen. Sie liegen innerhalb je eines Fehler- 
rechtecks, das unsystematischen Eigentonänderungen bis zu 
etwa 15 Hertz und unsystematischen Dämpfungsänderungen 
bis zu etwa 60 sec! entspricht. 


Zusammenfassung 


In der vorliegenden Arbeit wird zunächst theoretisch unter- 
sucht, in welcher Weise ein Luftpolster mit absorbierender 
Abschlußwand auf eine Kolbenmembran zurückwirkt. Aus der 
Rückwirkung läßt sich der Absorptionskoeffizient der Wand 
berechnen. Unter gewissen einschränkenden Bedingungen lassen 
sich die für die Kolbenmembran abgeleiteten Resultate mit 
einer geringen Modifikation auf eine eingespannte Membran 
übertragen. Bei den Messungen, die an einem elektromagne- 
tischen Telephon ausgeführt wurden, machten sich Eigenton- 
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schwankungen der Telephonmembran, deren Ursache noch nicht 
festgestellt werden konnte, sehr störend bemerkbar. Es er. 
scheint prinzipiell durchaus möglich, auf dem Zusammenhang 
zwischen der Rückwirkung auf eine Membran und den aku- 
stischen Eigenschaften der Abschlußwand eine Methode zur 
Bestimmung des akustischen Absorptionskoeffizienten auf- 
zubauen. Elektromagnetische Telephone erweisen sich jedoch 
für diesen Zweck als ungeeignet. 


Wie uns Hr. Prof. Waetzmann freundlichst mitteilte, 
hat Hr. Dr. Meyer, Berlin ähnliche Untersuchungen mit ver- 
vollkommneten Hilfsmitteln ausgeführt, über die er nächstens 
berichten wird. Die vorliegende Arbeit wurde im Physikalischen 
Institut der Technischen Hochschule Breslau ausgeführt. 
Hrn. Prof. Waetzmann möchten wir auch an dieser Stelle 
für seine zahlreichen Ratschläge und sein förderndes Interesse 
unsern herzlichsten Dank sagen. Der Helmholtz-Gesellschaft 
sind wir für die Bereitstellung experimenteller Hilfsmittel zu 
großem Dank verpflichtet. 


(Eingegangen 28. November 1932) 


G. Je 
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Theorie 
der Leitfähigkeit polarisierbarer Medien. I 


Von George Jaffe 
(Mit 3 Figuren) 


Es fehlt nicht an Versuchen zur Deutung der Anomalien, 
die beim Stromdurchgang durch dielektrische Substanzen — 
feste Körper sowohl, wie Flüssigkeiten — auftreten. Die 
wichtigsten der theoretischen Vorstellungen kann man in zwei 
Gruppen sondern; nach der einen Auffassung handelt es sich 
um Besonderheiten des dielektrischen Verhaltens, nach der 
anderen sind die Erscheinungen auf Bewegung von freien 
Ionen zurückzuführen, so daß man es dann mit Anomalien 
der Leitung zu tun hat. Sehen wir von Maxwells Theorie 
geschichteter Halbleiter ab, so hat die erste Auffassung ihren 
Ausdruck in der Theorie von Pellat!) mit Verallgemeine- 
rungen von v. Schweidler?) und Wagner?) gefunden; neuer- 
dings ist auch die Debyesche Dipolvorstellung herangezogen 
worden.) Wenn nun auch diese Theorien formal die Er- 
scheinungen bis zu einem gewissen Grade beschreiben, so ver- 
sagen sie doch in einem Punkte vollständig. Nach den neueren 
experimentellen Ergebnissen, die Richardson?) und vor allem 
Joff& und seinen Schülern zu danken sind, besteht nämlich 
ein wesentlicher Zug der Erscheinungen in der Ausbildung 
von Polarisationsschichten, die den Elektroden mehr oder 
weniger nahe anliegen, und gerade davon geben die genannten 
Theorien keine Rechenschaft. 

Joffé und seine Mitarbeiter haben daher auch stets ihre 
Ergebnisse vom Standpunkt der Ionentheorie aus diskutiert. 
Bisher ist aber von diesem Gesichtspunkt aus noch keine 


1) H. Pellat, Compt. rend. 128. S. 1312. 1899; Ann. chim. phys. (7) 
18. S. 150. 1899; Journ. Be Phys. (3) 9. S. 313. 1900. 

2) E. v. Schweidler, Ann.d. Phys. 24. S. 711. 1907, daselbst eine 
Zusammenstellung der älteren Literatur. 

3) K. W. Wagner, Ann. d. Phys. 40. S. 817. 1913. 

4) S. Whitehead, Phil. Mag. 9. S. 865. 1930; G. Guében, Phil. 
Mag. 11. S. 405. 1931; 12. S. 197. 1931. 

5) S.W. Richardson, Proc. Roy. Soc. Lond. A.92. S.41, 101.1915; 
107. S. 101. 1925. 
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Theorie aufgestellt worden, die, von konkreten Annahmen aus- 
gehend, zu numerisch prüfbaren Resultaten gelangt; ja, es 
scheint dem Verfasser in der Literatur — soweit sie ihm be- 
kannt ist — noch nicht einmal mit Klarheit ausgesprochen zu 
sein, unter welchen Bedingungen nach der Ionentheorie über- 
haupt Polarisationen von solchem Betrage auftreten können, 
wie sie etwa an Kristallen beobachtet werden. Man beruft 
sich darauf, daß nach der üblichen Ionentheorie, die zwar für 
Gase ausgebildet, aber auch auf andere Medien anwendbar sei, 
Feldverzerrungen auftreten, aber es fehlt der Nachweis, daß 
auf diesem Wege die außerordentliche Abnahme des Auf- 
ladestroms oder gar Rückströme von der beobachteten Größe 
erklärt werden können. 

Diesem Mangel versucht die vorliegende Untersuchung ab- 
zuhelfen, und zwar geht der Verf. von dem Standpunkt aus, 
daß es unmöglich ist, die genannten Erscheinungen durch die 
Bewegung von Ionen zu erklären, die sich in allen Be- 
ziehungen ebenso verhalten, wie die sonst betrachteten. 
Die Anomalien der Leitung werden vielmehr auf die Existenz 
besonderer Ionen zurückgeführt; diese Ionen sollen hin- 
sichtlich ihrer Diffusion, Wanderung im Felde und Wieder- 
vereinigung den üblichen Gesetzen gehorchen, sie sollen aber 
nicht befähigt sein, das Medium an den Grenzen zu verlassen, 
oder ihre Ladungen dort abzugeben. Mathematisch findet diese 
Eigenschaft in einer besonderen Form der Grenzbedingungen 
ihren Ausdruck. Es ist von vornherein einleuchtend, daß so be- 
schaffene Ionen unter der Wirkung eineselektrischen Feldes zwar 
gegen die Elektroden anlaufen, sich aber dort stauen werden, 
so daß Polarisationen von beliebigem Betrage entstehen können, 
die schließlich das Fließen eines Stromes ganz verhindern. Wird 
dann das Feld abgeschaltet, so beginnt eine Strömung der Ionen 
in umgekehrter Richtung als Wirkung ihres Konzentrations- 
gefälles und ihrer gegenseitigen Abstoßung, bis wieder Ruhe 
eingetreten ist. 

Es ist also zu erwarten, daß auf diesem Wege die 
Anomalien der Leitung im Zusammenhang mit der Aus- 
bildung der Polarisationsschichten gedeutet werden können; 
dagegen ist es nicht möglich, durch solche Ionen die Existenz 
eines stationären Stromes zu erklären. Will man also der 
experimentellen Tatsache gerecht werden, daß der „Auflade- 
strom“ nicht auf Null abklingt, sondern in einen stationären 
„Reststrom“ übergeht, so muß man die Existenz von weiteren 
Ionen fordern, die den sonst üblichen Grenzbedingungen ge- 
nügen. Die Auffassung, daß man zur Beschreibung der 
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Anomalien und des Reststroms nicht mit einer Sorte von 
Ionen jeden Zeichens auskommt, ist nicht neu und wird durch 
experimentelle Tatsachen nahegelegt. Beim Studium des Elek- 
trizitätsdurchgangs durch flüssige Dielektrika unter der Wirkung. 
von Radiumstrahlen ist der Verf. schon vor mehr als 20 Jahren 
zu der Auffassung gelangt, daß der stationäre Leitungs- 
vorgang und die Anomalien auf Ionen verschiedener Art zu- 
rückzuführen seien), und diese Auffassung erhielt eine starke 
Stütze, als es später gelang, durch systematische Reinigung 
von Hexan die Anomalien auf ein Minimum herunterzudrücken.?) 
Ebenso scheint nach neuesten Ergebnissen das Auftreten der 
Polarisation in manchen Kristallen an die Existenz von Ver- 
unreinigungen geknüpft. Dafür spricht sehr stark die Tat- 
sache, daß kanadischer Kalkspat, dessen Leitfähigkeit 10000 mal 
kleiner ist als die des isländischen oder russischen Kalkspats, 
im Gegensatz zu den letzteren gar keine Polarisation zeigt), 
ferner die Tatsache, daß es gelingt, im kanadischen Kalkspat 
durch Einführung von Verunreinigungen Polarisationserschei- 
nungen zu erzwingen‘) oder sie in Salpeter auf die gleiche 
Weise erheblich zu verstärken.?) 

Wir überschreiten also kaum die Grenze dessen, was 
empirisch geboten ist, wenn wir zwei Sorten von Trägern 
heranziehen®); hypothetisch — und unseres Wissens neu — 
ist dagegen der Versuch, diese beiden Sorten durch ihr Ver- 
halten an der Grenze zu unterscheiden’), während man früher 
wohl eher an Verschiedenheiten in der Beweglichkeit, Diffusion 
usw. gedacht hat. Natürlich handelt es sich bei der von uns 
eingeführten Unterscheidung letzten Endes auch nur um Ver- 
schiedenheiten in der Geschwindigkeit des Austauschs an den 
Elektroden; aber die Auffassung, daß ein Ion bei der Be- 
rührung mit der Elektrode entweder sofort oder gar nicht 
seine Ladung an diese abgibt®), stellt eine Idealisierung dar, 
die leicht mathematischer Formulierung fähig ist. 


1) G. Jaffé, Ann. d. Phys. 25. S. 257. 1908. 
= G. Jaffé, Ann. d. Phys. 28. S. 326. 1909; man vgl. insbesondere 


8. 352 ff. 
3) A. Joffé, D. Rojanski, K.Sinelnikov, Ztschr. f. Phys. 66. 
8. 143. 1930. 
4) A. Joffé, D. Rojanski, K. Sinelnikow, a. a. 0. S. 164. 
5) B. Hochberg u. V. Joffé, Ztschr. f. Phys. 66. S. 172. 1930. 
6) Auch von A. Smekal (Ztschr. f. Phys. 56. S. 579. 1929) wird 
mit Nachdruck der Standpunkt vertreten, daß die Polarisation auf die 
Anwesenheit von Verunreinigungen zurückzuführen ist. 
7) Näheres darüber ist auf S. 221 ausgeführt. 
8) bzw. selbst in die Elektrode übertritt. 
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Nach alledem wird im folgenden der Versuch gemacht, 
das vollständige Verhalten polarisierbarer Medien bei Strom- 
durchgang mit Hilfe von zwei Sorten von Ionen (beiden 
Zeichens) zu beschreiben. Um jedoch die Rechnung zunächst 
nicht unnötig zu komplizieren, werden wir in den beiden 
ersten Paragraphen vom Reststrom absehen und nur die 
Existenz von solchen Ionen annehmen, die das Medium nicht 
verlassen können; dabei ist § 1 dem einfacheren stationären 
Problem und § 2 dem Studium zeitlich veränderlicher Vor- 
gänge gewidmet. Es wird sich zeigen, daß durch die vor- 
getragene Theorie die Anomalien der Leitung unter der 
Bildung von Polarisationsschichten qualitativ vollständig er- 
faßt werden; auch quantitativ ist die Übereinstimmung in 
manchen Fällen so vollständig, daß sich die Werte der Ionen- 
konstanten daraus ermitteln lassen. 

Im $3 wird dann die Annahme von normalen Leitungs- 
ionen hinzugenommen. Als wichtigstes Resultat ergibt sich 
hier die Tatsache, daß in erster Näherung — und zwar unter 
quantitativ angebbaren Bedingungen — die Anomalien der 
Leitung und der Reststrom sich einfach superponieren; die 
Feldverteilung ist dann die gleiche wie im §1, und der 
Strom berechnet sich auf Grund der so vorgegebenen Feld- 
verteilung. Damit ergibt sich die Aufgabe, die Leitfähigkeit 
eines ionisierten Mediums bei „erzwungener“ Feldverteilung 
zu behandeln; diese Aufgabe wird im § 4 durch sukzessive 
Approximation behandelt und für niedrige und hohe Feldwerte 
gelöst. — Der letzte Paragraph ($ 5) enthält den Vergleich der 
Theorie mit der Erfahrung. 

Unsere Theorie ist phänomenologisch in dem Sinne, daß 
die Gesetze der Diffusion, Bewegung im Felde und Wieder- 
vereinigung als gegeben angesehen werden; die betreffenden 
Konstanten werden als bekannt vorausgesetzt, und die Theorie 
sagt daher auch nichts über ihre Abhängigkeit von der Tem- 
peratur oder anderen Parametern aus. Die zu behandelnden 
partiellen Differentialgleichungen sind nicht linear; hinsicht- 
lich der Methode sind wir daher auf sukzessive Approximation 
angewiesen, nur im §1 ist eine strenge Integration durch- 
geführt. Die Ansätze sind aber so gemacht, daß erforder- 
lichenfalls die Rechnung weitergetrieben werden kann, ohne daß 
prinzipielle mathematische Schwierigkeiten zu erwarten sind. 

Zum Schluß dieser Einleitung möchten wir noch be- 
merken, daß man von unserer Theorie natürlich nicht eine 
Erklärung aller Einzelheiten erwarten darf, die sich an 
polarisierbaren Medien beobachten lassen. Einerseits ist die 
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Methode nur in erster Näherung durchgeführt, andererseits 
handelt es sich bestenfalls um ein stark idealisiertes Modell. 
Es kann kaum fraglich sein, daß außer den hier betrachteten 
Ionenvorgängen auch Inhomogenitäten der Medien („Locker- . 
stellen“ und dergl.) gelegentlich eine wichtige Rolle spielen’), 
und die Vorgänge, bei denen Elektronen beteiligt sind oder 
die bei hohen Werten der Feldstärke dem Durchschlag vor- 
auszugehen pflegen, liegen ganz außerhalb des Rahmens dieser 
Untersuchung. Es war nicht die Absicht des Verf.s, eine 
vollständige und abschließende Theorie aufzustellen, sondern 
vielmehr, einen Weg bis zu Ende zu verfolgen und zu unter- 
suchen, wie weit man auf Grund ganz bestimmt formulierter 
Annahmen in der Deutung des empirischen Materials ge- 
langen kann. 


I. Es gibt nur Ionen erster Art 


In diesem ersten Abschnitt machen wir die Voraus- 
setzung, daß nur je eine Art von positiven und negativen 
Elektrizitätsträgern, die wir „Ionen erster Art“ nennen wollen, 
vorhanden ist. Diese Ionen erster Art sollen bezüglich ihrer 
Wanderung im Felde, ihrer Diffusion und Rekombination den 
üblichen Gesetzen gehorchen, sie sollen sich aber von den 
sonst meistens betrachteten Ionen dadurch unterscheiden, daß 
sie durch die Grenzfläche des betrachteten Mediums hindurch 
keine Ladungen transportieren können, weder in der Form, 
daß sie selbst samt ihren Ladungen durch die Grenzfläche 
hindurchtreten, noch in der Form, daß sie dort ihre Ladung 
abgeben und dadurch aufhören Ionen zu sein. Die Diffe- 
rentialgleichungen des Problems bleiben also die üblichen, 
nur die Grenzbedingungen sind verändert. 

Wir betrachten unter Beschränkung auf das ebene Pro- 
blem eine Schicht des Mediums, die sich von der Ebene z = 0 
bis zur Ebene x =! erstrecken möge; & sei seine Dielektri- 
zitätskonstante, und die Richtung des elektrischen Feldes €, =h 
falle mit der positiven z-Achse zusammen. Wir bezeichnen 
mit p und n die Dichte der positiven bzw. negativen Volum- 
ladung, mit q die pro Zeit- und Volumeinheit erzeugte Elek- 
trizitiitsmenge eines jeden Zeichens; k, k’ und D, D’ stellen 
die spezifischen Geschwindigkeiten und die Diffusionskoeffi- 
zienten der positiven bzw. negativen Ionen dar. Verstehen 
wir endlich noch unter @ den Wiedervereinigungskoeffizienten 


1) Man vgl. dazu A. Smekal, Ztschr. f. techn. Phys. 8. S. 561. 1927 
und die dort zitierten früheren Arbeiten. 
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(bezogen auf Ladungsdichte, nicht auf Ionenzahlen), so lautet 
das System von Differentialgleichungen, wenn wir zeitlich ver- 
änderliche Zustände mit in Betracht ziehen 


0° 
0° 
(2) 4K -An), 
oh 4 
(3) 


Die Randbedingungen sind so zu formulieren, daB sowohl 
der von den positiven, wie der von den negativen Ionen ge- 
tragene Teil des Stroms an beiden Elektroden verschwindet; 
es ist also zu fordern 


dp 


(4) 


x=0 und c=]; 


auBerdem ist die Potentialdifferenz V zwischen den Elektroden 
vorgeschrieben 


(5) = v. 
0 


Handelt es sich um einen zeitlich veränderlichen Vor- 
gang, so muß ein Anfangszustand vorgegeben sein und q undV 
sind als gegebene Zeitfunktionen zu denken. 


§ 1. Stationärer Zustand 


Betrachten wir zunächst zeitlich nicht veränderliche Zu- 
stände und streichen demgemäß die linken Seiten der Diffe- 
rentialgleichungen (1) und (2), so gelangen wir, wenn noch q 
als räumlich konstant vorausgesetzt wird, zu dem gleichen 
System von Differentialgleichungen, das der Verf. früher unter 
anderen Randbedingungen durch sukzessive Approximation be- 
handelt hat.’) Es ist bemerkenswert, daß sich bei den jetzigen 
Randbedingungen die Integration streng durchführen läßt, wenn 
nur die eine Voraussetzung 


(6) 


1) G. Jaffé, Ann. d. Phys. 48. S. 249. 1914; 75. S. 391. 1924. 
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eingeführt wird. Die Vereinfachung des Problems beruht 
darauf, daß der Gesamtstrom 


‚6 
t+hkn+kp), 


der ein erstes Integral der gekiirzten Gleichungen (1) und (2) 
darstellt, als Folge der neuen Randbedingungen (4) im statio- 
nären Zustand verschwindet: die Ionen nehmen unter der 
Wirkung des Feldes nur eine neue Verteilung an, aber ein 
Strom kann nicht fließen. 

Was nun die Einschränkung (6) betrifft, so ist sie eine 
Folge der Townsendschen Relation’), wenn die positiven und 
negativen Ionen gleiche Wertigkeit haben; wollte man aber 
Ionen verschiedener Wertigkeit berücksichtigen, so wäre auch 
das einfache Rekombinationsgesetz aufzugeben. Wir lassen die 
Frage offen, inwieweit die Townsendsche Relation, deren 
Gültigkeit aus einfachen kinetischen Prinzipien folgt, außer auf 
Gase und dielektrische Flüssigkeiten auch auf feste Isolatoren 
anwendbar ist, und beschränken uns in diesem Paragraphen 
auf solche Fälle, in denen (6) zutrifft. Damit ist übrigens der 
Fall keineswegs ausgeschlossen, daß k und D oder k’ und D’ 
gleichzeitig verschwinden, so daß nur eine Ionenart beweglich 
ist. Der allgemeinere Fall, daß (6) nicht zutrifft, läßt sich 
nach dem Näherungsverfahren behandeln, das wir im nächsten 
Paragraphen auf das nicht-stationäre Problem anwenden werden. 

Um die Übersicht über den Einfluß der maßgebenden 
Parameter zu erleichtern, führen wir dimensionslose Variable 
ein und setzen 


(8) 


dadurch nehmen die Differentialgleichungen für den stationären 
Fall, unter TEEN von (6), die folgende Form an 


(9) 


ad? d 
(10) ae + 


i=D 


(11) 


1) Man vgl. etwa Handb. d. Phys. 22. S. 363. 366. 1926. 


tet 

| 
V 
q 
er 
n 
In 
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Um dieses System von Differentialgleichungen zu integri 
machen wir zunächst die Annahme, daß Dati: 
(12) l—-av=0 

ist; das wäre z. B. dann der Fall, wenn q und a verschwinden, d.h. 
wenn weder Bildung noch Rekombination von Ionen stattfindet. Von 
dieser Einschränkung werden wir uns hinterher wieder befreien. Unter 
der Annahme (12) lassen sich die Gl. (9) und (10) einmal integrieren und 
ergeben dann bei Berücksichtigung der Randbedingung (4) 


(13) 
d 
(14) 


Aus (11), (13) und (14) lassen sich leicht die folgenden Gleichungen 
herleiten 


i 

d 1 dn 

(16) iat a): 
von denen die zweite integriert 

7 
(17) aa 95" +7 
liefert, wenn die Integrationskonstante mit c,/4 bezeichnet wird. Er- 
weitert man (17) mit 7 at , 80 kann man nochmals integrieren und er- 


hält unter Einführung einer weiteren Konstante c,/8 


d 1 — 
bzw. 
+c 
& bestimmt sich also als elliptisches Integral in Funktion von 7, oder 
letzteres als elliptische Funktion von & Denkt man sich diese aus- 


gewertet, so ergeben sich x und » aus der Kombination von (15) mit (17) 
und (11) mit (18): 


1 
(20) +$ +Vi+ar+a), 
1 2 Cy 


Bisher haben wir an der Annahme festgehalten, daß weder Bildung 
noch Rekombination von Ionen stattfindet; unter dieser Voraussetzung 
wären die Gesamtzahlen der positiven und negativen Ionen als gegeben 
anzusehen, und aus diesen beiden Daten zusammen mit der Forderung (5) 
bestimmen sich dann die drei Integrationskonstanten c, bis c,.') 


1) Wir gehen auf die Konstantenbestimmung für diesen Fall nicht 
ein, jedoch ist die Rechnung im nächsten Paragraphen in erster Näherung 
durchgeführt [Formeln (84), (84a) und (88)]. 


Von d 
kalisel 


(27) 


Man e 


handeı 
Wurze 


Stelle 
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wir at 

D 


richtu: 
nicht | 
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Nun überzeugen wir uns aber hinterher, daß bei unserer Lösung das 
Produkt 1 » » konstant herauskommt; in der Tat folgt aus (20) und (21) 


und dieselbe Tatsache ergibt sich auch, wenn man (13) und (14) un- 


mittelbar in der Form 
(23) 


integriert.‘) Unterwirft man also die Konstanten c, und ¢, der Be- 
dingung 

1 
16 | 4 «) 
so wird damit der Gl.(12) genügt und wir haben damit in (19), (20) 
und (21) die vollständige Lösung des Systems (9), (10) und (11) vor uns. 
Die Konstantenbestimmung verläuft aber jetzt anders, als in dem oben 
skizzierten Fall. Findet nämlich Bildung und Wiedervereinigung von 
Ionen statt, so sind nicht mehr die Gesamtzahlen der positiven und 
negativen Ionen als bekannt anzusehen, sie bestimmen sich vielmehr 
erst aus der Lösung der Aufgabe. Dagegen ist leicht einzusehen, dab 
dann die Differenz der insgesamt vorhandenen positiven und negativen 
Ladungen vorgegeben werden muß; denn denkt man sich den stationären 
Zustand aus irgendeinem Anfangszustand hervorgegangen, so bleibt der 
ursprünglich vorhandene LadungsiiberschuB des einen Zeichens be- 
stehen, weil immer gleich viel positive und negative Ladungen ent- 
stehen und vergehen (und ein Verlust an den Elektroden nicht statt- 
findet!). Diese vorzugebende Differenz der positiven und negativen Ge- 
samtladung liefert zu (24) eine zweite Bestimmungsgleichung für die 
Konstanten c, bis c,, während die dritte wieder in (5) gegeben ist. 

Wir schreiten nunmehr zur Vervollständigung unserer Lösung 
durch die Konstantenbestimmung, wollen uns aber von vornherein auf 
den einfacheren Fall beschränken, daß kein Überschuß von Ionen des 
einen Zeichens vorhanden ist. Bedeute 7, eine Wurzel des Polynoms 
in (18), so kann man 


(25) Cy = — +6 N”) 
setzen und findet dann aus (24) und (25) 
(26) — 2n,°. 


Von diesen beiden Werten führt nur das obere Vorzeichen zu physi- 


‚ kalisch brauchbaren Lösungen. Daher läßt sich (18) in die Form bringen 


(2%) = + 8 — 

Man erkennt, daß fiir 4,2 <8 nur die eine positiv reelle Wurzel 7, vor- 
handen ist, im Falle 7,? > 8 noch die zweite Yr,? — 8; dieser zweiten 
Wurzel entspricht aber kein reeller Wert von & Es gibt also nur eine 
Stelle &= &,, bei der 7=n, positiv ist und dy/d& verschwindet.) 


1) Die Konstanten ce, und c, sind natürlich nicht mehr unabhängig, 
sondern ergeben sich gemäß (20), (21) aus c,, c, und c,; darum haben 
wir auch die Integration von n vorangestellt. 

__ 2) Das Wurzelzeichen in (27) ist bei der vorgeschriebenen Feld- 
richtung positiv zu verstehen, weil man sonst der Randbedingung (5) 
nicht genügen kann. 
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Nun gehen die Gl. (11), (13), (14) ineinander über, wenn man & durch -£ 
ersetzt und a mit » vertauscht; daraus folgt, daß 7 bezüglich der 
Stelle § = &, symmetrisch verläuft und daß a und » durch Spiegelung 
an der Geraden & = &, auseinander hervorgehen. 

Diese Eigenschaften gelten noch allgemein, nimmt man aber die 
Forderung hinzu, daß die Gesamtladung 


l 
é 
(28) Q= [w- mar - (ho - 20) 


verschwindet, so muß h(l) = h (0) und daher auch 7 an beiden Elektroden 
gleich sein. Das kann — für beliebige Werte von ! — nur dann zu- 
treffen, wenn & mit der Mitte des Elektrodenabstands übereinstimmt; 
ist andererseits ( von 0 verschieden, so ist § aus der Forderung zu 
bestimmen, daß @ einen vorgegebenen Wert hat. Da wir den letzteren 
Fall ausschließen, wollen wir für den Rest dieses Paragraphen den 


0-Punkt unseres Koordinatensystems nach x = a verlegen. 


Für die weiteren Rechnungen sind die beiden Fälle zu unter- 
scheiden, daß 7,” kleiner oder größer als 8 ist. Wir behandeln zunächst 
den Fall 7,2=5. Durch die Substitutionen 


2 
(29) = (089; 3- = x? 
geht (27) integriert in 


dq 

30 
V 1 — x' sin* 


über; dabei ist über die Integrationskonstante bereits so verfügt, dab 
der Minimalwert 7 = „, an der Stelle & = 0 liegt. 
Die Beziehung (30) läßt sich nunmehr umkehren 

31 $338 

und das so gefundene n ist zur Bestimmung der letzten Kon- 
stanten nach £ zu integrieren. Das unbestimmte Integral läßt 
sich ausführen 


& 


32 
(82) dnV2E— nsny2E 
wobei 
2 
(33) x? +, =8 


gesetzt ist. Bezeichnet man die zu / gehörige „reduzierte 
Dicke“ mit 0 


die | 
aus 

weil 
sein 


(37) 


stitut 
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so erfordert die Bedingung (5) wegen der Symmetrie von 7 
unter Berücksichtigung von (8) 


(35) 


sn 
2 
» Vv’, miss, 
dn — — x, —— 
2 y2 


die den Zusammenhang zwischen ö und y, darstellt. 

Die Ionendichten erhält man jetzt aus (20), (21) oder 
aus (23) und (32); in letzterem Falle ist c, =c, = 1 zu setzen, 
weil nach (12) a-»=1 ist und außerdem fiir §=0 a, =», 
sein muß. So ergibt sich 


dny2E— x,sny2E 


Wir wenden uns nunmehr dem Fall 7?>=S zu. Durch die Sub- 
stitutionen 


(37) 
v 


0 
und Integration findet man aus (27) 
(39) = K mE 3 
V1 —#'sin® 2 


die Integrationskonstante ist wieder so bestimmt, daß 7 fiir § = & den 
Wert 7, annimmt, und K ist das vollständige elliptische Integral erster 
Gattung. Durch Umkehrung der Beziehung (39) und analoge Rech- 
nungen wie oben findet man statt (31), (32), (36), (37) die entsprechenden 
Gleichungen 


(31’) 


der 
lung 
6/2 
0 
Aus (32) und (35) erhilt man die Beziehung | 
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(37) 


Im Grenzfalle 7,? = 8 gehen die beiden Lösungssysteme ineinander 
über; dann ist x = 0, x, = 1 und man erhält 


No 
(31”) = —— ny? = 8; 
7 cos V2 & ’ 
der zugehörige Wert von ö bestimmt sich aus 
Vv’ 


(36”) si n? = 8. 


V2 e” 41’ 

Unsere Formeln lassen erkennen, daß die Lösung außer 
von der unabhängigen Variablen & nur von den beiden Para- 
metern ö und V’, d.h. von der „reduzierten Dicke“ (34) und 
einer „reduzierten Potentialdifferenz“ [vgl. (35)] abhängt. Für 
die numerische Berechnung muß bei vorgegebenem ö und V’ 
zunächst der Minimalwert der Feldstärke 7, errechnet werden. 
Dafür stehen die transzendenten Gl. (36) und (36’) zur Ver- 
fügung; es muß aber betont werden, daß diese Gleichungen 
sich nicht nach 7, auflösen lassen, weil die Funktionen sn 
und dn den Wert von 7, noch implizit über den Modul x 
[vgl. (29) bzw. (38)] enthalten. Will man also die streng 
gültigen Formeln auswerten, so muß man zuerst bei ge- 
gebenem V’ aus (36) bzw. (36°) ö als Funktion von 2, be- 
rechnen und diese Beziehung graphisch umkehren. 

Wir gehen auf diese Rechnungen nicht ein, sondern be- 
schränken uns auf denjenigen Grenzfall, der das größte 
physikalische Interesse bietet. Bei fast allen Messungen über 
polarisierbare Medien liegen die numerischen Verhältnisse so, 
daß die Gebiete freier Raumladungen auf relativ dünne, den 
Elektroden anliegende Schichten beschränkt sind. Dieser 
Grenzfall tritt dann ein, wenn der Parameter ö eine große 
Zahl ist, und das drückt sich in unseren Formeln so aus, 
daß x, ->0 und «<-> 1 geht. Dann erhält man nach be- 
kannten Grenzbeziehungen aus (31) und (37) 
V2E 
1- x,? Sin? V2E 
1+x Sin V2 


41 


(40) n= 
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während x, und 7, sich gemäß (36) aus 


berechnen.) 
Für große Werte von ö — und zwar schon für Werte 


ö ‘ 
von etwa ——- = 5 an — wird also x, sehr klein, und dann 


sind die Formeln (40), 
(41) gültig. Was die 7 
Abhingigkeit der Funk- 0-40 
tionen 7, a und » von p22 
V’ anbetrifft, so er- 
kennt man aus (40) und ; 
(41), daß sie sich schon 75 2-09 
für mäßige Werte von 
V’ einem Grenzverlauf 
annähern, da 4, und x, 
dann nur noch sehr 
wenig von V’’ abhängen. 


In Figg. 1 und 2 
ist der Verlauf von 7 
und a für den Para- ) 


meterwert ö=Y2-5und 


/ v:04 
einige Werte von V’ 4 7 
gezeichnet; als Abszisse 7, 

ist Y2& gewählt, sc 

daß die Elektroden bei 7 Z 

Y2£=+5 liegen. Von mez 

der symmetrischen 7- 7 2 3 4 5 
Kurve ist nur der Zweig 125 


: > . — und Fig. 1. Die reduzierte Feldstärke 7 als Funktion 
er Verlauf von » er- yon V2 & für verschiedene Werte von V’ 


1) Es ist erforderlich, die im Text gegebene zweite Näherung zu 
benutzen, wenn die Formeln für beliebige V’ gelten sollen; begnügt 
man sich mit der Näherung n = Coj V2&, so muß x, = — 
klein sein. Durch Einsetzen der Näherungslösungen (40), (41) in die 
Differentialgleichungen überzeugt man sich, daß der Fehler von der 
Ordnung n,? Coj? Y2 € ist; die Niherung versagt also nur in unmittel- 
barer Nachbarschaft der Elektroden für große V’. 
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gibt sich aus demjenigen von a durch Spiegelung an der Ordi- 
natenachse. 

Denkt man sich die drei Kurven für 7, a und » bei ge- 
gebenem ö-Wert auf dasselbe Blatt gezeichnet, so gelten diese 
Kurven auch noch für den Fall, daß die Gesamtladung nicht 


Z v:05 
v:022 
ea 


u 
4 


% 


KW 


4 5 


0 7 2 J 


Fig. 2. Die reduzierte Ionendichte x als Funktion von V2 & 
fiir verschiedene Werte von V’ 


verschwindet, aber die Lage der Elektroden ist dann noch zu 
bestimmen‘); man muß ein Intervall von der Breite ¥2-d so 
lange nach links oder rechts verschieben, bis der vorgeschriebene 
Wert von Q wiedergegeben wird. 

Der Grenzfall x, -> 0 und x—» 1 tritt in den Formeln (31), 
(36’) und (37’) ein, wenn 8/V’ eine kleine Zahl ist; es ergibt 
sich dann 


(40') 
(41’) %= = ens, 

ö 
(42) = 


1) Das folgt daraus, daß man dann lediglich in den Formeln (30) 
bis (42) & durch & - &, zu ersetzen braucht; &, bestimmt sich aus der 
schon auf $. 226 angegebenen Forderung. 
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Das Feld ist näherungsweise homogen, die Näherung gilt 
aber nur dann im ganzen Intervall 0 = §= 7 wenn außer as 
auch V’ eine kleine Zahl ist. 


§ 2. Nichtstationäre Zustände 


Bei der Behandlung zeitlich veränderlicher Vorgänge sind 
wir auf sukzessive Approximation angewiesen. Einen geeig- 
neten Ausgangspunkt für den Ansatz gewinnen wir, wenn wir 
Bildung und Wiedervereinigung der Ionen vernachlässigen und 
dem Feld den Wert h,= V/l erteilen. In erster Näherung 
werden also die Gesamtzahlen der positiven und negativen 
Ionen als gegeben angesehen, und der Vorgang wird so be- 
handelt, als ob die Ionen sich in einem homogenen Felde be- 
wegten; die Rückwirkung der Ionenverteilung auf das Feld 
wird dann erst hinterher berechnet. Beides erreichen wir, 
indem wir die Unbekannten nach Potenzen von q": entwickeln 
und folgenden Ansatz machen: 

(43) neng"+tng 
+h, q+ +... 

Die Grenzbedingungen sind dann dahin zu erweitern, 
daß die Gl. (4) für jede Näherung p,, n,, h, gelten müssen, 
während das Integral (5) für alle h, ü=1, 2, 3...) ver- 
schwinden muß. 

Geht man mit dem Ansatz (43) in die Gl. (1), (2), (3) 
hinein), und vergleichen die Koeffizienten gleich hoher Po- 
tenzen von q's, so erhalten wir für die erste Näherung das 
System von Differentialgleichungen: 


Op 0? OP, _ 
On, , On, vp, Om _ 
(45) kh ax 0, 
ah 4 
6) 


1) Wir sehen von der Einführung der reduzierten Variablen (8) 
ab, weil diese erst bei Mitberiicksichtigung von Bildung und Wieder- 
vereinigung eine Vereinfachung herbeiführen. Substituiert man aber (8) 
in das System (1), (2), (8), ohne von (6) Gebrauch zu machen, und wählt 
!=Yget als reduzierte Zeit, so erkennt man, daß x, », 7 nur von &,t’ 
und folgenden Parametern abhängen: der reduzierten Dicke 5, dem 
reduzierten Potential V’ (Gl. (34) und (35)], der Langevinschen Zahl 
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p, und n, berechnen sich also in dieser Näherung unabhängig 

voneinander, während sich h, hinterher durch eine Quadratur 

bestimmt. Die Gültigkeit der Relation (6) wird nicht mehr 
vorausgesetzt. 


Die Gl. (44) und (45) lassen sich nach Methoden behandeln, die 
denen der Wärmeleitungstheorie nachgebildet sind, und zwar läßt sich 
sowohl die Methode der ,,Wirmepole“, wie diejenige der „Eigenlösungen“ 
übertragen. Um zu Formeln zu gelangen, die für alle Werte der Para- 
meter eine bequeme Berechnung gestatten, empfiehlt es sich, beide Be- 
rechnungsarten durchzuführen. 

Zunächst verifiziert man leicht, daß 


(e— i— 
V4at V4ar 
Lösungen der Gl. (44) bzw. (45) sind, wenn 
Dt=t, D’i=7, 
(48) > _H’ 


gesetzt wird. Die Lösungen (47) entsprechen einem „Wärmepol“ der 
Wärmeleitungstheorie mit dem Unterschied, daß sich der Mittelpunkt 
der Verteilungen (47) mit der Geschwindigkeit kh, bzw. k’h, nach der 
Seite der positiven bzw. negativen x verschiebt. 

Wir behandeln zunächst nur die Gl. (44). Ohne Rücksicht auf die 
Randbedingungen ist 


+00 


Hr? 

1 

(49) = [ (iy di 
Pr ane f 


—co 
eine Lösung von (44), die sich für ¢ = 0 auf p, = f(x) reduziert. 

Den Randbedingungen geniigen wir nunmehr durch Anwendung 
des Prinzips der Spiegelung. Unsere Randbedingung (4) ist analog zu 
behandeln, wie in der Wärmeleitungstheorie die Grenzbedingung, dab 
Abkühlung nach außen stattfindet. Nehmen wir vorläufig an, daß der 
Grenzbedingung nur an der einen Ebene x = (0 zu genügen sei, s0 
findet man durch völlig analoge Rechnungen wie in der Wärmetheorie') 
folgendes: 


0 co 
50) p= we +f rae 
— oo 0 


ist eine Lösung von (44), die sich für ¢ = 0 auf 


(51) mn=f@, 0zr<o 


1) Man vgl. etwa Riemann-Weber, Die partiellen Differential- 
gleichungen der mathematischen Physik. 5. Aufl. II. Bd. $ 39. 
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reduziert und die der Randbedingung 


dp, 
genügt, wenn F, (4) aus der Differentialgleichung 
(53) =e * (a) — HF 
bestimmt und die Integrationskonstante durch die Forderung 
(54) F, (— 0) =f) 


festgelegt wird. So ergibt sich 

(55) Fadi), 0si<o. 
0 


Die genannten Eigenschaften der Lésung (50) sind leicht direkt zu 
verifizieren. 

In völlig analoger Weise ergibt sich, durch Spiegelung an x = |, 
daß die Lösung 


ı 
(6) py = — [ Rove dit, 
Ant 
— oo 


die sich für ¢ = 0 auf 
(57) -x<rsl 
reduziert, der Grenzbedingung 


R dp, 
(58) _,=0 


genügt, wenn F, (4) aus der Differentialgleichung 

dF, _ 

HF,(j) =e di + Hg (21-2) 
mit der Bedingung 

(60) gl—0 =F, +0) 

bestimmt wird. Das fiihrt auf 


(59) 


(61) F,@ =e"@-) (ger-m+H[gel- mar), Isi<om. 


Um nun den Randbedingungen an den beiden Grenzen x = 0 und 
x = zugleich zu genügen, müssen wir ein Verfahren anwenden, das 
auf unendlich viele Spiegelungen der eben dargelegten Art — ab- 
wechselnd an x = 0 und r = 1 — hinausläuft. 

Der Anfangszustand sei in der Weise gegeben, daß für ¢ = 0 
(62) 0<rsl, 
vorgeschrieben ist. Wir versuchen nun zwei Funktionen 
(63) F,@, Fw, Isı<e, 
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so zu bestimmen, daß 
ı 


If... 

di + [rwe di 
nt 

(64) 


co 


+ F F, (We 4: au} 


l 
den beiden Grenzbedingungen (52) und (58) geniigt. 

Es ist fiir das Folgende zweckmäßig, die Funktion f(A) je nach 
Bedarf zu F, (A) oder zu F, (4) zu rechnen, und diese letzteren Funk- 
tionen folgendermaßen zu unterteilen: 


F,ü)= —(m— 11, m=1,2,3---. 
66) Osist, 
mlsSis(m+1)l, m=1,2,3--- 


Fassen wir die so definierte Funktion F, (A) als bekannt auf, so 
liefert die Grenzbedingung (52) gemäß (53) die Forderung 


67) dF(-b HR), 


da 


und sieht man umgekehrt die Funktion (65) als bekannt an, so erfordert 
die zweite Grenzbedingung (58) nach (59) die Giiltigkeit von 


Eliminiert man zwischen (67) und (68) F, (4), so ergibt sich 


+ HF,(-) = | 


(68) 


(69) + HF,(—A’ ime aF, + HF,21- 3); 
di di 

diese Gleichung kann man als Differentialgleichung für die Funktion 

(70) wa F 21-2 


auffassen. Integriert man sie, so resultiert für F, (4) die folgende 
Periodizitätseigenschaft 


(71) F,(—@ +20) =e" (- + 
dabei ist C, eine noch verfiigbare Integrationskonstante. 

In véllig analoger Weise findet man 
(72) F, =e”! Fa 2) + 


Die beiden Funktionalrelationen (71) und (72) sind ee gleich- 
bedeutend (wie man erkennt, wenn man in (71) die Substitution 


”=— (442) ausführt, so daß man F, (4) und F, (a) als einheitliche 
Funktion auffassen kann, die im ganzen Bereich — © <i< +@ dem 
gleichen Bildungsgesetz genügt. 

Durch die beiden Beziehungen (71), (72) kann man nun die F,” 
und F,” rekurrent berechnen, wenn man den Funktionsverlauf für ein 


(75) 


G. Je 
Inter 

1 

an d 
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mitte 

zeich 

(73) 

ein. 

(74) 

4 
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im 
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(76) 

fach 

Man 
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Intervall von der Breite 2! kennt; die jeweils auftretenden Konstanten 
6,” und C,” sind dabei so zu bestimmen, daß die Fortsetzung sich stetig 
an das vorhergehende Intervall anschließt. Als Ausgangsintervall ver- 
wendet man 1=1=1 oder 0=1= 21; dabei bestimmt sich F,'(d) un- 
mittelbar aus (55), oder F,1(4) aus (61). ‘ 

Wir wählen 0 —> 27 als Grundintervall und führen die Be- 
zeichnung: 

21-4 


ı 


ein. Setzen wir außerdem 


ı 
(74) A= [fway, 
v 
so läßt sich das Ergebnis der Rechnung folgendermaßen darstellen: 
m—1 
(a) =e + 2m) — AHe"* >} 
0 
-2mlzis-2m-Nl, 
m—1 


0 
- 


0 
2ml=A= (2m + 
m—1 


0 


m=1, 2,3--- 


Denkt man sich (73) und (75) in (64) eingesetzt, so haben wir die 
allgemeine Lösung der Differentialgleichung (44) vor uns, die sich für 
t=(0 im Intervall O<A=1 auf f(x) reduziert und den Grenzbedin- 
gungen genügt. Ein Blick auf die Gl. (44), (45), (47) und (48) zeigt, 

ß man die analoge Lösung von (45) erhält, wenn man r mit 7’ und 
H mit — H’ vertauscht. 

Wir wenden diese allgemeinen Formeln auf den einfachsten Fall 

an, dab für ¢ = 0 p konstant ist. 


(76) t= 0, 0=si=sl. 


In diesem Falle stellt sich das Ergebnis der Rechnung etwas ein- 
facher dar, wenn man das Grundintervall — !<=i = +1 zugrunde legt. 
Man findet aus den allgemeinen Formeln (73), (75) oder unmittelbar aus 
(55), (71), (72): 


*h 
k- 
rt 
(75) 
| 


(77) 


F, (4) und F, (A) als eine Funktion zusammenzufassen und 
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m—1 
F?™+! = Hl + ml) fl + Hü+2m 
0 
—2m+bisis-2ml, 
—1 
2m Hi Hl 
F, (4) = —CHle 
v 
(a) = Ce"*11 + HA 
-l=ijéz=0, 
(4) = 
m—1 
v 


m=1, 2, 


Für die Durchführung der Integrationen in (64) —— es a 
olgendermaßen 


in vier Teile zu zerlegen: 


(78) < 


= 
0, 1, 2, 
m 
@,” (4) = 
2m-bisi=s2ml, m=1, 2, 3--- 
m=O, 1, 2, 3--- 


(r—4— 


Multipliziert man jede dieser Teilfunktionen mit : 


integriert über ihren Gültigkeitsbereich und summiert dann über die 
zulässigen Werte von m, so gelangt man zu dem Schlußresultat 


4 
(79) = >: 
1 


G. Je 


mit 


(80) 


(81) 


und 
die { 


(82) 


(83) 


und 


2mlı, 


en 
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mit der Bedeutung: 


+00 
Dn 
—oo 
+00 ie 
-00 7 
(80) | 
—oo 
oo 
V4r 
Dabei bedeutet ® (x) das Gauss sche Fehlerintegral'): 
(81) fe "as 
Va 
0 


und 4, bis 4,, sowie das gleich benötigte 4,, sind Abkürzungen für 
die folgenden Differenzen: 


Y4r V4r 
(82) Yır 
(mi+z+Hr% 
4, 4r 4r 
(@m+Di-z+ Hn 


Fiir das Folgende notieren wir noch: 


+00 
bs - fins, + 4, 
Ou 2 Var 
— oo 
ea 1 
Hml 
m € H + —— 
2 2 { 4 Var 


1) Man vgl. E. Jahnke u. F. Emde, Funktionentafeln. Leipzig 
und Berlin, 1909, S. 31 ff. 


(83) 


| 

| 
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Man kann verifizieren, daß (79) der Differentialgleichung (44) 
genügt; ebenso erkennt man leicht, daß (83) den Grenzbedin- 
gungen (52) und (58) genügt. Für {= 0 verschwinden alle 
Glieder mit Ausnahme des Gliedes m = 0 in P,, und dieses 
liefert p, = C, wie es sein soll. Für t=oo werden alle 
Glieder zu 0, außer denen der Summe P,; diese liefern 


(S4) = t=O. 


Die Lösung von (45) erhält man aus (79) durch Vertauschung 
von C, H, r mit C’,— H’, r’, wenn C’ die konstante Anfangs- 
konzentration der negativen Ladung ist. Für t =o bleiben 
dann nur die Glieder in P, übrig und liefern 

Die Endwerte (84) und (85) sind natürlich in Überein- 
stimmung mit der strengen Lösung (23), wenn man in dieser 
das Feld h als konstant behandelt. 

Unsere Grenzbedingungen verhindern im stationären Zu- 
stand das Fließen eines Stromes (vgl. S. 223); solange aber 
Bewegungen der Ionen im Inneren des Dielektrikums statt- 
finden, werden Ladungen auf den Elektroden influenziert; die 
Änderungen dieser Ladungen sind einem Strom äquivalent, 
den wir nun noch berechnen wollen, weil er am leichtesten 
der Beobachtung zugänglich ist. 

Die Flächendichte der wahren Elektrizität auf einem 
Leiter, der an ein Dielektrikum angrenzt, ist durch 


(S4a) n,g"= 


e €@ 


gegeben, wenn g das Potential und n die Richtung der ins 
Innere weisenden Normale ist; also haben wir an der Anode 
z=0: 


Es ist daher zunächst die Feldverteilung zu berechnen; 
in der von uns behandelten Näherung finden wir aus (43) 
und (46): 
h =hy+ hg" 


(86) 


hy=*2 —n)de + h,*; 
0 


é 
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dabei ist h,* eine Funktion von t, die den zeitlichen Verlauf 
von h, für x= 0 darstellt. Sie bestimmt sich aus der Forde- 
rung (5) (man vgl. S. 231): 


l € 
4n 1 
87) =— =" fasfin 
0 


Denkt man sich die Lösungen p,, n, eingesetzt und die 
Quadraturen ausgeführt, so stellen (86) und (87) den zeitlich 
veränderlichen Feldverlauf beim Aufladungsvorgang dar. Wir 
begnügen uns damit, das Ergebnis der Rechnung für den 
stationären Endzustand [durch Einsetzen von (84), (84a)] an- 
zugeben. Es wird: 


(88) 


und damit haben wir — allerdings nur in erster Näherung — 
die Feldverteilung für den in $ 1 nicht zu Ende geführten 
Fall, daß die Konzentrationen der positiven und negativen 
Ionen fest vorgegeben sind und weder Bildung noch Wieder- 
vereinigung von Ionen stattfindet. 


Nun wenden wir uns wieder der Berechnung der zeitlich 
veränderlichen Influenzladung auf der Anode zu. Da h,= 7 
ist, ergibt sich aus (85), (86), (87): 


(89) o- fag 


_ Nun interessiert uns weniger diese Ladung selbst, als 
> Bee Änderung, die den beobachtbaren Strom I 
stellt: 


d "le d 
0 


Die Berechnung von (90) läßt sich vereinfachen, indem man 
die Differentiation nach der Zeit am Integranden ausführt 
und die Differentialgleichungen (44), (45) berücksichtigt; so 
findet man 


44) 
lin- 
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el 
Ae 
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I=1,+I1_, 
kD ö 
I,=—* _ Hp,)dz, 
(91) + l : ) 


Setzt man (83) in die Formel für Ry ein und führt noch die 
dimensionslosen Variablen 


(92) 


ein!), so ergibt sich das Schlußergebnis in der Form: 


+00 
V4T 


2ıT!+2m-1l+rT 2m—-1+vT 

2 A ver) 
93 20-142 T 
(93) —(2v-14+ 2m+vT)® Vit 


4 


_(2m+vT7? 
4? 


Die entsprechende Formel fiir J_ erhält man jetzt, indem man 
C,D, T,v mit C’, D’, T’, —v’ vertauscht und [wegen (91)] das 
Zeichen umdreht. 

Für kleine Werte von T, und damit von t, liefert nur 
der Summand m = 0 Beiträge, und zwar wird 


4 4) 

(99) I, = 1-5 Th — eT +...|, T<1. 
Der von den positiven Ionen herrührende Strom I, be- 

ginnt also mit dem Werte 


(95) 


CDr 


=Ckh,, 


l 


1) Wir unterscheiden v von dem im § 1 (S. 227 ff.) benutzten V’, 
weil dort die Voraussetzung (6) gemacht wurde. 
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wie nicht anders zu erwarten ist, da C die konstante Aus- 
gangsdichte der positiven Ladung ist. Betrachtet man an 
Stelle von I, den Quotienten rig ,/I%, so hängt dieser 
nur von der *Zeitvariablen T und dem’ Parameter v ab. Ent- 
sprechendes gilt für R = I_/I°. 

Man erkennt leicht, daB I, fiir groBe Werte von T ver- 
schwindet, denn solange T>2mv ist, verschwindet in (93) die 
geschweifte Klammer (wegen lim (x) = 1); ist aber T ver- 


gleichbar mit 2mv, so ist m negativ und groß, so daß der 
Faktor e”” den betreffenden Summanden zu null macht. Da- 
gegen erweist es sich als nicht durchführbar, unmittelbar aus 
(93) eine asymptotische Formel für große Werte von T ab- 
zuleiten. Wir greifen darum auf die Differentialgleichung (44) 
zurück und behandeln sie nunmehr nach der Methode der 
Partikularlösungen. 


In bekannter Weise findet man, dab 
He 


tt 
Pi, =e [an cos "7" 


partikuläre Integrale der Gl. (44) sind, wenn 
en 


(97) iD += ‚also + T 


gesetzt wird. Die p,,, geniigen auch noch den Randbedingungen (4) 


und werden dadurch zu Eigenfunktionen unseres Problems, wenn die 
Bedingung 


(98) Bn 


+ 


v 
2ma 


erfüllt ist. Das System (96) ist aber nicht vollständig, sondern es kommt 


noch die stationäre Lösung 
(99) Pro = Ao 


hinzu. Die Superposition der Lösungen (96) und (99) 
oo 
(100) p, = A,e”* + [4 cos + Bu sin 


läßt sich nun folgendermaßen einem vorgegebenen Anfangszustand (62) 
anpassen. Integriert man zunächst 


(101) f(z) = + 3 [4- cos T + Bn sin” 0sısl, 
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nach x von 0 bis 1, so verschwinden wegen (98) alle Glieder bis auf 
das erste, und es ergibt sich 


l 
Hf dx 
0 
(102) A, = 


Danach wird 
He H x 


(103) F(a) = Aye 


oo 
- > eo Bu sin ™ 
1 


Um hieraus die Koeffizienten A„, B,, in Übereinstimmung mit (98) zu 
bestimmen, zerlegen wir F(x) in zwei Teile 


(104) F(x) = g(x) + u(@), 


deren einer g(r) als gerade Funktion und deren zweiter u(r) als un- 
gerade Funktion über 0 hinaus von x = 0 bis « = — 1 fortgesetzt werden 
soll. Daher muß werden 


co 
(105) g(x) = >in A, €08 ‚u@)= >" B,, sin = 
1 1 


und wegen (98) 
H 
(106) w (a) = (x). 


Wir erhalten also fiir u(x) eine lineare Differentialgleichung mit 
konstanten Koeffizienten 


H H 


die stets integriert werden kann. Die dabei auftretende Konstante ist 
so zu bestimmen, daß 


(108) 
0 


(107) (a) + 


wird, denn die Fourierentwicklung von g (x) darf nach (105) kein Absolut- 
glied enthalten. 

Damit ist unser Problem vollständig gelöst, denn sobald g (x) und 
u (x) gemäß (104), (107), (108) bestimmt sind, ergeben sich die A, und B„ 
als die Fourierkoeffizienten der Funktionen g (x) bzw. u (a). 

Wir führen die Rechnung wieder für die konstante Anfangs- 
verteilung (76) durch. Es wird, mit Hl=v, 


(109) 
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und damit erhält man für die Koeffizienten in (100) 


Cr 
q'*(e” 1) 
(110) 2mn 2Cv 1- (Dre 2 
v qm? n® 1+ v 
Zmn 


Geht man mit dieser Lösung wieder in (91) hinein, so 
ergibt sich für den Aufladestrom I, das Schlußresultat 


m v 


| \sis 
| 

Diese Formel ist mit (93) völlig gleichwertig. Beide 
gelten für alle Werte der Parameter, doch sind sie für sehr 
verschiedene Parameterbereiche bequem anwendbar und er- 
gänzen sich in dieser Hinsicht. Man erkennt ohne weiteres, 
daß für große Werte von T (und zwar schon von T—1 an) 
die Formel (111) sich auf das erste Summenglied reduziert; 
der Abfall erfolgt also dann exponentiell und zwar für alle 
Werte von v. 

Die numerische Rechnung zeigt, daß Formel (111) für 
Parameterwerte, die kleiner als etwa v = 10 sind bis zu etwa 
T~0,1 herab sehr bequem anwendbar ist. Dagegen erfordert 
sie für noch kleinere Werte von 7, namentlich bei großen 
v-Werten, die Berücksichtigung von vielen Gliedern. Gerade 
unter diesen Umständen führt aber die Formel (93) — bzw. 
die daraus abgeleitete Formel (94) — schnell zum Ziel. In 
der oberen Hälfte der Fig. 3 sind die Werte von R= 1 für 
einige Parameterwerte als Funktion von T = = gezeichnet. 
Man erkennt, daß die Kurven sich für kleine v einer Grenz- 
lage nähern. In der Tat folgt aus (111) und (95) 


8 = g Cater 
(112) lim R= 


(ll) I, = 


oo vw 


diese Kurve hat selbst keine physikalische Bedeutung, da fiir 
v= 0 keine Ionenbewegung mehr stattfindet, aber sie gibt für 
0 <1 die Verhältnisse hinreichend genau wieder; bei v = 0,5 
beträgt die Abweichung nur noch etwa 0,5°/,. 

16* 
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Für große Werte von v und 1/7 läßt sich die weit. 
läufige Formel (93) erheblich vereinfachen, indem man die 


1 


95 


025 N 


U=4 


v:0 
v:2 


-025 


0% 


Fig. 3. Der relative Ladungsstrom R u. der relative Entladungsstrom R 
als Funktionen der Zeitvariablen 7 für verschiedene Werte von v 


asymptotischen Formeln für ®(z) heranzieht. Unter der 
Annahme 


(1133) 


und der weiteren Einschränkung, daß vT weder nahe an 1 
noch nahe an 0 liegt, so daß auch noch 


aus, 
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= 2 
(113b) >1 wud 


ist, kann man sich auf das Summenglied m = 0 beschränken 
und findet dann 


( _ ten? 
CDv 2 y= +7 
(114) om 


Unter den Annahmen (113a) und (113b) ist also der Verlauf 
im wesentlichen geradlinig 


(14a) 1, = - oT] 


und die ganze I,-Kurve setzt sich dann näherungsweise aus 
3 Stücken zusammen: einem parabolischen Stück (94) in der 
Nachbarschaft von T = 0, einem exponentiellen Stück für 
große T und einem geradlinigen Stück zwischen beiden. Das 
geradlinige Stück verläuft um so steiler, je größer v ist. 

Der Verlauf der R’-Kurven wird ebenfalls durch Fig. 3 
wiedergegeben, wenn man T,v durch T’,v’ ersetzt denkt.') 
Der von den negativen Ionen herrührende Strom J_ [vgl. (91)] 
hat natürlich das gleiche Vorzeichen wie I a) Sind die Kon- 
stanten der beiden Ionensorten verschieden, so dürfen die 
Kurven für 7, und I_ natürlich erst dann superponiert werden, 
wenn sie auf eine gemeinsame Abszisse umgerechnet sind. 
Dadurch kann eine groBe Mannigfaltigkeit von Kurvenformen 
entstehen. Sind aber die Konstanten k und D für beide 
lonensorten gleich, so gilt die Darstellung Fig. 3 auch für 
(dd, +1) 


+ 1°) mit den gemeinsamen 


den relativen Gesamtstrom 


Werten von v und T. 

In derselben Weise, in der wir hier den Anfangszustand 
p=C behandelt haben, kann man mit Hilfe der allgemeinen 
Formeln (64) und (75) bzw. (100) auch den Fall erledigen, daß 
als Anfangszustand diejenige Ionenverteilung vorgegeben wird, 
die sich unter der Wirkung eines Potentials V, ausgebildet 
hat, und daß nunmehr ein davon verschiedenes Potential V, 


1) Denn die > in (111) ist eine gerade Funktion von v. 
2) Denn führt man in (111) die auf S. 240 gegebene Substitution 
aus, so zeigt sich, daß I_ das gleiche Zeichen hat, wie /,. 
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angelegt wird. Die bekannten Anomalien, daß bei Erhöhung 
des Potentials wiederum ein zeitlich abnehmender Strom gleicher 
Richtung fließt, bei Erniedrigung des Potentials aber zunächst 
ein Strom entgegengesetzter Richtung auftritt, würden dann 
zwanglos ihre Deutung finden. Das läßt sich ohne alle Rech- 
nung voraussehen. Nach (85) ist nämlich die gesamte Influenz- 
ladung, die sich auf der Anode ausgebildet hat und die also 
als zeitliches Stromintegral tg Idt beobachtet werden kann, bis 
auf den Faklor e/4r durch den Wert der Feldstärke h,., 
gegeben. Dieser Wert ist nun aus (31) bzw. (31’) berechenbar 
und wächst natürlich mit V (man vgl. Fig. 1), daher bewirkt 
Erhöhung von V ein positives, Erniedrigung von V ein nega- 
tives Stromintegral von angebbarem Betrage. 

Wir beschränken uns darauf, den einfachen Fall zu be- 
handeln, daß (84) bzw. (84a) als Anfangszustand gegeben sind 
und nunmehr zur Zeit {= 0 das Feld abgeschaltet wird. In 
der von uns behandelten ersten Näherung ist dann h, =0, 
die Differentialgleichungen (44), (45) reduzieren sich auf die 


Wärmeleitungsgleichung und die Grenzbedingungen auf 2 = 0, 


Sm = 0 fir In bekannter Weise findet man 


oo 
1 2C (— mae 
(115) =C + >} aay e 7008 


v 


nebst einer entsprechenden Gleichung für n,, die durch Ver- 
tauschung von C,T,v mit (”,T’, —v’ entsteht. An die Stelle 


der Gil. (91) tritt jetzt für den Entladungsstrom I: 
I= + ) 
(116) (p, — 9,0), 


(4, — 2, 0); 


daraus ergibt sich 


co 
o it 
mit einer entsprechenden Formel für /_, die man erhält, wenn 
man C,D,v durch C’, D’,(—v) ersetzt und das Zeichen um- 
dreht. 
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Für t=0 ergibt sich 


oo 
T%__4CD 


1+ 
vr 
die Summation läßt sich mit Hilfe der Formel 
l—-e” 
ausführen'), und man findet 
(120) Ty =— ckh,. 


Dieser Wert ist entgegengesetzt gleich dem Anfangswert (95). 
Läßt man also ein Potential V = h,l solange wirken, bis sich 
der stationäre Zustand ausgebildet hat und schaltet dann das 
Potential ab, so setzt der Entladestrom in entgegengesetzter 
Riehtung mit dem gleichen Wert ein, mit dem der Ladestrom 
begonnen hatte, falls die Ionen gleichmäßig verteilt waren. 


Für v = 0 erhält man aus (117) fiir R = a den Wert (112) 


+ 
mit entgegengesetztem Zeichen; fiir kleine Parameterwerte 
ist also die Entladungskurve ein genaues Spiegelbild der 
Ladungskurve. Für Werte v>1 besteht diese Beziehung 
nicht mehr, wie aus Fig. 3 hervorgeht, in deren unterer Hälfte R 
für die gleichen Parameterwerte gezeichnet ist, wie in der 
oberen Hälfte der Verlauf von R. Man erkennt daraus, daß 
zwar anfangs der Entladungsstrom schneller abfällt, als der 
Ladestrom, daß aber dennoch — namentlich bei größeren 
v-Werten — der Ladungsvorgang schneller beendigt ist, als 
der Entladungsvorgang. 

Es sind noch einige Worte über den Gültigkeitsbereich 
unserer ersten Näherung (44) bis (46) zu sagen. Wir sind 
zu ihr gelangt unter der zweifachen Annahme, daß erstens 
Bildung und Rekombination vernachlässigt werden können und 
daß zweitens das Feld, in dem die Ionen sich bewegen, als 
homogen angesehen werden darf. Die erstere Annahme ist 
dann gerechtfertigt, wenn die Zahl der positiven und negativen 


1) Man erhält diese Formel, indem man die Anfangsverteilung 
v et zT 


y (x) = 


| 
in eine Fourier-Reihe entwickelt und dann y (l) — y (0) bildet. 
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Ionen durch irgendwelche andere physikalische Bedingungen 
vorgeschrieben ist, man denke etwa an die Verhältnisse bei 
der Lösung eines sehr verdünnten Elektrolyten. Aber auch 
wenn Bildung und Wiedervereinigung stattfinden, tritt ihr 
Einfluß gegenüber der Diffusion um so stärker zurück, je kleiner 
der Plattenabstand ist.) Auf der anderen Seite ist auch im 
stationären Endzustand die Inhomogenität des Feldes um so 
geringer, je kleiner die reduzierte Dicke d, und demnach der 
Plattenabstand 1, ist. (Man vgl. 8. 230/231.) Unsere Lösung wird 
also eine um so bessere Approximation an die strenge Be- 
handlung darstellen, je kleiner der Plattenabstand ist; dann 
sogar auch noch für beliebig lange Zeiten t. 

Ist andererseits J groß, so kann unsere Lösung erster 
Näherung nicht mehr für beliebig lange Werte von t korrekt 
sein, weil dann das Feld im Endzustand sehr unhomogen ist 
und diese Inhomogenität sich im Laufe der Zeit ausbildet. 
Es läßt sich auch voraussehen, in welchem Sinne die Ab- 
weichungen unserer Formeln liegen werden. Da die Ausbildung 
der Inhomogenität des Feldes der Bewegung der Ionen ent- 
gegenwirkt, wird bei der strengen Lösung der Strom I schneller 
absinken, als nach der ersten Näherung. 

Für große Werte von ö kann also von unseren Formeln 
nicht mehr erwartet werden, als daß sie den Verlauf — zwar 
anfangs richtig — aber weiterhin nur qualitativ wiedergeben. 
Dann müssen also die weiteren Näherungen (43) herangezogen 
werden. Die Differentialgleichungen für die höheren Nähe- 
rungen p,,n, werden alle von der Form (44), (45), nur daß die 
Gleichungen nicht mehr homogen sind, sondern daß auf den 
rechten Seiten vorhergehend bestimmte Funktionen von x und t 
stehen; das zugehörige h, berechnet sich jedesmal hinterher 
aus einer mit (46) gleichlautenden Gleichung. Wir sehen davon 
ab, an die Behandlung der zweiten Näherung heranzugehen, 
solange so wenig Daten, wie bisher, bekannt sind, die eine 
quantitative Prüfung der Theorie ermöglichen würden. Mathe- 
matisch begegnet die Behandlung der höheren Näherungen 
keinen wesentlichen Schwierigkeiten; mit Hilfe der Eigen- 
funktionen (96), (99) läßt sie sich nach der Methode durch- 
führen, die aus der Wellenmechanik hinreichend bekannt ist. 


1) P. Langevin, Journ. de Phys. (4) 4. S. 322. 1905. 
(Fortsetzung folgt.) 


(Eingegangen 21. Oktober 1932) 
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